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HAUPTAUFSÄTZE 


Berechnung der Reibung und Tragkraft eines endlich breiten 
Gleitschuhes auf ebener Gleitbahn. 


Von W. Frössel in Göttingen. 


Die vorliegende Arbeit behandelt den dreidimensionalen Strömungsvorgang, 
wie er zwischen einem endlich breiten ebenen oder balligen Gleitschuh und 
einer ebenen Gleitbahn auftritt, wenn bei guter Schmierung zwischen beiden 
eine relative (Geschwindigkeit herrscht. In der weiteren Auswertung der 
entwickelten Theorie wurden insbesondere die Trag- und Reibungskräfte für 
die genannten Gleitsysteme berechnet. Die Ergebnisse sind zur Erleichterung 
der Anwendung in Kurven dargestellt. 


I. Einleitung. 


Das Ziel der hydrodynamischen Schmierungstheorie ist die vollständige Beherrschung 
der Strömungszustände, insbesondere der Trag- und Reibungszustände im Schmierspalt eines 
Gleitlagers oder einer Führung mit endlicher Breite bei gegebenen Lager-, Lauf- und Schmie- 
rungsgrößen. Dieses bedingt die mathematische Behandlung des dreidimensionalen Strömungs- 
vorganges im Spaltgebiet, die mit einer zweidimensionalen Druckv erteilung verbunden ist, wie 
sie von Michell') für einen schmalen ebenen Gleitschuh (Breite < Länge), von Vogelpohl?) 
für den allgemeinen Gleitschuh nach dem Prinzip des Minimalsatzes der Lagerreibung und 
von Nahme?°) für den in der Breite einseitig begrenzten balligen Gleitschuh durchgeführt ist. 


Die vorliegende Arbeit gibt eine weitere Lösung des genannten Problems an, die ins- 
besondere für Gleitschuhe mit Breite = Länge konvergiert. Ihre Anwendung in der Praxis ist 
durch die weitgehende zahlenmäßige Auswertung der Funktionen und Darstellung in Zahlen- 
tafeln und Diagrammen wesentlich erleichtert worden. Aber auch hier beschränkt sich wie 
in den schon erwähnten Arbeiten die theoretische Behandlung auf den Gleitschuh auf ebener 
Gleitbahn. 


II. Entwicklung der Formeln. 


a) Allgemeines. Der Strömungsv organg der dreidimensionalen Spaltströmung zwischen 
Gleitschuh und Gleitbahn erfordert zu seiner Berechnung gegenüber der zweidimensionalen 
eine bedeutende Erweiterung der theoretischen Behandlung. In erster Linie macht sich das 

1) A. G. Michell: Z. Math. Phys. Bd. 52 (1905), S. 123. — Ostwalds Klassiker Nr. 218, Leipzig 1927, S. 202. 


2) G. Vogelpohl: VDI-Forsch.-Heft 386, Berlin 1937. 
3) R. Nahme: Ing.-Arch. Bd. 11 (1940), S. 191. 
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durch die endliche Begrenzung des Gleitschuhes in der Breite verursachte seitliche Abfließen 
des Schmiermittels, verbunden mit einem seitlichen Druckabfall, geltend, dessen Auswirkung 
sich auf die Tragfähigkeit und Reibung des Gleitschuhes fortpflanzt und durch einen erhöhten 
Arbeitsaufwand ausgeglichen werden muß. 

Wie in fast allen Strömungsfragen kommen auch hier zuerst die Navier-Stokesschen 
Bewegungsgleichungen für Flüssigkeiten in Betracht, die auch für veränderliche Zähigkeit in 
Vogelpohl?) Gl. I ausführlich aufgeführt sind. Um aber bestimmte Lösungen zu erhalten, 
muß ihre große Vielseitigkeit zunächst durch Vernachlässigung unwichtiger oder sehr kleiner 
Glieder stark eingeengt werden. Dabei scheiden die Trägheitsglieder infolge der überwiegen- 
den Zähigkeit des Schmiermittels und der im Spalt herrschenden laminaren und stationären 
Strömung aus. Eine weitere Vereinfachung der Bewegungsgleichungen wird durch die An- 
nahme einer im ganzen Spaltgebiet konstanten Schmiermittelzähigkeit erreicht. 

Werden die Koordinaten x, y, z so gelegt, daß die x-Achse mit der Bewegungsrichtung 
der Gleitbahn, die y-Achse mit der Spalthöhe und die 2-Achse mit der Gleitschuhbreite zu- 
sammenfällt, so darf die mittlere der entsprechenden Geschwindigkeitskomponenten u, v, w 
eines Flüssigkeitselementes als verschwindend klein und die auftretenden Drücke p in dieser 
y-Richtung als unveränderlich angesehen werden. Unter Berücksichtigung dieser Annahmen 
vereinfachen und verringern sich die Navier-Stokesschen Gleichungen auf die beiden nach- 
stehenden 
op "u du 0 

N | —— 2 


+33 #55 
dx 0a? 099 02 
2 BE re SE 
dp mw %w  O’w 
lt at) 0 
02 oe! Oy’ dr. 


die aber wegen der schwierigen Klammerausdrücke einer Lösung noch wenig zugänglich sind 
und deshalb auf weitere Vereinfachungen untersucht werden müssen. Letzteres geschieht 
durch Abschätzen der Größenordnungen der einzelnen Klanımerglieder, wozu für y die Größen- 
ordnung des Spaltes h, für x die der Gleitschuhlänge ! und für z die der Gleitschuhbreite b 
gesetzt wird, wobei b und ! groß gegen h sind. Der Vergleich der zweiten Ableitungen 
0” 1 0? 1 0° 6, re \ um 
le” Erde” r und Fre = zeigt, daß die nach «x und z abgeleiteten Glieder gegen das 
nach y abgeleitete verschwindend klein und daher vernachlässigbar sind. Die Gl. (la) gehen 
somit über in 

0P u 


day 
‚(db, 
%p_ "wm 0 
d2  NdyT 
die mit der Kontinuitätsgleichung 
du om 
=> ı 
EM ne nn re 
a . R dp dp 
und den Randbedingungen y=0, u=U; y=h, w=0 unter Konstanthaltung von Sg und zZ 
für veränderliches y% auf die allgemeine Differentialgleichung dieses Strömungsvorganges 
d° op, Bdhdop 6nTdh 
ce Pı Re N N 
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führen. Hierin bedeutet h die Spalthöhe, deren Veränderung durch die Gestalt des Gleit- 
schuhes und der Gleitbahn vorgegeben ist und die nur eine Funktion von x sein soll, U die 
Geschwindigkeit der Gleitbahn und 7 die Zähigkeit des Schmiermittels. Zur Beschränkung 
auf die für die praktische Anwendung am wichtigsten Spaltverläufe sind für die weitere 
mathematische Entwicklung der ebene und der ballige rechteckige Schmierkeil ausgewählt 
worden. 

b) Der ebene Schmierkeil*). Der ebene Schmierkeil wird durch die innere Be- 
grenzung des ebenen Gleitschuhes und der ebenen Gleitbahn gebildet (Bild 1). Die unend- 


#4) O. Reynolds: Phil. Trans. Roy. Soc. .Lond. (A) Bd. 177 (1886), S. 157. Ostwalds Klassiker Nr. 218. 

*) Nach Fertigstellung dieser Arbeit wurde mir eine Abhandlung von M. Muskat, F. Morgan und M. W.Meres 
unter dem Titel „The Lubrieation of Plane Sliders of Finite Width“ erschienen im Journal Applied Physies, Bd. 11 
(1940), Nr. 3, S. 208 bis 219, bekannt, die dasselbe Problem des ebenen Schmierkeils mit gleichem Lösungsansatz in anderer 
Sehreibfornm enthält. Trotzdem erscheint es nützlich, mit Rücksicht auf die späteren zahlenmäßigen Auswertungen den 
etwas anders gearteten Entwicklungsgang gekürzt wiederzugeben. 
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lich breit und lang gedachte Gleitbahn bewegt sich in der angegebenen Pfeilrichtung mit der 
Geschwindigkeit U unter dem örtlich festgehaltenen Gleitschuh entlang und reißt dabei das 
vor dem Gleitschuh befindliche Schmiermittel 
mit in den Schmierspalt hinein. Dabei ent- y 
stehen unter Mitwirkung der in der Bewegungs- 
richtung veränderlichen Spaltform die im fol- 
genden gerechneten Trag- und Reibungsgrößen. 
Bild 1 zeigt einen Schnitt durch die Längsachse 
des Gleitsystems. Die Spaltbeziehung für den 
ebenen Schmierkeil lautet bei den kleinen vor- 
kommenden Winkeln ö 






Gleitbahn 





khe—2)0 . .. . 


worin die Länge a vom Gleitschuhanfang bis Pr 

zum Schnitt der beiden verlängerten Gleit- Bild. Schematische Darstellung des ebenen Gleitsystems 
flächen eine für die Gleitschuhstellung charak- 

teristische Größe darstellt und x der Abstand vom Spaitanfang im Sinne der Gleitbahn- 
bewegung bedeutet. Für beide Werte ist die Gleitschuhlänge als Maßeinheit zu nehmen (vgl. 
Bild 1). Mit diesem Verlauf erhält die Differentialgleichung (3) die nachstehende Gestalt 


0’p 0’p 3 dp 69 U 
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Als Lösung wird dafür eine Funktion in der Form 


es 
Cos 4, 2\ . \”/Cosi„_,z Cos4,z\ { 
er Cte+ Cos 4, GB Cosi.p)fn(®) er 
==] 
oder ausgeschrieben 
Cos 4, 2 Cos4,2 Cosi, z 
p Cos 4,b fota) + Cos4,b LCos4,b f@) 


! : (6a) 
‚[jCos4,2 Cos4,2\. 
| - \eta)+:-: 


Cos4,b Cos4,b 


) 
angesetzt, die als neue Größen die halbe dimensionslose Gleitschuhbreite b 5. die noch 
zu bestimmenden Werte A; und die noch unbekannten Funktionen f; enthält. Während von 
den vier Grenzbedingungen, die für die Spaltränder bei z=+b und <=0 und 1 den gleichen 
Spaltdruck p= Außendruck = Null vorschreiben, die beiden der Seiteiikanten bei z=+b 
durch die rechte Seite der Gl. (6) erfüllt sind, müssen die beiden übrigen für die Vorder- 
und Hinterkante bei e=0 und 1 durch richtige Wahl der f;-Funktionen nach befriedigt 
werden. 


Werden nun die ersten und zweiten Ableitungen nach x und z von Gl. (6) in die 
Differentialgleichung (5) eingesetzt, so wird diese erfüllt, wenn auch ihre linke Seite ver- 
schwindet, d. h. wenn jede der daraus entnommenen, nachstehend aufgeführten Ausdrücke 
verschwindet. Diese sind gleichzeitig die Differentialgleichungen für die f;- Funktionen und 
heißen somit 
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Die erste dieser Gleichungen ist die Differentialgleichung für den unendlich breiten, ebenen 
Schmierkeil mit zweidimensionaler Strömung. Sie hat die Lösung 


i 6nT zil—i) (8) 
— Ts o 9). 
I %(2a—1) (a— x) 
j - a R Co84„_,2, . . 
Bei den übrigen Gleichungen ist der Faktor Cosi „n= 1,2,...) stets verschieden von 
ıN—1 
Null, so daß für das geforderte Verschwinden der linken Seiten die Klammerausdrücke sorgen 
müssen. Diese lassen durch die gleichgebaute Form der ersten und zweiten Hälfte sofort 
den Wert f„=f„_, als partikuläre Lösung erkennen. Für die allgemeine Lösung wird nun 
die Beziehung 


h=fn 1 IAn-ı 


angesetzt, die nach Einsetzen der entsprechenden Ableitungen die Klammerausdrücke der 
Gl. (7) auf die einfachere Form 


[7 ’ PR} Be ( 
n-ı "3 _ ; 9 en ie Vi 


bringt, die die Lösung 


Au. Fi An (0 — x) + B, 1 N, An ‚(a x)} 
a x 
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hat (n=1,2,...), worin A„_,, B„-, und /„_, noch verfügbare Konstanten sind und J, und N, 
aus bekannten Tabellen’) für Zylinderfunktionen (Besselsche Funktionen) entnommen werden 
können. Wie die Funktionen f müssen auch die Funktionen g (für <=0 und 1) verschwinden. 
Das bedeutet homogene Randbedingungen, die für die Lösungen einer homogenen Differential- 
gleichung wie Gl. (9) bekanntlich nur Eigenfunktionen zulassen. In solchen Fällen erfährt 
die numerische Auswertung eine wesentliche Erleichterung, wenn die allgemeine Lösung der 
Gl. (10) durch eine ihrer Konstanten dividiert wird. Sie heißt dann 


_ In-ı _Jı An-, (a — + Cn-ı Nıtin-ı (a —8)} 


G = (il). 
re 7 a—ıa 
C„-, sowie 4„_, lassen sich allein mit den schon angegebenen Grenzbedingungen @,„_,—=0 
für e=0 und 1 berechnen. Ihre Bestimmungsgleichungen lauten dann 
I) x 7 N 
für -—=0 J YA 45 + En-ı N, tin-ı 0} Br (12a) 
. a . . . . . . . . . ic , 
für =1 J,; \An- ‚(aD + Cu, N tina, =0 (12b) 
. m 1 . . . . . — . 
Die Ermittlung der Eigenwerte /„_, folgt dann aus 
FG a) 
„ / ı An - rm F 
J {i,_,(a—-1)} -— NG, -WM=0 Ww=1%..)e. . . . (iM. 


f 64...) 


Für die Anwendung der A-Werte wird der kleinste positive Eigenwert mit /„_,=4, und die 
weiteren in positiver Richtung aufeinander folgenden mit den laufenden Zahlen 1,2,3... als 
Indizes bezeichnet. Diese so gefundenen Eigenwerte bestimmen zusammen mit den Kon- 
stanten C„_, die Grundformen der Kurven von Gl. (11) und ihre Schnittpunkte mit der x-Achse, 
wobei die Anzahl der Schnittpunkte in dem betrachteten Bereich O< x<S1 durch die Ordnungs- 
zahl n— 1 der Gl. (11) angegeben wird. 


5) Jahnke-Emde: Funktionstafeln mit Formeln und Kurven. 3. Aufl. Leipzig und Berlin 1938. Watson: 
A Treatise on the theory of Bessel funetions, Cambridge 1922. 
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Um die in Gl. (6) stillschweigend vorausgesetzte Konvergenz der Glieder zu er- 
reichen, müssen die höheren Funktionswerte 




















i h REG a j 
von fa =fn-ı — In, In dem ganzen betrachte- 9 j 
ten Bereich O0<Sr<s1 gegen Null gehen, ru nn IN 
d. h. die 9-Funktionen müssen sich in dm „h- _ i 
genannten Gebiet möglichst gut den f- Funk- h 
tionen anpassen, damit ihre Differenzen klein 20 


werden. Es gilt also mit dem in Gl. (11) noch ’ 
verfügbaren Faktor A„_, die in der Kurven- . 
grundform @,„_, bereits angegebenen Funk- 
tionswerte derartig zu verändern, daß die 
folgenden aus f„_, — 4n-., neu entstehenden 
f-Funktionen stets bessere Anpassungsfähig- 
keiten für die folgenden g-Funktionen bieten. 
Bei der Auswertung dieser Methode zeigt es 
sich, daß zum Anpassen an jede weitere f- 
Funktion immer gerade die mit dem nächst- 
höheren 4-Wert gebildete g-Funktion geeignet 
ist. Daraus ergibt sich durch Betrachtung 
von Gl. (6), daß ihre Konvergenz durch die 
nach den höheren Gliedern steigenden A-Werte 
in den Klammerausdrücken bedeutend ge- 
steigert wird. Das Bild 2 zeigt als Beispiel # 
das Annähern der höheren Funktionswerte f; 
an Null für die Gleitschuhlage a—=1,25. Für * 
andere häufig vorkommende Gleitschuhlagen 
sind die Funktionswerte f; sowie die zuge- 
hörigen A;-Werte in gleicher Weise ermittelt gzl- 
und in der Zahlentafel 1 zusammengestellt. 
Aus den in der Zahlentafel 1 dargestellten 47 
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Größen ist der gemeinsame Faktor I schon „° 
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herausgezogen worden. Mit dieser Maßnahme 


. . Bild 2. Verlauf der Funktionen f; und 9; beim ebenen 
lautet nun die Druckgleichung i ” 2. allein Saure: 


Gleitsystem (a =1,25). 
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mit der die Druckverteilung für jede gegebene Gleitschuhstellung a (vgl. Bild 1) über die 
ganze Gleitschuhfläche errechnet werden kann. Die Rechnung gestaltet sich für die in der 
Zahlentafel 1 aufgeführten Werte von a besonders einfach, da hierfür die benötigten f;- und 
/-Werte unmittelbar daraus entnehmbar und in die Gl. (14) einzusetzen sind. Zu beachten 
ist, daß für die Größen z und b die Gleitschuhlänge als Maßeinheit genommen werden muß. 

Die Gesamttragkraft P des ebenen Gleitschuhes ist durch graphische Integration der 
Gl. (14) nach x und z erhältlich. Sie liefert bei Division mit der tragenden Gleitschuhfläche, 
die wegen der Wahl der Gleitschuhlänge als Maßeinheit zahlenmäßig durch die Gleitschuh- 
breite B=2b ausgedrückt wird, den mittleren Druck oder besser die Tragkraft pro Breiten- 
einheit in dimensionsloser Form 


b 


1 
pe | e Ne )Jarae—m, Be EAN 


06 


die aus Gleichgewichtsgründen gleich der zu tragenden Last pro Breiteneinheit sein muß. 
Aus der Division der Gl. (14) und (15) folgt bei gleichzeitiger Beseitigung der oftmals un- 
bequemen Größen n, U und ö die dimensionslose Druckverteilung p/p,, die für jede Gleit- 
schuhstellung a eine andere Gestalt besitzt, sonst aber von n, U und ö unabhängig ist. Sie 
stellt also die Druckverteilung im Schmierspalt für den Belastungsdruck p,=1 dar und ergibt 
in einfachster Weise durch Multiplizieren mit jedem anderen gewünschten Belastungsdruck p, 
die entsprechende Druckverteilung. 

Auch der für die punktförmige Anbringung der Last wichtige Angriffspunkt I, (Bild 1) 
läßt sich für jede Gleitschuhstellung a aus der zugehörigen Druckverteilung durch Ermittlung 
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des Druckpunktes, dessen Lage aus Gleichgewichtsgründen mit der des Angriffspunktes der 
Last zusammenfallen muß, nach dem Momentensatze zu 
1b 


1 6nVÜ r 
ve her] \dirdz. N . (16) 


Tl) 


bestimmen, worin m, aus Gl. (15) zu entnehmen ist. 


Die zwischen einem Gleitschuh und einer Gleitbahn auftretende Reibung wird aus der 
Schubspannungsverteilung an der Gleitbahn nach dem bekannten Ansatz 


U dph‘) 


nn de? (19) 


T= 


erhalten. Die Umformung mit der Spaltbeziehung des ebenen Gleitschuhes G]. (4) lautet dann 


U dp(a— .r)ö 


> 


(a x)o dr > ER: 


T- Y) 


. dp... u nn : 
Den noch unbekannten Druckgradienten d £ liefert die Differentiation der Druckgleichung (14) mit 
dp 6nU \ Cos 4, 2 Cos4,2 Cosi, z 


(n \ 
de | Na+ Hr... . (8, 


Cos/,b Cos ),b) 


If) + 


Cos4,b 
i a i me h B ö 6nT a+r—2axr 
worin die zugehörigen f}-Werte mit der differenzierten Gl. (8) fj, = = "Ba-T)Ka ); er- 
. ö” 2a— — x)‘ 

rechnet und in der Zahlentafel 1 zu finden sind. Ebenso wurden die übrigen f};-Werte ver- 
mittels f,=f},-ı — 9), ,, die die Differentiation der früher angewandten Beziehung f, = f„_, — In-ı 
darstellt, gefunden und eingetragen. Dazu war g,_, mit den bekannten Differentialformeln 
der Besselschen Funktionen Gl. (10) zu bestimmen, die für den vorliegenden Fall umgeformt 
lauten 





2 v; AT : 
‚ a 2J, in ‚ (a 7), An-ı 4; Al 
In-ı = An-ı (a r)? m „Jo \An ‚a — x)s| 
2X, {d„-,‚(a—a) 4 
N, \An-ı aur er F 
+ B.-: Fr en N, ti„_, (a — x)\ 





Mit der Abkürzung F(x,2) für den geschweiften Klammerausdruck in Gl. (18) heißt diese 


6nU ’ u. ; 4 ö 
“: — 1 - F(x,2). Die Gl. (17a) schreibt sich damit 


de 6 A 
4 nÜU | 1 er Fa ; | ru 
T > 3: F (ıx,2)- (a #), (19a) 
und nach Elimination von ö mit Gl. (15) 
1 
+3-F(x,2)-(a--.x) 
T a— x 
= - . (19b). 
YnUpr v6m, 


Die mittlere Schubspannung r,, die wegen der eingeführten neuen Maßeinheit !=1 (Bild }) 
gleichzeitig die Reibungskraft pro Breiteneinheit darstellt, wird nun durch Integration der 
Gl. (19) und Division mit der Gleitschuhfläche (zahlenmäßig gleich der Gleitschuhbreite B= 2b) 
erhalten zu 


ıb 
GG < WB 2 Wr Mar 90: 
"1.:7 7 \ \ rt (2,2). (a ef A — 
00 
oder 


2 ER ‘ | ) 
— ® x.2)-( x g ER e  —  \ 
VnUp, by6 1 \ la „+3F (x, 2): (a eo); dxd (20b) 
00 


Die rechte Seite von Gl. (20a) und (20b) enthält nur Abmessungen des Gleitschuhes und die 
Gleitschuhstellung a. Sie ist daher für einen gegebenen Gleitschuh mit vorgeschriebener 
Stellung a konstant. Die Veränderung der mittleren Schubspannung r, hängt somit von der 
Zähigkeit des Schmiermittels n7, der Gleitbahngeschwindigkeit 7, dem Winkel ö oder dem 
äußeren Belastungsdruck p, ab. 
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ec) Der ballige Schmierkeil. Der ballige 
Schmierkeil unterscheidet sich vom ebenen durch die 
gekrümmte Gleitfläche des Gleitschuhes (Bild 3). Die 
theoretische Entwicklung des hydrodynamischen 
Schmierungsvorganges folgt im Prinzip der früheren. 
Trotzdem erscheint es angebracht, wegen des anders- 
gearteten Charakters des Spaltverlaufes die Herleitung 
der benötigten Funktionen durchzuführen und zalılen- 
mäßig für die spätere Verwendung vorzubereiten. 
Bild 3 zeigt einen Schnitt durch die Längsachse des Bild 3. Schematische Darstellung des balligen 
Gleitsystems. Die Arbeitsweise ist wie beim ebenen Gleitsystems. 
Gleitsystem. 

Die charakteristische Größe dieses Schmierkeils ist der engste Spalt h,, Bild 3, dessen 
Lage durch die vom Mittelpunkt des Gleitschuhkreises aus auf die Gleitbahn gefällte Normale 
gegeben ist. Der Schnittpunkt dieser Normalen mit der Gleitbahn wird als Ursprung des 
x, y,2-Koordinatensystems gewählt. Als Spalt kann der senkrecht zur Gleitbahn gemessene 
Abstand zwischen Gleitbahn und Gleitschuh genommen werden, solange x,<r’ ist (vgl. 
Bild 3). Die Spaltbeziehung lautet dann R=h,+ r’ — Yr’?— x? (21) und bei Vernachlässigung 
der höheren Glieder einer Reihenentwicklung 





z 
= h, .- 5 > De Ba A Bud AB TER SE 2 Se En SR BE Zen (22). 


Mit diesem Spaltverlauf erhält die allgemeine Differentialgleichung (3) die Form 
0’p 0’p | 6x op rs 


; ; x Ba 1.2732 1 
02? 0° ( ». dx F 
1 +3 - 77 Je r'h, (1 4 Ip m h}-2r'h, 


Hierin ist wie früher z die Koordinate in der Breitenrichtung mit dem Ursprung in der Mitte. 
Für die weitere Auswertung ist es vorteilhaft, die dimensionslosen Größen 


“ * ai 
= und [= 


V2r'h, 


einzuführen und damit Gl. (23) umzuschreiben in 


V2r'h, 


op, op 6£E dp 12nUE ‚v2 8 


= = E 5 . A 
actaetrireads 4 m on 
Die Benutzung der gleichen Ansatzform . wie bei Gl. (6) 
Cos u,£ DR \ Cos un-,° Cos u„£ 
sv  Cos u,£n, 7°) s u ‚co Cos Eu Pal)... 


»—=ı 


führt nach Einsetzen der erforderlichen Ableitungen in Gl. (24) auf ähnliche gewöhnliche 
Differentialgleichungen wie für den ebenen Gleitschuh, nämlich 
6£ 127 U, € 
i —: 9% 4 2r'- == 
r) Ptir: Po 1: yZr Ars) 
Cos u, ( re bE u 6E 


b) Cos u,» Mtjzemtm mn —- "jr8% 9) 0, 


Cos u, © 6£ 
ft | [7 IS )= 0, £ r (26), 


ce : — "x ea re ae en 
Cosa, u tige rr / 


allgemein geschrieben 


6£ 


(it td ı Yan — Yn-ı ir&Pn-ı — di ‚9u-,)=0 


n=1233..) 


Cos un-, 8 


Cos un, & 





worin ö, die durch Y2r’h, dimensionslos gemachte halbe Breite B/2 des Gleitschuhes_ be- 
deutet. Zur Unterscheidung ist hier für die noch unbekannte Funktion statt f nun @ ge- 
schrieben. Die Grenzbedingungen für die beiden Seitenkanten {= + {, des Schmierspaltes sind 
wieder p= Außendruck = Null. Sie werden schon durch die Klammerausdrücke der Gl. (25) 
befriedigt. 
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Die Lösung der Differentialgleichung (26a) lautet 


_„TyprhnJ3E—-d, 3 t 1:8#°+5) ERDE ‚| 
— h: \4 Ur ip + aretigs +c, 5 archg Ei +-cı} - 





(27). 


— 

-H 
Irre 
I 


Sie gibt die Druckverteilung des zweidimensionalen Strömungsvorganges für den unendlich 
breiten balligen Gleitschuh an, die an der Eintrittskante &, mit der Grenzbedingung 9, = Außen- 
druck =0O beginnt, deren Ende aber wegen der hinter der engsten Spaltstelle folgenden Er- 
weiterung nicht durch das Gleitschuhende festgelegt, sondern durch den Strömungszustand in 
diesem Spaltteil gegeben ist. Wie mehrfache Versuche zeigten, können Unterdrücke in dem 
erweiterten Schmierspalt durch das Vorhandensein von Luft oder Gasen in allen gebräuch- 
lichen Schmiermitteln nicht dauernd bestehen?) ®). Schon nach kurzer Betriebsdauer zerfällt 
der anfängliche Unterdruck, und der Außendruck wandert von der Hinterkante des Gleit- 
schuhes aus bis an den hinteren Rand der Überdruckzone vor. Dadurch ist aber nach 
W.Stieber’) bedingt, daß die Druckkurve an dieser mit £, bezeichneten Stelle mit horizon- 
taler Tangente in die Abszisse einmünden und gleichzeitig beendigt sein muß, was bedeutet, 
daß nunmehr zwei weitere Grenzbedingungen zo und 9,=0 bei £=£, gegeben sind. 
&, liegt noch nicht fest, folglich können die drei aufgestellten Grenzbedingungen durch die 
Gl. (27) mit ihren zwei willkürlichen Konstanten erfüllt werden, wenn £&, noch eine Funktion 
von £, ist, d. h. daß £&, nur ganz bestimmte von &, abhängige Werte annehmen kann. Mit 
dieser Bedingung ist nun die vollkommene Klärung der wirksamen Gleitschuhfläche erreicht. 


gr = Di i A 1 
Die Differentiation der Gl. (27) liefert nun mit der Grenzbedingung _ 


; d£ 
die Konstante e, =—6&. Die zweite Konstante ce, wird erhalten, wenn Gl. (27) für die 
Grenzbedingung 9,=0 bei E=£, angewandt wird. Damit lautet die Gleichung für die Druck- 
verteilung 


3 HU ABB) E45) Aa 3) & (+58) 
-7 EP A+E) 3 +&) | 


| 
+(1- 38) (aretg & — arctg | | 


Yo 
(28). 


Da nach Obigem auch bei E=£, 9,=0 sein muß, liefert der eckige Klammerausdruck von 
Gl. (28) die zur Bestimmung von £, benötigte Beziehung 


&1-35)-8,(1+58) 8A-35)—-&,(1+55) 
(1+ 5)? (1+£)? aaa ie 5, A 


+(1-38) (arctg &, — arctg &) =0 
In der Zahlentafel 2 ist die Abhängigkeit £&, von £, dargestellt. 


17 
Irre 





Ir 
... 


Zahlentafel 2. Werte von &, in Abhängigkeit von £.. 





vr 
- 
Kr 
Ir 
Ir 
In 
7 





2 1 2 1 o 
- 1,0 0,35787 — 24 0,45559 — 6,0 0,47355 
— 12 0,38794 — 2,8 0,46203 7,0 0,47413 
— 14 0,40944 -3,2 0,46593 — 80 0,47445 
— 1,6 0,42510 — 3,6 0,46842 — 10,0 0.47478 
1,8 0,43640 — 4,0 0,47012 00 0,47513 

— 220 0,44470 — 5,0 0,47246 


Bei den übrigen in Gl. (26) dargestellten Differentialgleichungen erfolgt das geforderte Ver- 
schwinden der linken Seiten durch Nullsetzen der Klammerausdrücke. Dafür werden wieder 
wie früher Lösungen mit der Ansatzform 


On HPn-ı , DOn-ı . al 


in Anwendung gebracht, die mit ihren Ableitungen die Suiginpitiiteihie eg 
(26b,c) usw. in die folgenden 
ur 
Hate +a-,.,=0 .., s=LR8: ME 
verwandeln. 


6) W. Frössel: Forsch. Ing.-Wes. Bd. Nr. 6 (1938), S. 261 bis 278. 
?) W. Stieber: Das Schwimmlager, Hydrodynamische Theorie des Gleitlagers. Berlin 1933. 
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An den durch die Gl. (29) festgelegten wirksamen Gleitschuhgrenzen &, und £,°) muß 
jedes Glied der Gl. (25) für sich Null werden, wenn ihre ganze rechte Seite verschwinden 
soll. Das kann aber nur erreicht werden, wenn jedes einzelne Glied der Funktion 9, = un-ı 
— @,„_, verschwindet. Damit hat die homogene Differentialgleichung (31) zwei homogene 
Randbedingungen zu erfüllen, die bekanntlich stets auf Lösungen von Eigenfunktionen mit 
den Eigenwerten «; führen. Die dritte Randbedingung, die einen tangentiellen Einlauf der 
Druckkurve in die Abszisse bei £=£, verlangt, kann durch die obige Differentialgleichung 
nieht mehr erfüllt werden. Hierzu wäre eine weitere von «; unabhängige Konstante erforder- 
lich, die auch in dem Ansatz Gl. (25) zum Ausdruck kommen müßte. Damit würden aber die 
Lösungsschwierigkeiten stark anwachsen und die aufzuwendende Arbeit für eine Lösung un- 
geheuer vermehrt werden, obwohl die erreichte Verbesserung selbst bei schmalen Gleitschuhen 
kaum merklich sein dürfte. Aus diesem Grunde blieb die Nichterfüllung der dritten Rand- 
bedingung bei der folgenden Auswertung unbeachtet. 

Die Auffindung einer allgemeinen Lösung der Differentialgleichung (31) ist mit Schwierig- 
keiten verbunden, deshalb wird eine numerische Auswertung nach dem bekannten Runge- 
Kutta-Verfahren vorgenommen, was zwangsläufig eine Einschränkung der Lösungen auf be- 
stimmte ausgewählte Grenzen &, zur Folge hat. Von den drei Konstanten einer aus der 
Differentialgleichung (31) hervorgegangenen Lösung kann wegen der schon erwähnten homo- 
genen Randbedingungen die eine als Faktor herausgezogen werden und zunächst bei der 
numerischen Rechnung unberücksichtigt bleiben. Die zweite wird durch den Anfangswert 
der nach dem angegebenen Verfahren errechneten Kurve bei &, ohne weiteres Zutun fest- 
gelegt. Dagegen ist die Ermittlung der dritten, die den Eigenwert «; darstellt, durch die 
an einer anderen Stelle bei £, vorgeschriebene Randbedingung sehr erschwert. Dieser Wert 
kann nur durch Probieren gefunden werden, 
was eine mehrmalige Errechnung der Kurve 
von einer Grenze bis zur anderen erforderlich 
macht. Aus zwei dem vorgeschriebenen Rand- 
wert dicht benachbarten Kurven wird dann 
die gesuchte Funktion sowie der zugehörige 
‚ a Wert durch Interpolation oder Extrapolation 
bestimmt. Das Bild 4 zeigt vier Muster für 






































die Bestimmung der Funktionen . " und 
n 1 

der u;-Werte bei den gegebenen Grenzen 

&=—2 und £,=0,44470. Wie leicht erkenn- 





bar ist, müssen die Rechenschritte des an- 
gewandten Verfahrens aus Genauigkeits- 
gründen mit steigenden w;- Werten immer 
kleiner werden, wodurch die aufzuwendende 
Arbeit beträchtlich erhöht wird. Zu jedem 
errechneten Punkt ergibt sich aus der Rech- 
nung gleichzeitig die Richtungstangente der 
Kurve an dieser Stelle. 

Das angewandte Runge-Kutta-Ver- 
fahren?) zur numerischen Bestimmung der 
Funktionen verlangt die Aufteilung der ge- 






























er gebenen Differentialgleichung zweiter Ordnung 
u — tar ; 
re Nick (Gl. (31)) in zwei erster Ordnung. Dazu kann 
473 2 i a 2 
1 FRE wen allgemein, wenn keine bessere Aufteilung mög- 
E = ANNIE lich ist, gesetzt werden 
\ zT 
ie | ZN he yo .. . .-. (88) 
N . . « 
dd: ach Gl. (31 
gi urn ua Yy—i— und damit nach Gl. (31) 
a & | 6: 
Bild 4. Lösungskurven der Differentialgleichung (31) y,= y, = | ng y; u’ er (32b), 
zur Bestimmung der «;-Werte, .® 


worin y wieder eine Funktion von & ist. 


8) Die Größe 52 ist nach der genannten Beziehung auf der ganzen Gleitschuhbreite konstant, was streng ge 
nommen für den endlich breiten Gleitschuh nieht mehr genau zutrifft. Für diesen müßte die wirksame Gleitfläche am 
Ende durch eine Randkurve begrenzt sein, also &3 erneut als Veränderliche von { angesehen werden, auf der die Grenz- 


ö 


h) 
bedingungen p=0 und = erfüllt sind, wobei 3 r den Differentialquotienten nach der Kurvennormale darstellt, Die 


ck 
unüberwindlich erscheinenden Schwierigkeiten, die einer solehen Lösung entgegenstehen, zwingen aber dazu, die erste 
vereinfachte Annahme beizubehalten. 
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Zur Ausführung des Verfahrens wurden die bekannten Formeln vierter Ordnung benutzt. 
Gegeben ist der Anfangswert der erwarteten Kurve o=0 bei &,=—.2. Ihre Steigung kann 
an dieser Stelle wegen der Homogenität der gegebenen Differentialgleichung (31) und der 
homogenen Grenzbedingungen beliebig gewählt werden und ist mit y, =1 angesetzt worden. 
Der Sprung k wurde zu 0,2 und «; probeweise zu 4,13 angenommen. Die fertiggerechnete 
Kurve ist in Bild 4 eingezeichnet. Sie trifft die «&-Achse noch nicht in dem gewünschten 
Punkt &,, so daß eine weitere Kurve mit einem besseren «,-Wert zu rechnen ist. Genügt 
auch diese noch nicht den Anforderungen, so wird aus beiden durch Interpolation die richtige 
ermittelt. Ebenso wird bei der Ermittlung der höheren Funktionen verfahren, die durch 
höhere u;-Werte gekennzeichnet sind und deren Kurven die «-Achse in dem betrachteten 
Gebiet ein- oder mehrmalig schneiden (Bild 4). 

Die Genauigkeit des angewandten numerischen Verfahrens hängt von der gewählten 
Sprunggröße-h und von der Krümmung der gesuchten Kurve ab. Die Fehlerabschätzung kann 
nach der bekannten Art erfolgen, nach der an der zu prüfenden Kurvenstelle eine zweimalige 
Durchrechnung des Verfahrens notwendig ist und zwar 1. in zwei Schritten mit dem Sprung h 
und 2. in einem Schritt mit dem Sprung 2h. Der Felıler zwischen der mit dem Sprung h 
gefundenen numerischen und der genauen Kurve beträgt dann näherungsweise '/ıs der Differenz 
der beiden Ergebnisse. Eine bequeme und ständige Überprüfung der Genauigkeit sowie eine 
Kontrolle für Rechenfehler liefert während des Rechnens der dauernde Vergleich des letzten 
Zwischenwertes ® + k“ mit dem in seiner Spalte folgenden Hauptwerte + Aw, die schon 
beide annähernd gleich sein müssen, wenn der Sprung h richtig gewählt wurde. Plötzlich 
auftretende größere Differenzen zwischen diesen beiden Werten lassen gewöhnlich auf einen 
Rechenfehler schließen. Aber auch bei stark anwachsenden Kurvenkrümmungen können im 
Laufe der Rechnung steigende Unterschiede auftreten, die dann eine Verkleinerung des 
Sprunges h erforderlich machen. Da jeder Schritt eine in sich abgeschlossene Rechnung 
darstellt, ist eine beliebige Sprungänderung nach jedem Schritt möglich. 

Die weitere Au gabe ist nun nach der anfangs gemachten Ansatzforım (Gl. (30)) aus 


jeder vorhergehenden Funktion, angefangen mit 9, e L. (28)), durch Hinzufügen einer geeigneten 

















anderen nacheinander neue Funktionen 9, = —- ®, M=P—W.:.... = Pn-ı — On, ZU 
finden, deren Funktionswerte mit steigenden In- f 
dizes in dem ganzen Bereich zwischen &, und £&, 7% ——T 
immer mehr gegen Null gehen. Zu diesem Zweck a 

w, j 
sind die schon bekannten Funktionen ©"! vor- —lL_ 
bereitet, deren Grundkurven durch Multiplizieren | 
mit der bisher noch beliebig verfügbaren Kon- | ——— N — 
stanten a„_, entsprechend umgeformt werden 25 = P; 
können. Als Beispiel zeigt Bild 5 die Entwick- pr 
lung der aufeinander folgenden Funktionen für ER Se 5 ı 
&=—2. Da eine derartige Rechnung einen be- | | | Fe 


trächtlichen Arbeitsaufwand erforderlich macht, 
sind zur Erleichterung der praktischen Anwen- 
dung auch die Funktionen für 43, =—1; —2; 
-4 und —6 ermittelt und in der Zahle ntafel 3 











\ 
8 
BE > Der in Gl. (28) auftretende KUN Ri 
& 
Ausdruck 2 jr Y2r’h, ist bei der bisherigen INNE 
Entwicklung her rausgezogen worden und kann #4, 7 REM - v ” 


nun bei der endgültigen Aufstellung der Druck- 


E z % Bild 5. Verlauf der Funktionen 9; und %; beim 
gleichung als Faktor zugesetzt werden. Damit wird ; 


balligen Gleitschuh („= — 2). 


3nU 2 / Cos u, Cos un_ı Cos u„£ E 
v| 4 y v2 " hs en (1 Vai 9 ( + >. | = > P )o n(S) ? (33), 


Cos HoCH © 0S Hn-ı I los Unsh) 





Pe | 
worin {, die halbe dimensionslose Gleitschuhbreite bedeutet. Aus der nun gegebenen Druck- 
verteilung folgt durch Integration die Tragkraft des Gleitschuhes, die zweckmäßig auf die 
Einheit der Gleitschuhbreite bezogen, dann die Tragkraft pro Breiteneinheit P’ in der dimen- 


sionslosen Form £& 
| 3 7 I 4 1 \ 

2 h, h . 

1 


liefert. s 


%) Runge und König: Numerisches Rechnen. Berlin 1924. — Schultz: PURE ENE zur praktischen Ma- 
thematik, Leipzig und Berlin, Göschen, Bd. 1110 (1937). Ein anderes Verfahren, das die Zerlegung der Differential- 
gleic hung 2 2. Ordnung in zwei l, Ordnung umgeht, ist vonZurmühlangegeben. Z.angew. Math. Mech. Bd. 20 (194v), S. 104. 
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Die Lage des Druckpunktes bzw. des Angriffspunktes der äußeren Last (Bild 3) errechnet 
sich mit dem Momentensatz zu 
1 Sa <h 
mn Cy \ \ | 
Ss 0 


tr 


o 


i 3 nÜ DR u. Wr 
& +» h: V2r'h,)d&d? ER 


re Mi -; 
worin &, = 75, h, Ist. 
Zur Ermittlung der Reibung wird die Schubspannungsverteilung benötigt, die nach dem 
schon benutzten Ansatz (Gl. (17)) unter Einbeziehung der Spaltgleichung (22) geschrieben 
werden kann 


(36). 


a dp dp ds 


. . v Fr = X wi 
Mit der Einführung der dimensionslosen Größe = 5,7, oder z—; =&h, und = 
” v2r'ı Ir dx dE de 


1 1 
= ar 3: wird dann 
- o - 
U dp 1 h,1+#) 


h, 1I+&) döy2rh, 2 
Die Differentiation der Gl. (33) liefert dazu 


E - u Cos un_,£ Cosu,f\ ‚| 

„+ ; Pos - er » I Pnf 
Pa 0S uUn-,%5 Losunc J 
n=] 


Cos u, 


dp 3 nt 
dE 4 M 





5! 
2rh 1 ; > 
l o\\ Cos u, Cn 
worin die @,-Werte die Ableitungen der Funktionswerte g„ bedeuten und die Tangenten an 
den Funktionskurven darstellen, die die angewandte numerische Rechnung — mit- 
liefert. Nur die Gewinnung der 9,- Werte erfolgt mit der differenzierten Gl. 28) 


#._ 2 v7) U ER | = ds & 
9=6b h: y:ı h, Te ETE 
Ihre zahlenmäßigen Werte sind in der Zahlentafel 3 aufgeführt. Mit der Abkürzung ®(£,{) 


. . dp 
für den geschweiften Klammerausdruck in 7 heißt die Gl. (36a) nun 


BE, A ES £2 al 36 
T= h ı+2 "2 (il +& )- D(E, 27 . . . . » . . . (36b). 


Die Reibung pro Breiteneinheit 7’ erbringt wieder die durch {, dividierte Integration der 
Gl. (36b). Sie heißt dann 


se 


Fr FR 
nU V; u: | g (1 +E)D(E, rasat Be. Ua A 
;\ 


- 


h, sb. 


sı 


III. Auswertung der Theorie. 


Es ist zweckmäßig, einen allgemeinen Überblick über die Wirkung der Gleitschuhlänge, 
-breite und -stellung auf die Trag- und Reibungseigenschaften des betrachteten Gleitsystems 
zu gewinnen. Während der Einfluß der Gleitschuhlänge und -stellung durch die in der Zahlen- 
tafel 1 angegebenen Funktionswerte f; und ; bereits gegeben ist, kommt die Beeinflussung 
infolge der Gleitschuhbreite unter mittelbarer Mitwirkung der Gleitschuhstellung durch die 
ArWerte erst bei der Errechnung der Klammerausdrücke in den Gl. (14) und (33) zustande. 
Die Klammerausdrücke sind als Beispiel für den ebenen Gleitschuh mit a=1,25 und den 
Breiten b=1,5, 1 und 0,5 in Bild 6 dargestellt. Sie lassen erkennen, daß der erste Klammer- 
ausdruck den Hauptanteil zu der Breitenwirkung beiträgt, der durch die nächstfolgenden je 
nach der vorgegebenen Breite mehr oder weniger verbessert wird. 

Die Kenntnis über die Druck- und Schubspannungsverteilung auf der Gleitfläche ist im 
allgemeinen nicht sonderlich wichtig. Die Verteilungen haben im wesentlichen alle die im 
Bild 7a und b für das Beispiel des ebenen Gleitsystems (a=1,25 und b=1) und die im 
Bild 7e und d für das Beispiel des balligen Gleitsystems (&,=—2 und (,=2) gezeigte Ge- 
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hr Bild 6. Einfluß der 


70 
a=-J25 | Gleitschuhbreiteinder 
Reehnung nach den 
08 fr Gl. (14) und (33 
ee -05 (a 1,25; b 1.5; 


t und 0,5). 
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Bild 7a. Schubspannungs- 
verteilung über der ebenen 
Gleitschuhfläche 
(.=1,8; b=1). 









Bild 7b. Druckverteilung N AR: 
über der ebenen s j 
Gleitschuhfläche N : 
(“=1,85; b=I) N 
KIER) N — a 
% 





Bild 7d. Druckverteilung 
über der balligen Gleit- 





schuhfläche (Sı = -?; <h=?). 
ENE Bild 7e. Schubspannungsverteilung über der 
balligen Gleitschuhfläche („= —?; p=?). aum \- 


stalt. Dagegen ermöglichen die durch die Integration obiger Verteilungen ermittelten 'Trag- 
und Reibungsgrößen eine leichte Beurteilung der Trag- und Reibungseigenschaften solcher 
Gleitsysteme. Es ist deshalb für eine schnellere und bequemere Bearbeitung eines Gleit- 
systems vorteilhaft, die obengenannten Größen für verschiedene Parameter b bzw. {, aus- 
zurechnen und graphisch darzustellen. Zugleich läßt diese Darstellung auch die Abhängigkeit 
der erwähnten Eigenschaften von Veränderungen am Gleitsystem sofort erkennen. 
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Bild 8a (links). Trag- BR; 
kraft je Breitenein- 35 1 
heit bei einem ebenen Ill 

Gleitsystem. 
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Für den ebenen Gleitschuh sind die Trag- und Reibungsgrößen nach den Gl. (15) 
und (20a)'°) errechnet und in den Bildern Sa und Sb in Abhängigkeit von der charakteristi- 
schen Größe a aufgetragen worden. Sie stellen, da die Gleitschuhlänge schon früher als Maß- 
einheit gewählt wurde, zugleich Werte für die Breiteneinheit des Gleitschulies dar. Als Para- 
meter für die Kurvenschar tritt hierbei die dimensionslose halbe Gleitschuhbreite b -5, auf, 
wo b wegen der zur Längsachse des Gleitsystems symmetrischen Druck- und Schubspannungs- 
verteilung von der Breitenmitte zählt. Die äußeren Kurven der Scharen im Bild Sa und Sb 
geben die Werte für die extremen Gleitschuhbreiten b=® und O0 an und begrenzen somit 
das in Betracht kommende Kurvengebiet. Die Beziehung für die größten Tragkräfte bei b=x 
kann aus der Integration der Gl. (8) gewonnen werden und lautet 

[6nU_,, @ 2 
ee Bet 
Die kleinsten Tragkräfte bei der Breite b=0 sind für alle a-Werte identisch Null. Die Kurve 
fällt mit der Abszisse zusammen. 
Die Integration der Gl. (17a) ergibt die Formeln zur Berechnung der größten und 


en ’ a. en R 

kleinsten Reibungswerte, wenn der Differentialquotient Fr entsprechend eliminiert wird. 
- = 

dp _6nTV | 1 2a(a—|) 


Jazu gilt für b=» Jer Ausdruck 
Dazu gilt für b=» Jer Au druck , Fi un 


)- Mit diesem Aus- 
druck heißt die Berechnungsformel 


Dr a 6 
| : E In, ı Ss —I 


a WR i dp x 
Für sehr schmale Gleitschuhe mit b>0 geht auch je >(0, und es wird 
2. 
ml ' a 
Tr| =—In 
r| D) 7 


I 
Das negative Vorzeichen bedeutet nur, daß die Schubspannung der positiven x-Achse entgegen- 
gesetzt gerichtet ist. Die letzten beiden Grenzkurven vereinigen sich bei a=1 und a=x 
U U ” \ a u 
und verlaufen von | a = nach tr 5 =». Sie umfassen den verhältnismäßig schmalen 
Bereich, in dem die ganze Breitenwirkung des Gleitsystems auf die Reib'ngskräfte zur Geltung 
kommt. Für große Werte a laufen sämtliche Kurven gegen Null, was aber nicht gleich- 
bedeutend sein muß mit r,—0. In diesem Falle geht auch ö>0, und die Reibung bleibt in 
dieser Form unbestimmt. Die wirkliche Lösung ergibt sich hierfür aus der Gl. (17) mit 


d I 
P_g und h=const zu 7, — 1 
dx h 


10) Das negative Vorzeichen der Gleichung gibt an, daß die Reibung der positiven x-Achse entgegengesetzt ge- 
richtet ist. Es blieb daher bei der Auswertung unberücksichtigt. 
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Für den balligen Gleitschuh lassen sich die Trag- und Reibungsgrößen, die in Überein- 
stimmung mit denen des ebenen Gleitschuhes auch wieder für die Breiteneinheit gelten sollen, 
in Abhängigkeit von der dimensionslosen Lagerlänge &, nach den Gl. (34) und (37) bestimmen, 
worin die dimensionslose Breite £, als Parameter auftritt. Die so entstandene Kurvenschar 
> 5 3 nU 
für P 5 h, 
breite £,=» ermittelte Kurve begrenzt, die für 2, >» sich asymptotisch dem Höchstwert 


r’ zeigt Bild 9a. Sie wird nach oben durch die bei unendlicher Gleitschuh- 


D 


3 nl . . . . \ 
pP |5 . r' — 1,6317 nähert. Ihre Bestimmungsgleichung ist aus der Integration der Gl. (28 
- 0 


hervorgegangen und lautet 


s T £3 2 2 £ LE 2 23 
PI5" "=1+( ad 35 (1455) +5: 
2 h, = . (1-+ &)’ 1+& 
+(1—-38)- £&,(arctg &, —arctg$,). 
x e , Re er . - ‚[3 U, . : 
Die unterste Grenze dieses Gebietes fällt wieder mit P’/z nt 0 mit der Abszisse zu- 
sammen. 
T 120 - ]%4 
* 1 + & a + 77V gen 
| | | 9agir \ a 
momen me tl En di ne } 4 4 | ia 
br ae | [ | 









12 - 124 
| PB zur! III | 7 er’ 

Ze > N co 
g8 N 6 

















ı 2 
un - = $ : = m - v 
6 um 5 4 6 u ” 4 3 & 2 
Bild 9a. Tragkraft je Breiteneinheit Bild 9b. Reibung je Breiteneinheit 
bei einem balligen Gleitsystem. bei einem balligen Gleitsystem. 


+) ' 
In Bild 9b wird die entsprechende Kurvenschar für die Reibungsgrößen T/n I | ; 
I, 


gezeigt. Auch hier wurde das negative Vorzeichen der Gl. (37) wie schon beim ebenen Gleit- 
schuh nur als Richtungszeichen der Reibung gewertet. Die beiden Begrenzungen sind wieder 
wie beim ebenen Gleitschuh durch die Breiten {,=0 und & angebbar. Ihre Beziehungen 
werden aus den Integrationen der Gl. (36a) erhalten und zwar für den unendlich breiten 


: > & dp de 
Gleitschuh mit Hilfe von rim ar aus der Gl. (28) zu 
Ip 3 + & 
me T er _ lo; m BEN Soc E areto £ : £ E S|l 
T/nl l u. 2 1,5 &}) (aretg &, — arctg &,) — > (& +8) 
. s „dp 
und für den sehr schmalen Gleitschuh £,>0 mit 1. >0 zu 
Br ‚ı/2r 
I int R — — (arctg &, — arctg $,). 
h, ji k 


Die letzte Gleichung gilt nur unter der Annahme, daß die zwischen &, und &, bestehende 
Beziehung, die ursprünglich für den unendlich breiten Gleitschuh abgeleitet und später auch 
für den endlich breiten übernommen wurde, auch noch für die Breite £&,=0 anerkannt wird. 


RE i . dp s a 5 ’ ’ ) 
In Wirklichkeit besteht bei dz>o0 zwischen &, und £&, keine Beziehung mehr. Damit verliert 
diese Grenzkurve ihre praktische Bedeutung und hat nur Sinn als theoretische Grenze des 


r . . N . ER . . y r 2 r' 
Kurvengebietes. Beide Grenzkurven beginnen bei & —=0 mit der Reibungsgröße T’| 17 I | z 06 
I, 


: pP 
una streben für £, >» asymptotisch gegen die festen Werte T’/ı, U | z — 3,6411 und 2,0144. 
7 


y)* 
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Der Abszissenwert &,—=0 ist gegeben, wenn x, = Gleitschuhlänge selbst Null wird oder wenn 
h,>x geht. Der erste Fall hat keinen praktischen Sinn, wogegen der zweite einem parallelen 


, ., dp £ ES FO j 
Spalt nahekommt, für den bekanntlich T —0 und h=const ist. Hierfür bleibt T’ bei An- 
wendung der obigen Formel unbestimmt. Die richtige Lösung ergibt sich wieder aus der 
U 


einfachen Schubspannungsbeziehung r,—= 7 


Die beiden Diagrammpaare nach den Bildern 8a und b und 9a und b ermöglichen nun 
ohne große Rechnung und Zeitaufwand die Gleiteigenschaften eines gegebenen Gleitsystems 
zu ermitteln. Insbesondere läßt sich hiermit der Einfluß einer Spaltveränderung auf die Trag- 
und Reibungskräfte feststellen. 

Beim ebenen Gleitsystem bewirkt eine Spaltveräuderung bei festgehaltenem «a eine 
Veränlerung des Winkels ö und umgekehrt bei festem ö eine Veränderung der Größe a. 
Während der Winkel ö als Faktor in die Rechnung eingeht und die Diagramme dabei un- 
berührt läßt, ist die Auswirkung von a auf die Rechnung aus den Diagrammen zu entnehmen. 
Alle Trag- und Reibungszustände des gewählten Gleitsystems bewegen sich dann auf den 
durch die dimensionslose Breitenabmessung b= B/21 vorgeschriebenen Kurven, worin B die 
ganze Gleitschuhbreite und I die Gleitschuhlänge bedeutet. Da die nunmehr mit a und b ermittel- 
baren Kurvenwerte, wie schon früher erwähnt, die Trag- und Reibungsgrößen pro Breiten- 
einheit angeben, errechnen sich daraus die absoluten Trag- und Reibungskräfte durch Multi- 

ner ei n 6nU n U Br 
plizieren mit der Breite B und den Werten Fe bzw. j: Hierbei ist nur zu beachten, 
daß bei zahlenmäßiger Auswertung für 7, U und B gleichartige Dimensionen in Anwendung 
gebracht werden müssen. 

Das Diagramm Bild Sa läßt erkennen, daß die Tragkraft des ebenen Gleitschuhes ins 
Unendliche wächst, wenn a gegen 1 abfällt, d. h. wenn der Austrittsspalt gegen Null geht. 
Diese Bedingung kann jedoch technisch infolge unvermeidbarer Oberflächenunebenheiten nicht 
ganz erfüllt werden. Gleichzeitig würde aber auch die Reibung über alle Grenzen ansteigen, 
so daß dieser Fall praktisch nicht durchführbar ist. 

Beim balligen Gleitsystem macht sich eine Spaltveränderung, gekennzeichnet durch den 
engsten Spalt h,, außer bei der dimensionslosen Gleitschuhlänge 5, = 


B 
2y2rh, 
der dimensionslosen Trag- und Reibungsgrößen für gegebene Gleitschuhabmessungen aus den 
Diagrammen Bild 9a und 9b ist daher im Gegensatz zum ebenen Gleitschuh mit einer 
Wanderung von einer Kurve zur anderen verbunden. Obwohl die Kurvenwerte für & =» 
endlich bleiben, werden die Trag- und Reibungskräfte selbst auch hierbei unendlich groß, da 
gleichzeitig h,>0O geht. 

Zur Bestätigung der vorstehend entwickelten Theorie sind Versuche durchgeführt worden, 
deren Beschreibung und Ergebnisse zusammen mit einem Auszug aus dieser Abhandlung in 
der Zeitschrift Forschung auf dem Gebiete des Ingenieurwesens'') erscheinen werden. Die 
Versuche zeigen in den Grenzen der Meßgenauigkeit gute Übereinstimmung mit der Theorie. 


2, : 
- ; auch gleich- 
verh, ” 


zeitig bei der dimensionslosen Gleitschuhbreite £,= bemerkbar. Das Aufsuchen 


IV. Zusammenfassung. 

In der vorliegenden Arbeit wird der dreidimensionale Strömungsvorgang untersucht, der 
sich in dem Sehmiermittel mit konstanter Zähigkeit zwischen einem endlich breiten Gleit- 
schuh und einer Gleitbahn einstellt, wenn zwischen beiden eine relative Geschwindigkeit 
parallel zur Gleitbahn herrscht. Die Untersuchung erstreckt sich auf zwei Gleitsysteme, auf 
den ebenen, der aus der ebenen Gleitbahn und dem unter einem kleinen Winkel dazu ge- 
neigten ebenen Gleitschuh gebildet wird, und auf den balligen, bei dem die Gleitfläche des 
Gleitschuhes konvex ist. Das Ergebnis ist die Berechnung der Trag- und Reibungseigen- 
schaften der genannten Gleitsysteme, die zwecks bequemer Anwendung in der Praxis auch 
zahlenmäßig ausgewertet und in den Diagrammen Bild Sa und b sowie Bild 9a und b dar- 
gestellt sind. Dadurch ist es leicht möglich, die Gleiteigenschaften der unter verschiedenen 
Bedingungen arbeitenden Gleitsysteme miteinander zu vergleichen. 

Herrn Prof. Dr. L. Prandtl sage ich für die Anregung zu dieser Arbeit und für die 
Unterstützung bei ihrer Durchführung meinen herzlichsten Dank, ebenso Herrn Dozent Dr. 
H. Görtler für die liebenswürdigen Beratungen bei den theoretischen Entwicklungen. 306 


11) W. Frössel: Reibungswiderstand und Tragkraft eines Gleitschuhes endlicher Breite, Forsch. Ing.-Wes. 
Bd. 13 (1942) H. 2, S. 65 bis 75. 
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Zur Theorie der kompressiblen Potentialströmungen Il. 


Intermediäre und singuläre Integrale der Grundgleichung der ebenen 
adiabatisch kompressiblen Potentialströmung. 


Von H. Behrbohm und M. Pinl (Messerschmitt-A.G.) in Augsburg. 


Mit Hilfe einer auf @. Boole zurückgehenden Methode gelingt die Inte- 
gration der beiden „integrablen Kombinationen“ der zur Grundgleichung 
der ebenen adiabatisch kompressiblen Potentialströmung gehörigen beiden 
Charakteristikensysteme. Damit sind zwei-partielle Differentialgleichungen 
erster Ordnung gewonnen. deren jede Lösung auch der Grundgleichung 
(zweiter Ordnung) für das Geschwindigkeitspotential des adiabatisch kom 
pressiblen ebenen Strömungsproblems genügt (intermediäre Integrale 

Die Integration dieser Differentialgleichungen erster Ordnung wird 
im nichtsingulären Falle nach Cauchys Charakteristikentheorie. im sin- 
qulären Falle nach Lagranyges Theorie des vollständigen Integrals durch- 
geführt. In beiden Fällen sind die Integralflächen z (x, y), weiche das (re- 
schwindigkeitspotential repräsentieren, Torsen. Im singulären Falle ver- 
schwindet die Schallgeschwindigkeit im Strömungsfeld, dagegen nicht im all- 
gemeineren nichtsingulären Falle. Für spezielle Konstantenwahl gehen die 
Resultate der Untersuchung in mehr oder weniger bekannte Ergebnisse aus 
der Theorie der Minimalflächen über. 


$ 1. Einleitung. 


Das Geschwindigkeitspotential 2=2z (x, y) einer ebenen adiabatisch kompressiblen Strö- 
mung mit der x, y-Ebene als Strömungsebene genügt der quasilinearen partiellen Differential- 
gleichung zweiter Ordnung 


A+elpP +) - pr - 2pqs+ß+tulp+g) Meer, ,. ch, 
wenn wie üblich zur Abkürzung 
02 02 0’ z 0’ 2 Ö’2 
p= q= = 


ig tw Sn m 
geschrieben wird, und wenn ferner mit Verwendung irgendwelcher Normierungswerte a, 
und ®w, für die Schallgeschwindigkeit a(x,y) und den Geschwindigkeitsbetrag m (x, y) im 
Strömungsfeld 

a 1—x 
0, / .) 


gesetzt wird; x ist dabei der adiabatische Exponent, so daß unter gewöhnlichen physikalischen 
Bedingungen u negativ ist. Eine Herleitung der Differentialgleichung (1,1) findet man in 
unserer unter [1]') genannten Arbeit. Wie dort führen wir auch hier die Bezeichnungen 
A=i+u(p +). C=it+lu (pt) 
ein und setzen noch 
R=A-p’=C+g, S=—p»qg, T=A—-d’=C+p, 


so daß (1,1) auch in der Gestalt 


R(p,9r+2S(p,Q)s+T(p.N)t=V0 .. . 12 
geschrieben werden kann. 


Zur Bebandlung der Grundgleichung (1,1) oder (1,2) werden zumeist Methoden voran- 
gestellt, welche die quasilineare Gleichung durch eine lineare ersetzen. Dazu empfiehlt sich 
entweder Legendres Berührungstransformation [2] oder eine Linearisierung durch eine ge- 
eignet verallgemeinerte Minkowskische Stützfunktion [3]. Demgegenüber versuchen wir 
hier einen direkten Weg zu gehen, wie er der klassischen Charakteristikentheorie [4], [5], 
[6], [7] entspricht. 


1) Die Ziffern in den eckigen Klammern beziehen sich auf das Schriftenverzeichnis am Ende der Arbeit. 
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Wir betrachten also die beiden Charakteristikensysteme ’) 


dz-pdx— gqdy=0, dy—r,de=0, dp+r,dq=V. . . . .(1,3a), 
dz—pdx gdy=0, dy-n,de=0, dp+rdgq=0. . . . .(1,3b), 
die durch die Lösungen 
) 7 ’ 
EEE u EIER. a 
( Ei 2, ( u 
der quadratischen Gleichung 
EEE es ne. 2 


näher bestimmt werden. Wir setzen zunächst 4 C==0 voraus. Dann sind die beiden Systeme 
voneinander verschieden. 


Für /=u=1 entsteht aus (1,1) die Differentialgleichung 
I+M)r—2pgas+(i+p’)i=0). . . PLN IT .- (4,6) 
der gewöhnlichen Minimalflächen. Ferner gilt in diesem Falle 


De ER en 2. 25205 Alu Aa en a kml. Hd Me 
A=1+pP+rd’, »- [+@ FE 7 i+g f 


Diese Spezialisierung ist für unsere Untersuchung charakteristisch: durch sie gehen alle A 
und « enthaltenden Formeln in bekannte Beziehungen aus der Theorie der Minimalflächen 
über [3]. Wir werden in $ 6 auf diesen Zusammenhang näher eingehen. 


$ 2. Die integrablen Kombinationen der Charakteristikensysteme. 

Nach (1,4) sind die Wurzeln », und », der quadratischen Gl. (1,5) Funktionen von p 
und g allein. Mithin ergeben sich aus (1,3a) und (1,3b) die beiden integrablen Kombinationen 
dp ö „„Pa+rVyPF—l(CHPICHN) 
dg P =» C+q - 

In beiden Fällen folgt 
p’(C+)—2ppga=—-(C+P), ep" +1)=(ar — pP. 


Zur Integration dieser gewöhnlichen Differentialgleichung erster Ordnung setzen wir*) 


. 6, Dy > ) 
p=osng(o), q=0C0S9@ (0); 0° = p” = u, o= arctg 1 
N U e & & q 


und bekommen 
dy\? 
== or| / . 


do 





— [4 +(u Deli+e (gr) 


do 


Vp?+q? 


BEE 2d 
9= arete 7 = E+4# \ V re —i)g ® + eonst 
q 0 


} A+ u0? 
Vr+t 
nach geeigneter Wahl der Anfangswerte p, und q,*). Das sind aber zwei partielle Differential- 
gleichungen erster Ordnung mit der Eigenschaft, daß wie jetzt gezeigt werden soll — jede 


ihrer Lösungen auch Lösung unserer Grundgleichung (1,1) ist. 
Sei nämlich 2 (x, ,y) eine Lösung der partiellen Differentialgleichung 


Vp?+g? 
i P _.[1/Ä+(u—1)o?do 
arcte — _M = onst — 5) 
F (p.g) = arctg g i \ | 17:3 s +eonst=0. . . .... 1). 
Vote 


2) Vgl. [4], Bd. I, Kap. Il, $ 50, S. 100. 
») Vgl. [7], S. 1%. 
% Es sei etwa(p + a) | # +u (pP? + A) 0. 














XUM 


Z. aungew. Math. Mech. 


Bd. 21 Nr.6 Dez. 1911 Behrbohm u. Pinl, Zur Theorie der kompressiblen Potentialströmungen 343 





Dann bestehen identisch in x,y (was wir durch die den Formeln beigefügten geschweiften 
Klammern fx, ,} andeuten) die folgenden Relationen°): 


F(p, g) v0, Fi. r+ F, Ss v0, F,s - F\, t v0, (F, F, ==0 rt 8° 0) (x, y, 


und es gilt somit 


F, F, 
(A-pP)r+A- Ma s=0, (A-N nn s+A— PMt=0 ,Y) - . (2,2). 
F F, ini 
Durch Addition der beiden Identitäten (2,2) folgt 
| F. F,| 
(A—p?) r+1(A PP) +A- Ms +A—- Mt=0 1, (2,3) 
| F, F, b 
Nun gilt, falls wir noch zur Abkürzung + W= + - setzen, 
‚ ..qFipW ; ptiqW PR, yFipW 
F, : ,.. = Fer - - . Ix,yX. 
P+q P+q F, ptigM ne; 


Hieraus ergibt sich nach einer leichten algebraischen Umformung für den Koeffizienten der 
Ableitung s in (2,3) 
qFipW 


F F, ptiqW 
) : 3 
ptiq W 


2 4 2 2 
(A P?) +(A-—-) F (A— p?) ipW 


p 


x 


(Ag 2p4 12,4}. 


+14 


Damit ist also in der Tat gezeigt: die beiden (Scharen von) partiellen Differential- 
gleichungen erster Ordnung 











Vor+a: 
p . i+(u—1)o?do 
arctg i RT Ui” - +-const = 0 . (2,4a), 
> q \ | A-+ u0° ee 
Vri+as 
Vp?+a? 
u do 
arctg ! +i\ | m 1) e IE rel. . . . (24b) 
a, 7 u0° 0 
V» + 











definieren zwei intermediäre Integralscharen der Grundgleichung (1,1). Wie 
man durch partielle Differentiation nach 4 leicht feststellt, entspricht dabei (2,4a) dem 
System (1,3a) und (2,4b) dem System (1,3b). 


Bis jetzt hatten wir immer AC==0 vorausgesetzt. Für AC=0 fallen die beiden 
Charakteristikensysteme (1,3a) und (1,3b) zusammen. Insbesondere ergibt sich für A=0, 
C=0: 

BEER 1 2.24: 0 Be dim nA ann 


und daher für jede Lösung 2 (x, y) dieser partiellen Differentialgleichung erster Ordnung 


i/+u(p +) 0, pr+qgs=0, ps+gt 0. ri +—=B. 


R(p,9r+2S(p,Qs+T(p,Qt=0 (2,9% °). 
Daher gilt: die Torsen (25) sind singuläre intermediäre Integrale der Grund- 
gleichung (1,1). 


Demgegenüber führt die Differentialgleichung C=4+(u— 1)(p+q)=V nicht auf 
intermediäre Integrale von (1,1). Während die Differentialgleichungen (2,4a b) der nichtsingu- 
lären intermediären Integrale von einer additiven Konstanten abhängen, ist die Gl. (2,5) der 
singulären intermediären Integrale durch die Konstanten 4 und « völlig bestimmt. 


5) Partielle Ableitungen nach den Argumenten p, 9, ©, 7, #, v deuten wir durch Indizes an, z. B.: F Fp Bi 
X%ı0,Xzr; Ableitungen nach den Argumenten 7 und « zeigen wir öfters auch dureh Striche an, z. B.: w’, f’, o*, u® yeye 


6, Vgl. [1], S. 202 (der zuerst genannten Arbeit). 
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$ 3. Die nichtsingulären intermediären Integralflächen. 


Wir betrachten das Charakteristikensystem der partiellen Differentialgleichung erster 
Ordnung (2,1) und erhalten: 


0. . EN. 4 — ptigqW a 2 A Be 
do ul de ff" de °° P?+qf' do rt riW, 
(3,1). 
dp \ dy i 
re (pF,} F,)=®, 2 (4 F,+F,) 0 | 


Nach den beiden letzten Gleichungen dieses Systems sind p und 4 vom Kurvenparameter o 
unabhängig. Durch Integration erhalten wir entsprechend (für beliebige Differentialgleichungen 


F=0 der beiden Scharen) die vierparametrige Streifenschar ’”) 


yFipW \ ptiqW | 
(ne): +06 r ze. U =(, oO Pr} > 
Y ( ‚rc P+q y(o Cs ’ P+q (32). 
se), EFioW, PO ae 


längs welcher die Funktion F zufolge (3,1) bekanntlich einen konstanten Wert annimmt 
dF 

(do 
schreiben, gewinnen wir die zu der bestimmten Differentialgleichung (2,1) aus den beiden 
Scharen jeweils zugehörige dreiparametrige Schar „charakteristischer Streifen“. Ihr Träger ist 
der Geradenkomplex 


0). Indem wir längs (3,2) speziell das Verschwinden dieses konstanten Wertes vor- 


alien 3 ee 


wo die Vektoren x, 4, 3 durch die folgenden Komponenten gegeben sind: 


ai Bi >; = Yır Yaı Yar > 6, , 6, “| 
ER EN iq # ipW ptig m. zw (3,4). 
> Ba 2,77 | pP+q ) pP+r FR | | 


Die Integralflächen von (2,1) sind spezielle Regelflächen im Geradenkomplex (3,3). Um sie 
zu bestimmen, entnehmen wir der Streifenschar (3,2) die nur von dem einen Parameter 7 ab- 
hängige Streifenschar 


ey M)+o5ln, 4=fl, 4=I9M: 1=ulM, =, Bull). . (85) 
und betrachten die „Anfangskurve* 
(0; r)=Yir). 


Die Anfangskurve Y(r) ist Träger eines „Anfangsstreifens“, wofern für die Ableitungen ce). (k = 1,2,3) 
die Relation 
dee) +glN ec; (?) ea AT 


besteht. Sie ist darüber hinaus Träger eines „Integralstreifens“, wenn für f und g die Relation 
Vr+o® 
F (f(r).g (r)) = arctg ö Fil | 


VR+9 


7 2do er 
e + (m 1)  — —- cont=0 A 
A-+ uo” 0 
besteht, und dieser Integralstreifen ist weder ein charakteristischer Streifen noch ein Fokal- 


streifen, falls 
gsifW, ftigW 


OT) = (a, Yr — Yo) P+g C, P+g ci =E0 N 


gilt. Diese letzte Ungleichung garantiert ferner für die Matrix 
%s Mor 2a 
I I 8% 


in einer gewissen Umgebung von o=0 den Rang 2. Mithin sind wir berechtigt, unter unseren 
Voraussetzungen der einparametrigen Streifenschar (3,5) die Trägerregelfläche x (o,r) zuzu- 
ordnen, und die Cauchysche Charakteristikentheorie liefert uns das Resultat: genügt die 


?, Vgl. [6], Kap. II, $ 3, S. 66. 
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Raumkurve y(r) der Bedingung (3,6) und befriedigen die Funktionen f(r) und 
g(r) die Bedingungen (3,7) und (3,8), so sind die diese Anfangsbedingungen 
erfüllenden nichtsingulären intermediären Integralflächen der Grund- 
gleichung (1,1) durch die Regelflächen 





gaFifndWin 
Fn+u’() 

‚dzisoWwn TREE i+ (u DH) ... (89 

P+r(n i+ uf? +) 


2,(0,1T)=c,(T) FioW(r) 


x, (o,r)=c,()+o 


2, (9,1) =c,(?) 











gegeben. 

In der Tat gestattet die Voraussetzung (3,8), die Funktionen z(o,r), p(r), q(r) als 
Funktionen 2 (x, 9), p (x, 9), q(x, 9) darzustellen, und die allgemeine Theorie zeigt, daß die 
so entstehenden Funktionen p (x, Y), 4(x,y) mit den partiellen Ableitungen z, und z, der so 
entstandenen Funktion 2 (x, y) identisch sind®). Die Regelflächen (3,9) sind also Lösungen 
des partiellen Differentialproblems erster Ordnung und nach den Ergebnissen von $ 2 auch 
solche des partiellen Differentialproblems zweiter Ordnung. Diesen Sachverhalt wollen wir 
noch durch die folgenden geometrischen Betrachtungen bekräftigen. 


$ 4. Die Gaußsche Krümmung K und die mittlere Relativkrümmung H*. 


Die Differentialgleichung F=0 enthält die Veränderlichen «&, y, z nicht explizit. Ihre 
Integralflächen sind daher notwendig Torsen. Eine zusätzliche Bestätigung dieser Tatsache 
erhalten wir durch das Verschwinden der Gaußschen Krümmung K (o,r) der Regelflächen (3,9). 
Wir berechnen sie als Quotient der Diskriminanten der zweiten und ersten metrischen Fun- 
damentalform der Flächen (3,9) ®): 

K h,, h;; Je hi. h (Ko, Kr. Yo) 


ae } Aka, Er, Lor) _ _(ko,&r, &rr) 
>, = =, A = RR 
Jr I2a — Ihe V 911 922 "I: 


Von " Vouden—- Ih 


Die Diskriminante 9,, 91: — 9, der ersten Fundamentalform verschwindet in unserem Falle 
nieht identisch, da wir 1+p?+q?=20 1x, y} voraussetzen. Da der Flächenortsvektor x (o, r) 
nach (3,9) von o linear abhängt, verschwindet h,, und wir erhalten für die Diskriminante 


h,,h h’, der zweiten Fundamentalform den Ausdruck 


11 22 


_ (Eos Er, Ko’ __ Ge drt ori __ rd de)’ en 


2 2 2 
Jıı I22 Ir Yı1 Isa ir Jı1 I22 — Iı2 
Für die weitere Rechnung führen wir die Abkürzungen 


grifW no _ftiaW 


der Ig +g 2, (= FiW, 


[77 


mn,()=+i(Wz,—-Wz), m. ))=2,2,-2%,23. m()=+i(Wz —W’z,) 
ein. Dann ergibt sich für die Determinante (3,,3,9,) mit Rücksicht auf (3,4) und (3,6) 
Gn3 9 elti, W—-Wa)+feiz,- 23,2) +lti(Wai—2,W’)+g(e 2, — 2,2)]e;- 


Verstehen wir jetzt unter p,, Ps, ?, die Komponenten des Flächennormalvektors x, X L,, so gilt 


Pı=H%gotar Korksos Pa okır "Karkıor Pa Ko Kern Kırfaos 

> Ds 

p= Pi R u=- P3 

Ps Ps 

Nun ist nach Voraussetzung (3,1) 
Po — PP: I, P: pP a 
— = EePıo PıPso_g, = 4,=f Pre PrPseo_g 
Ps Ps 


*) Vgl. [6], Kap. Il, $ 3, S. 67 bis 68. 
9) Vgl. [3], Kap. IV, $ 41, 42, 44; S. 85 bis 92. 
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oder 
2, (m, -fm,)ec +l2, m, +m,(+iW+9g2,)]e;=0, | 
u (4,2). 
(2, f+iW)m,-+ 2, m,] ci + 2, (m, — 9gm,) ce, =0 | 
Diese beiden Bedingungen sind proportional, da ihre Determinante (bezüglich ec} und c;) ver- 
schwindet, und jede von ihnen ist mit dem Verschwinden der Determinante (4,1) äquivalent. 
Da (4,2) nach Voraussetzung besteht, verschwindet also notwendig (3,,3,9,) und daher auch 
die Gaußsche Krümmung X (o,r), wie dies für alle nichtisotropen Torsen charakteristisch ist. 

Die mittlere Relativkrümmung der Flächen (3,9) ist durch 


„_ (89,8) + (f, Lo: €) 


D) (£, Lo I.) (4,3) 


H 


gegeben, worin f einen willkürlichen Vektor bedeutet, für den nur (f,x,,X,)=E0 gilt. e(o,r) 
ist der Ortsvektor der zu x(o,r) gewählten Eichfläche bzw. der Relativnormalvektor der 
Fläche x (o,r) im Sinne der relativen Flächentheorie [8], [9], [10]. Das Bestehen der Grund- 
gleichung (1,1) ist mit 7*=0 äquivalent, wenn man den Relativnormalvektor e(o,r) gemäß 


1 1 


I u 


e={—-2u—1)p R+ulpe tn] °“, 2u- Def tue tm] °*, 
1 

A+ul#® +) "R+A- WW +HPN : : (44) 
wählt'°). Nach (3,1) ist der Vektor e mit den Komponenten (4,4) von o unabhängig. Daher 
verschwindet die Determinante (f,e,,X,). Jetzt wählen wir £f=£0,0,1}. Dann reduziert sich 
(f,X,,Xr) auf 

(oT) = X, Lyr — KyaKır- 

Nach Voraussetzung (3,8) verschwindet die Determinante A(o,r) nicht identisch. Für die Be- 
rechnung von (f, x,,e,) ergibt sich nach (4,4) und (3,9) 


) 
- aufien 1)yA 2 RR _1)p4 Er 
43 


T 





1 1 
2 1 5a en zu 5 
—— Pr as rw. u” \-Hifzig w» 4° ) (4,5) 
1 
ui) u ge ae ee 2 
„gs A A (f+oas)tziAWlgf-gPN) 3 u + 
Nach (3,7) und (2,1) gilt 
Vreun+s®:n 
et, a A+(u—1)o?’ do . 
F(f(r), g(r)) arg +i \ | 4 r r e 0 + const = 0 er}. 
Vr +9 
Diese Identität leiten wir einmal nach r ab und erhalten 
er -ti,, rer Ä+ru-DEH+N_. 1 bw rw u 4,6 
Pry Fi F+g | It) e) zw. ti IFFgg er: Eh, 


Damit ergibt sich für den Klammerausdruck der rechten Seite von (4,5): 

AU-WY- (PN Ä+u HN -B+u- DPF -PH+M=0 An. 
Somit verschwindet die Determinante (f, x,,e,) identisch und daher ebenso auch die mittlere 
Relativkrümmung H*. 

8 5. Die singulären intermediären Integralflächen. 


Zur Integration der Differentialgleichung (2,5) der singulären intermediären Integral- 
flächen empfiehlt sich die Lagrangesche Theorie der vollständigen Integrale‘), Der mit 


10) Vgl. [10], $ 2, S. 81. 
11) Vgl. [6], Kap. Il, $ 4, S. 70 bis 72. 
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(2,5) verknüpfte Mongesche Kurvenkomplex ergibt sich durch Elimination der Ableitungen 
p und q aus dem System: 


En: BL ‚ dz p4,+44 P+q 


Es rc Ari u "a et” A "1 A, p 


in der Gestalt der Mongeschen Gleichung 
DEREN > ee . (51). 
Wir betrachten sodann das vollständige Integral 


V (x, 9,2;a,,a,,a,(a,,a,))za,c+ta,y+ta,2+a, =; 


p= “Um > -, ala ta) tig —=0 


a, a, IE 
dA, Ad, | 


der Gleichung A=0, d. Iı. eine zweiparametrige Ebenenschar. Sie enthält die einparametrige 
Lösungsschar 


=] u, a,=i(l+u?), a,=—2 —- 9, «(a +a)+/a=0V, 
die wir zur integrallosen Darstellung des Kurvenkomplexes (5,1) verwenden. Indem wir für 
a, die willkürliche analytische Funktion —2iy (u) einführen, gewinnen wir diese integral- 


lose Darstellung aus dem System 


2 > a\ u . ‚ 
1—-u)z+il+u)y 2] ; u2z—=2iy(u) 


07 


u 
-2ucst2iu Yy 2] rn 2iy’ (u) 


2 r+2i Yy 2iy" (u) 


‚mit der reellen Systemdeterminante si] = durch Auflösung nach x, y, z in der Gestalt 


, ce BE: j nn .. ö 
= il UF 5 yzy—-uy-+ a A 2 (y uw”). . (5,3). 
Für die Darstellung gilt notwendig y’”’ (u)$0 tu}, wie man durch Berechnung der ersten 
Ableitungen x’, y', 2’ erkennt. Man erhält so alle Komplexkurven mit Ausnahme der so- 
genannten „diskriminierenden Lösungen “ ’*). Sie stellen die Fokalkurven der Charakteristiken 
von A=0 dar und sind als solche notwendig die Gratlinien der Integraltorsen (die Charakte- 
ristiken selbst sind die diskriminierenden Lösungen der Mongeschen Gleichung (5,1)). 


Die gesuchten singulären intermediären Integralfächen sind durch die 
Gleichungen 





. 1— u? A 
x, (u,v»)=xr(u)+ nr’ =ily u" a y") ey”, 


1+u? 1+uw 
.(uV)=yW)+try= y-uy+——y” +0 oy”, y"’E0 WW, |.. 69 
! : 7 re Sa AR 
z,(uv)=2W)tr!=-—i | Fe (y—uy') +1 N uva 
[2 


Ä / 











gegeben; sie hängen von einer willkürlichen analytischen Funktion y(w)ab, 
deren dritte Ableitung nicht identisch verschwindet. 


Zum Beweis berechnen wir die Komponenten 





Lou Kzu 2. i1/ A | ae zus Kıu # 1 7 N > or. Hr Cou ER 
— ( w)oy'"?, =5z (+u)oy'"”?, =invy 
7 r 
! Lv Xgv 


[X . 
2,0, Lv Lv Cıv 


12) Vgl. D. Hilbert: Math. Ann. Bd. 73 (1913), S. 93 bis 108. 
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des Normalvektors der Fläche (5,4). Bekanntlich gilt 


Lu KXzu ZıurXau l 
p= = =—;z = u”) 
Lo Var Xıv, Xov zu 
FzuyKıu Zur Xu Fe 2 
= - =—o | —(1+W) 
ee RE IR zu u 


und daher folgt 
A=i+tulP+P)=0 due, y}. 


Die Torsen (5,4) sind also tatsächlich Lösungen von A=0 und daher auch solche der Grund- 
gleichung (1,1). 


Um zu einer reellen Darstellung der singulären intermediären Integralflächen zu kommen, 
oO > . . ” = 
entnehmen wir dem vollständigen Integral (5,2) die einparametrige Lösungsschar 


a=+i] „ecosu, @, +i] > sinu, ,=+1l, a,=-yiw, ula+ta)+i/a=V fu}. 


Da «<0, sind die Koeffizienten a, und a, für reelle « jetzt reell. An Stelle von (5,3) erhalten 
wir die Darstellung 


.:2 / u (m: ’ [7 . u ’ . „ _ 7 
z=r+iV (sin uy' + cosuy”), y-i f (cos u y sıinuy), 2: =+tly+ty'). 


An Stelle von (5,4) bekommen wir die reellen Torsen 





2 u > ’ vorn 
z,(u,v)=+i | j [sin u y" + cosu y’”’ + v cosu (y’ + y'"")] 


.]/ u : Re ; (5,5) 
2,(uv)= i 1 7 [eos u y’ — sin u y v» sin u (y’ + y’"”)] 


2, (u, v)= +F[y+y” ++ y””)] 











mt A+u(P H)=0O fu,v,y}, u<0. 


8 6. Übergang zur Theorie der Minimalflächen. 


Wir berechnen zunächst das Quadrat des Vektors 3=3(r) nach (3,2), (3,4) und (3,5) und 
erhalten: 


. PFr23ifgW—- PW’+P+2ifgW—- W’—-W?’(f’+ 9) I+f’+g 


Ägsi +9 nr 
grifW . . i . ö A n a 
Da 2, = "rg: nicht identisch verschwindet, können wir den Vektor 3* — _- bilden und 
. 


erhalten für seine Komponenten 


ftigW ar .P+g)W 
grifw "Tr giW 


2, = 1 5 U == 


. 


Für die Ableitung 23’ ergibt sich 


(gFifwi 


”r 


Nun gilt nach (4,6) 
(| I =+iW(ff +yyY). 


Daher erhalten wir für 25’ schließlich 


9 





Bi nr ei 27 


’ (gFifW) 




















YUM 
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Wie wir einleitend erwähnt hatten, gelangen wir zur Theorie der Minimalflächen, wenn 
wir die Formeln unserer Untersuchung auf den Sonderfall = u«=1 spezialisieren. Aus (6,1) 
folgt dann zunächst #?—=0. Ferner wird jetzt 
; 1 „_+ff+gg) 
+W _—— = . + u — if 99 
i+f+gP 


+ 


yiıfıg 


Sodann folgt aus (6,2) das Verschwinden der Ableitung 25’, also 2} =const= — /‘ Mit 3 ist 
natürlich auch gleichzeitig 3* ein isotroper Vektor. Somit gilt 


_ 


v2,’ +3 +, =1+l”+23°’=0, 23=+iyl+T?=const; 


oder mit Rücksicht auf (3,1), (3,4) und (3,9) 


dx dy dz_ ß N — rn x 
ER Fat T’:+iy1+[?”, eh 9, =B=9. 


Für die isotropen Kurven unserer Integralflächen ergibt sich jetzt"): 
(29,d0+9,dN)dr=0, 
für die Asymptotenkurven: 


Ru, 


für die Krümmungskurven 


dr, -dodr, do 
0, 92 + Ie | =geh,.d’=0. 
wi. 0 se 


Im nichtsingulären Falle (g,, 9. — 91, =E0) gewinnen wir so das Resultat: die nichtsingu- 
lären intermediären Integraltorsen (3,9) der Grundgleichung (1,1) reduzieren 
sich für A=u=1 aufdie isotropen abwickelbaren") Zylinder 


u 


Ju 
m 


L(o,r) =y)+oz3ld), 20, YyzElId 23:,:,=1:—-l:ziylil+l?; K=0to,r. 
Diese Zylinder tragen zwei Scharen isotroper Kurven, dagegen nur je eine 
Schar von Asymptotenkurven und Krümmungskurven, die beide mit der 
einen (geradlinigen) Schar der isotropen Flächenkurven zusammenfallen 
und die Charakteristiken des Differentialproblems darstellen. 


Wie man nach (4,4) leicht überprüft, geht die Relativnormale e für A=u=1 in die 
Flächeneinheitsnormale n über und daher nach (4,3) die mittlere Relativkrümmung H* in (die 
mittlere Krümmung H. Wegen H=K=0 haben wir damit die bekannten abwickelbaren 
Mongeschen Minimalregelflächen '’) gewonnen. 


Im singulären Falle (g,, 9: 9. =0) gehen wir von der Differentialgleichung (2,5) aus 
und gewinnen für /=u=1: 


ee ea N ar 


und daher nach (5,1) die Mongesche Gleichung 
I+y’+2:’=0 


des isotropen Komplexes. Das vollständige Integral (5,2) ist jetzt für die Differentialgleichung (6,3) 
durch eine zweiparametrige isotrope Ebenenschar gegeben (a +a3 +a?}=0). Aus (5,3) ent- 
steht die bekannte integrallose Darstellung der Kurven des isotropen Komplexes '*), aus (5,4) 
und (5,5) ebenso die der isotropen Torsen. Da jetzt «>0, sind beide Darstellungen der 
Integralflächen (und diese selbst) komplex. 

13) Vgl. [3], Kap. IV, $ 41, 52, 46; S. 87, 106, 94. 

14) Für diese Abwickelbarkeit ist die Voraussetzung Yıı 922 - 9,0 wesentlich! 

15) Vgl. L. Berwald: Sitz.-Ber. Akad. München 1913, S. 143 bis 211. 

16) Vgl. [6], Kap. I, $ 28, S. 45 bis 46. 
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Die integrablen Kombinationen der Charakteristikensysteme (1,3a, b) führen in der Theorie 
der Minimalflächen, d. h. für die Differentialgleichung (1,6) bekanntlich auf ein Clairaut- 
sches Differentialgleichungspaar mit dem allgemeinen Integral '’) 


p= aqFiyl-+a? SE N  — 


das wir als einparametrige Geradenschar in der Ebene der Koordinaten p und 4 deuten. 
. Mi . . R ; E { 
Setzen wır in der Differentialgleichung (2,1) A=u=1, so entsteht 


V»: +q? 
I 1 1+p? B 
arctg P +3 \ = „+teonst= Fit lg + ® p rq Lo" 
1 Jevil+e vPr+Ff 


Vri+a 
Gehen wir zu zyklometrischen Funktionen über, so ergibt sich: 


+i 


vr+q 


" a ) 
"; = sin Jaretg L_a* 
| q 


p y 2 i 
arcetg — = Farcesın = + 
vP®+q 


Daraus folgt 


peosa* - qsn—=Fi 


und diese lineare Relation in p und g wird mit (6,4) identisch, wenn wir die Integrations- 
konstante a*=arctga setzen. Differenzieren wir (6,4) nach a, so bekommen wir aus 


a=Fig(l +9) und (6,4) die Relation (6,3), d. h. das singuläre Integral des oben erwähnten 
Clairautschen Gleichungspaares. Auch so erscheint daher unsere Terminologie singulärer 
Integrale gerechtfertigt '?). 

Zusatz bei der Korrektur: Die Integration der intermediären Differentialgleichungen (2, 4a, b) nach der 


Lagrangeschen Methode und eine Untersuchung über die physikalische Natur der intermediären Strömungen 
erscheint demnächst in einer Arbeit von H. Behrbohm. 
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Turbulente Ausbreitung eines Stromes erhitzter Luft. 
Von Wilhelm Schmidt in Bochum. 
(Aus dem Kaiser Wilhelm-Institut für Strömungsforschung.) 
Il. Teil*) 
Das rotationssymmetrische Problem. 
Wir denken uns die Luft, die wieder durch ihre Dichteunterschiede gegen die ruhende 
Luft in Bewegung gesetzt wird, und die keinen Anfangsimpuls mitbringt, von einer kreis- 
förmigen Wärmequelle aufsteigend. Die x-Achse legen wir in Richtung der Strömung, die 
y-Achse senkrecht dazu. 
Die Auftriebskraft, die die Strömung hervorruft, ist 
gop». 
Hier gilt genau wie im ebenen Fall die erste Frage der Abhängigkeit des Mischungsweges, 
der Geschwindigkeit und der Temperatur von der Entfernung von der Wärmequelle Der 
Mischungsweg erweist sich entsprechend der Prandtlschen Abschätzung 
Db db Tu 
=u 
dt dz 
db 


dz 


b=r:x und !l=cx 


konst., 


als dem Abstand vom Ausfluß und damit auch der Strahlbreite direkt proportional. 

Um die Abhängigkeit der Temperatur und der Geschwindigkeit vom Abstand von der 
Wärmequelle zu finden, machen wir die Ansätze 9=Dx”-" und w=W- x? mit den zunächst 
noch unbestimmten Exponenten m und p und den willkürlichen Konstanten D und W, die 
lediglich bewirken sollen, daß die Ansätze dimensionsrichtig sind. Die Wärmetransport- 
gleichung liefert uns nun einen ersten und der Impulssatz einen zweiten Zusammenhang 
zwischen den unbekannten Koeffizienten m und p. 

Wir machen wieder die Voraussetzung, daß der Wärmetransport durch die einzelnen 
Querschnitte konstant ist: 

war huiih.: : ». ,„, ne nt a 


In dieser Gleichung wird o durch ersetzt (s. S. 269) und die Ansätze für « und 9 werden 
=) C 00 k 
eingeführt, so daß wir als erste Bestimmungsgleichung für m und p erhalten 


re Eee ee 3 


Die zweite Bestimmungsgleichung ist eine Folgerung aus dem Impulssatz. Der Impuls wird, 
wenn auch hier o durch o, ersetzt wird, 
J=nb’ou’ rn .”r+?? 
und die Impulsänderung 
dJ 


) p + R 
dx 


nm I# 


Die Auftriebskraft ist 
ab pgodran"t? 


Der Impulssatz: Auftriebskraft = Impulsänderung liefert 


.-Mm+2 r?P +1 


4 
und daraus als zweite Bestimmungsgleichung für m und p 
EEE EEE ZELTEN EEE SEP 


Aus den beiden Bestimmungsgleichungen errechnen sich m und p: 


Ü 
ja 


m = und See REN EEE 5 


wo 
.. 
— 


Unter Einführung der Dimensionslosen n7 = = bekommen wir für die Geschwindigkeit und die 
Temperatur die Ansätze 


*) Vgl. den I. Teil in Z. angew. Math. Mech. Bd. 21 (1941) Heft 5 S. 365 bis 278. 
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EN 
ST 5. ra EAST SEE N „7,0578 TH 
x 


u 
9 = — On). . (39b). 


Zur Ermittlung der Funktionen «, v und © stehen wieder drei Gleichungen, die Bewegungs- 
gleichung, die Kontinuitätsgleichung und die Wärmetransportgleichung zur Verfügung. 


Bewegungsgleichung: 
ou du 

u ® 
dx 7] 


« 


dudu 
TEN: 


? . 
I yy day) 





> fi 1 2 Ö 
ii A on 


Kontinuitätsgleichung: 


d d 
etz, erN=0 (4), 
Wärmetransportgleichung: 
is) 99 du] 
u Fe Tre Ta Ar Y7 35 || er 4 © 


Die Geschwindigkeitskomponenten « und v sind durch die Kontinuitätsgleichung miteinander 
verknüpft. Aus « kann v bestimmt werden. Nach dem Ansatz (39a) ist 
ce 
u=x *'f(n) 
und 
uy=nz:f(n). 


Setzen wir 
Inf) dn=F (1), 


dann wird 


„fa)=F(n) 
und 

2 
uy= °F’ (n), 


demnach 
iR" (m) 
a. en A ET EU 
) 
Aus der Kontinuitätsgleichung gewinnt man mit diesem Ausdruck für « die Geschwin- 


digkeitskomponente ®. Es ist 


Ö (23 pr 2 Ei Eee 
.3 = —r . s 
Aprie F') ser ern. 


Durch Integration nach y folgt daraus 


1 e vr 2 4 \ 
oy=z | F'ndy-— z# s\ Fady 





Ar - 37) 
also 4 = {F — ) ee a en ar 


Wir multiplizieren nun die Wärmetransportgleichung mit #g und bezeichnen das 


Produkt #9 © mit ©,. 
Um nur Gleichungen zwischen dimensionslosen Größen zu haben, nehmen wir eine 


Dimensionsbefreiung vor, indem wir setzen 


F=UyOß®; F=U,y9B®; FF=U,y9sB®"; 
9,=0:y(O,P° 9; = 12 YO, pP 9; A = U} YO, B} ©"; 
Bi. A U, f > ®\ U? 
e=—- VB ®; oe.) 9,B|® 3 IE ’ = V (9,ß’®- 
an a’ Dr. eg 
“ 





XUM 


ETRWITTN 
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U, ist die in irgendeiner Entfernung s von der Wärmequelle in der Strahlmitte gemessene 
Geschwindigkeit. Im gleichen Punkt ist auch die Temperatur #, zu messen, mit der dann 
auch ©, bekannt ist. Mit 


) 2 od’ u ’ D"' Do’ 
.. ya ac °Y in zn? 
u U,Y Mn _2P” | 2 
dp Ss Bar a 
und 
’ rg <q [6] ec od" DD’ 4 
= U TOR 85 =) 


Fi: % 
lautet die Bewegungsgleichung 


5PP N? —5DPP"'n’— D"n’ 


+30 6 BD’ n+2 0" Dr BD 3-20 Dahl 
und mit 
9. BIBI, u ): 9 _UYOP,, 9 _UYVIOP 
x Pg 3. 23 °Y P9 P A. 5 Pu 8 
und 
FE n (9, Pr e ( D"” v 9.6, 
Ppg #\n »# 


die Wärmetransportgleichung 


3.5299, +5P OO ?+PMn+P' n—- PD’ OO ?—- dB" —PO)—=0. . (46a). 


Um den Faktor 3c? auf der linken Seite von (45a) und (46a) zu beseitigen, machen wir die 
Substitution 
n=an mit a=(3e) \ 
und erhalten das Differentialgleichungssystem 
5PP m —-5P PD’ m — D’M+60D"’n,+20"D" nm \ 
-30": 330° 20 0" =3/ Be} 9,1: | 
5P 9,1, +5P 9! m+P' 9/1, + PP’ 9m, —- PP’ Nm — 


. (45b), 


PO —-PN=U0. . . (6b). 
Wird noch Gl. (46b) mit (3c?)® multipliziert und das Produkt (3c?)’ ©, mit ©, bezeichnet, 
so nehmen unsere Differentialgleichungen die endgültige Form 


5B DD —5BD"’ Den +6D BD’, +2D D-30308 20 Dt — 31} 


I) 


9, (4de), 


>5®P 9,1 +5P 9,1 +P 9n,+®" On, - PB’ 9% -—- PB" Om —- P%=0. . . . .(dbe) 


an. Zu diesen Differentialgleichungen kommen noch die Randbedingungen hinzu. Wenn wir 
(lie dimensionslosen Geschwindigkeiten mit «, und », bezeichnen, lauten diese: 
1. für die Mitte des Strahls 
G,=1, 
u,=endlich d.h. ®’=0, 
u.=® ee. 98, 
® muß Null werden wie n?. Nur dann hat «, in der Mitte einen endlichen Wert. 


2. für den Strahlrand ,=9,r: 


9,=0, 

0,=0 (wegen der Stetigkeit des Anschlusses nach außen), 

.-— aa. eh 

du, [23 % . . 

y = d.h. ©” —=0 (ebenfalls wegen der Stetigkeit des Anschlusses). 
i 


Die Lösung des Differentialgleichungssystems in geschlossener Form ist nicht möglich, 
sondern muß durch Reihenansätze erfolgen, in denen die beim ebenen Problem erwähnte 
Singularität zu berücksichtigen ist. 
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Die Ansätze lauten demgemäß 
P=an+an?’+a,m+ta,;n, 
9 


9, = 


18 
’ ta,m +a,n, 


3 
=-I+b,n?+b, 1 +b,n?+b,i + b,n, 
7 3 “ . 13 v r 
-20,+75.a,n’ 9,1 +5; n?+8a,n + 


19 
?—... 
15 

? th... 


. (47a). 


| 


19 15 
5 4a,n, 


2 


u . 
Werden diese Ansätze in die Differentialgleichungen eingesetzt, dann ergeben sich für die 


Koeffizienten die folgenden Bestimmungsgleichungen 
441 


da’ + 8 


a, =, 


161 2205 e 
Ta, +  u4,= 35, 
2 17 12285 TR TERS 
Ya,d, 3 - a? + g 4,0+ 11254 =3b, 
689 >45 En 22815 
"Ta, g 4% +3276a,a, + ga a,=3b, 
13689, 
‚+ 2» ‚=Bb, 


272 u,a, —218a, u, 


1577 a,a, - 1289a, a, 


ı0a,b, 7 
7 a,b,+ - a,b, - a,b, T a,b, —=0 
2>a,b,+235a,b, +25a,b = a,b, - 225a,b, nn u=0 
— T a,b, ee a,b, T a,b, - 468a,b, = 0 
1053 a, b, — 1152 a, b, 


65 65 
= a,b, + ... “ee 5 14 b, 1 5 a, b, 
0a,5,+Wa,b, +Wa,b,+40a,b, 


+40 a,b, 


DS 


6, +5 


55) 5) 
5 a,b, +3 a, b,+za 


6669 a, b, — 2024 a, 


Daraus berechnen sich die ersten Koeffizienten 


= 


do; 


1,6333 a + 4,7 


= 


a, 


ey = — 


= 


„. 1Wwa;+59855 a, a, + 11520 a, a 
1145a,a, +3%01 a, a,+ 81225 a,a, + 106704 a, a, = 12b 


Be vo 


63 
- a, b, —u 


3l5a,b, - 720a,b, 


- 850 a,b, =0 
9% 3% 3591 4275 
,b, +5 SL b,+z 34 u— ga b, g 9% b 
16245 5543 
a,Dd, 4 b, =VU. 


b, R 


2 _ 258 a,—+t 184 a} 


ive+a, = 


7884 at + 1,775} + 0,7218. 
(6 +8 a)? 


6+8a? 


2,2105 a} + 10,2863 a! + 5,9497 a3 + 19,1838 a3 + 0, 39 


5160 a; - + 720 


8300 as 3 
63 (6 + 8 a?)? 


20 a, . 
s Re — 


yv6+8a}’ 
1444,539a; + 1315,539 a) + 1029,805 a} + 217,776 % 


(6+8a}): 


5 





Die übrigen Koeffizienten wurden nicht explizit dargestellt, sondern durch Einsetzen der be 


reits berechneten in die Bestimmungsgleichungen gewonnen 
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‘du, (7, R) F r > : i 
Aus w, (n,rR)=0 und d . —=0( erhält man für a, und n,r ungefähre Werte. Es 
1 

zeigte sich, daß a, in der Nähe von 0,54 zwischen 0,54 und 0,55 liegen muß. Ein verbesserter 


Wert von a, ist 0,5405. Für n,r ergab sich 2,1. Damit werden die Koeffizienten 


a, = 0,5405; a, = — 0,2749912: a,=0,0641625; 
a, = — 0,0085605;  a,=0,000814919; a,-= 0,000055042: a, = 0,00000308857 
en SZ a 3 1,2477; bb, = 0,6801786; b, = -- 0,222; 
b,=0,049%135;  b,—= — 0,00834187;  b, — 0001099833 . 


Die Konvergenz der Reilienentwicklungen um den Nullpunkt wird bei Annäherung an 
den Randpunkt schlecht. Wir können aber wieder eine Entwicklung um den Randpunkt an- 
geben, die gerade in dessen Umgebung gut verwendbar ist. Wir setzen: 


N,=NıR ı 


und erhalten als Lösungsansätze für ®, «, und © 


Pat tan tim timt 

I . 2: . 
u, u 34,1, td, y —dä,n \ | (47b). 
9,=b,1 +bm+b +b,m+ + 


Durch Koeffizientenvergleich findet man Bestimmungsgleichungen für die Koeffizienten. Daraus 
bestimmen sich diese wie folgt: 

5: _  Ba;+ttb,ntr 
ü,; a.=- ö,; aü,= — - 
o 3 +) 4 2 

2 * 160a,n,.R 


5625 ä, 2400 dä; b, n\ 27 432 I 7 R 
dl, : 


s 120000 @ddyı r 
RR L Bd, +2 L 65, +10 +, 
Pn nn 2 — = . > - r\ 4 
1 ’ VA, NıR “ IN ETTR 
r 514275 a2 b, + 858500 at Bi 7? pr — 32400 a4 bin) a + 91584 bi ni pr — 843760 ab, mir 
’ 1500000 @ ir 


Läßt man die Entwicklungen um den Nullpunkt und den Randpunkt in einem An- 
schlußpunkt (n, = 1,6) übereinstimmen, dann berechnen sich die Unbekannten der letzten Ent- 
wicklungen zu 


ä, = 0,5057; b, — 0), 093: 7 — 051. 
Mit diesen Werten erhalten wir für die Koeffizienten der Entwicklungen um den Randpunkt 
die Werte 

a, =05057; Aa, — 0,2107; a, =0,058063: A, = 0,00496 ; 


b,=1; b, = 0,0581: b, = 0,02836 . 


Der genauere Wert des Randpunktes 7,r ist 2,11. Die Lösungen der Differential- 
gleichungen erhalten damit die Form: 


7 18 
D = 0,5405 1? — 0,2749912 n,? + 0,0641625 7° — 0,00856053 7,° ++ -- | 


u, — 1.081 — 0.962469 7, + 0,32081 7? — 0,0556435 7° ++ -- + + 
9, = 1 — 1.247948 m: + 0,6801786 n} — 0,222 2 + 0,0499514 n} | 
bzw. a 
D = 0,5057 — 0,2107 7? + 0,058063 7} + 0,00496 7 + + 
u, = {0,6321 7} — 0,232252 7} — 0,0248 7} .) un x (4b). 
9,=n; + 0,0581 7} + 0,0236 7, +: 
23* 
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In Zahlentafel 5 sind die Werte von ®, «,, v, und ©, angegeben. Wenn der Maßstab der 
Geschwindigkeit wieder so geändert wiıd, daß 

















Lak; | I FL los- | | die Scheitel des Geschwindigkeits- und des 
Z ET x ö . 
P- u N | | Auftriebsprofils zusammenfallen, erhalten wir 

En: TER die in Bild 5 dargestellten Profile. 

N | m u 

Br nn Das Temperaturprofil sinkt schneller ab, 
= II als das Geschwindigkeitsprofil. Die Querge- 
__ beschwindighai? ———— gg emeratur be schwindigkeit hat ihr Maximum nicht am 
Rande selbst, sondern etwas vom Rande ent- 


fernt im Inneren des Stralıls. 


F + + 0,4 — + — + 4 
ur 2 | Berechnung der Druckunter- 
. | schiede. Auch hier war in erster Nähe- 
u AR rung vorausgesetzt worden, daß der Druck in 
i TER PEARNEE den Querschnitten senkrecht zur Strömungs- 
richtung konstant ist. In zweiter Näherung 
sollen die Druckunterschiede berechnet wer- 








„2. BE 06 4 .. # #4” den. Wir gehen aus von der zweiten Be- 
Bild 5. wegungsgleichung 


(49). 


N) r : N) " > 1 Ö “ Y\ 1 N) (0, # Y) GT) | dp 
O4 oe\dz y dy 2 


u pr rt 
In dieser Gleichung sind neben der scheinbaren Schubspannung die Normal- und Tangential- 
spannung berücksichtigt, die sich aus dem Tensor 
„ou 
oy 


o0y 


I RE ER nee. ° 


bestimmen. Die anderen neu hinzuzufügenden Spannungen können vernachlässigt werden 
Für o, und “+ ergeben sich aus dem Tensor die Größen 
ou dv du v 


j] —»7? — 91? 
wer dylay' =2z,y 


u F dv : " re . R 
6, wechselt, da es das Vorzeichen von er annimmt, zweimal sein Vorzeichen, o;. das das Vor- 
zeichen von ® hat, dagegen nur einmal. 


Zahlentafel5. 





1) 1 





Ze 9; e „ad v | .08i) ” | 081) 
0 1 0 0 1,081 (1) 0 (0) 
0,1 0,9612 0,00532 0,1078 1,0775 (0,9751) 0,018 (0,01665) 
02 0.8937 0.0207 0.1995 0,9974 (0,9227) 0,0273 (0,02529) 
0.3 0.808 0,0466 0,2785 0,9284 (0,8588) 0.027 0.02498) 
0.4 0,7244 0,07598 0,3432 0,8579 (0,7937) 0,026 (0,024) 

0,5 0.6345 0,1127 0,3892 0,7785 (0,7201) 0,0135  (0,01245) 
0,6 0,5468 0,1533 0,418 0,6967 (0,6445) 0,0077  (— 0,00712) 
0,7 0,463 0,1967 0,43185 0,617 (0,5708) 0,036 (— 0,033) 
0,8 0,356 0,2412 0,4303 0,5379 (0,4976) 0,07 (— 0,067) 
0,9 0,3154 0,2815 0,4145 0,4605 (0,426) — 0,1069  (— 0,099) 

1 0,2529 0.3219 0,3897 0,3897 (0,3605) 0,1508 (— 0,1395) 
1,1 0,19995 0.3595 0,356 0,3237 (0,2985) — 0,1886 (— 0,1745) 
1,2 0,1518 0,3926 0,3134 0,2612 (0,2416) 0,233 (0,2158) 
1.3 0,1128 0,415 0,2656 0,2043 (0,189) | — 0,2677 (0,2476) 
1,4 0.0808 0.44605 0,2202 0,1573 (0,1455) — 0,311 (— 0,2877) 
1,5 0.0563 0,4658 0,1729 0,1153 (0,1067) — 0,3447 (— 0,3188) 
1.6 0,0373 0,4808 0,1277 0,0798 (0,0738) — 0,373 (— 0,3451) 
1.7 0,02285 0,4913 0,0846 0,04978 (0,0461) 0.3956  (— 0,366) 
1.8 0.0121 . 0,4978 0.0524 0,0291 (0,027) | — 0,409 (— 0,3784) 
1,9 | 0.00512 0,5026 0,0255 | 0,0134 (0,0124) | — 0,416 (— 0,385) 
2 0.001832 0.5034 0.0073 0.00365 (0.0033) - 0,4122 _ (— 0,381) 
2,11 0 0,5057 | 0 0 (0) - 0,399 (— 0,369) 


ı) Die eingeklammerten Zahlen sind die auf die Scheitelhöhe 1 umgerechneten Geschwindigkeiten. 
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Durch Integration der zweiten Bewegungsgleichung nach y erhalten wir 


Y u Y Y u 
1 s ie of * EA, pr "o 
„[pr= etz \ u edy+\ m dy : 5 7, dy+ lo, +\ er \ z dy. 


1) 0 0 1) u 


Setzt man in diese Gleichung für «, ®, T/.,, %, und o+ ihre Werte ein, dann erhält man durch 
partielle Integration für die einzelnen Summanden die Größen: 


a .„ BP 35#\]" 
2m ei 5 a. 2.» D’: Be 4 | 
1 \ x» m | 3 N 9 7 | 1 
y Yy ” 7 
KT: p: ny(O,P 5 Be N" aloe 00DD | 
N 2. I) dy+\ r dy i : | age 18 77° J +35 \ : dy 9 \ j’ d n, 


1) 0 £ N) 0 








Y , 
er .yd 3 u5 (910 PT" DD’ 2PP" 20"’o" 10B® 
ET Bar gi \\9 „7 a 4 y” 7 3 
1) z 0 
3er2y(9,B#2 DD De Deren ELAR AA, 
n " R _ ar Fr n r | Fr 3) ] 
Bi (3 7 9 Bnh I) Tg) 
u [1 
Yy ” 
| EZ 3eD3y (9% 3° 20 KDD WDD' WE" WED 
o\y*Y 3 Eu Er uk zu ae De ee De’ Dh 
#3 gi 
u u 
3eUtiy/(9, OD 5 5er (20er Ssor BB 
” Usyı er E a ( -+5 5 +\(3 a EM 6 
Pr {9 7 3m sh J\9 on 3n 
0 
Y 4 n 
a a TR ce 2 d® 0bP’ WEL 
r—jos,y=+r un S nYy 7 = an 
o)y 9 „a IE 3 7° Du, ve 
0 a I) 


327?%V (0.5 I DD’ 1VODD 5 Pr “ BE »5®: 
——-t ser l sy I, 8) (B —--5 : 5 - \(<®: E : \d» 
x 13 N yn 3 no I 3n | 


h 
Be u;y (BEP 160, 10DR", 1091007 


. 2 er’. 3 ‘ 3 57 4 . 
n 3 3 n 3 nie 


1 N 
- (,#=+ s 
% x L 
Die allgemeine Gleichung der Druckverteilung wird wegen ihrer Unübersichtlichkeit nicht 
angegeben. 
Die auftretenden Integrale sind keine uneigentlichen Integrale. Die Unendlichkeits- 
stellen der Integranden im Nullpunkt heben sich in ihrer Vielfachheit gerade auf. 
Der Druckunterschied zwischen Strahlmitte und Strahlrand berechnet sich zu 


.17?5/(M A? 
Pm PR _ 0,3749. (309) | |} (9sB) 0. 


0 
. a 


Im Inneren des Strahls herrscht ein kleiner Überdruck. Die Druckdifferenz gegenüber der 
ruhenden Luft ist: 


Pm Pe _ 0,9957 (30); IV (PR 


9 # 
Die Druckunterschiede sind so klein, daß die Geschwindigkeiten nur unwesentlich geändert 
werden. 
Stromlinien. Beim rotationssymmetrischen Problem ist die Stromfunktion definiert 
durch die Flüssigkeitsmenge innerhalb eines Kreises mit dem Radius y; es ist daher 
1 dy 1 dy 


u= ‚= 
2nyodx 


2ayoy' 
In unserem Falle erhalten wir demnach als Gleichung der Stromlinien 


5 
N a a Be a RE 
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Für einen bestimmten Wert von y erhält man eine & | | R 
Stromlinie. . Bild 6 zeigt die Stromlinien für I | \ 
a: ) 5 





1 3 6 » m 23 
76’ 16’ 16’ 16’ 16’ 16° 
3 ss: 5353 5 6 78 ) ; 
ı Zu, zild 6. 
16° 16° 16° 16° 16° 16 ac J % 
X 


/ 





, 


n 


























KZEIE 


Die dazugehörigen Werte enthält Zahlentafel 6. 





Experimenteller Teil. 


Um festzustellen, wieweit die theoretischen Ergebnisse mit den tatsächlichen Verhält- 
nissen im Strahl übereinstimmen, wurden Geschwindigkeits- und Temperaturmessungen im 
runden Strahl vorgenommen. Die Geschwindigkeiten wurden nach der Hitzdrahtmethode er- 
mittelt. Bei diesem Meßverfahren wird die Messung der Geschwindigkeit auf die Messung 
der Temperaturerniedrigung zurückgeführt, die das strömende Gas in einem elektrisch ge- 
heizten Draht hervorruft. Mit der Temperaturerniedrigung ist eine Widerstandsänderung 
verbunden, deren Größe, aus einer Brückenschaltung ermittelt, ein Maß für die Größe der 
Geschwindigkeit ist. Die verwendete Sonde, eine Zweidrahtsonde nach Reichardt*), hatte 
zwei Hitzdrähte verschiedenen Durchmessers, wodurch die durch die Temperaturunterschiede 
der Luft im Strahl bedingten Widerstandsänderungen ausgeschaltet werden konnten. 

Die Lufttemperatur wurde mit einem Thermoelement gemessen. Zur Erzeugung der 
Strömung wurde eine Heizspirale verwandt, die über einen Transformator geheizt wurde. Um 
eine möglichst ungestörte Strömung zu erhalten, wurde die Spirale in ein kreisrundes Loch 


in einem Laboratoriumstisch so eingebaut (Bild 7), daß die Oberseite der Spirale mit der 


Tischebene gerade abschloß. Zum Ausmessen des Geschwindigkeitsfeldes wurde eine Lanze 
mit eingebauter Hitzdrahtsonde und eingebautem Thermoelement benutzt. 


%, Beschrieben bei W. Paeschke: Zeitschr. Physik der freien Atmosphäre Bd. 24 (1937), S. 167. 


Zahlentafel 6. 








1 3 6 10 15 21 
—_— ww m y zum y m == ww = 
v6 "Te v6 a ET 2 

N = Y F 4 Y LT Y T Y € Y F 7 Y 
0.2 0.645 0.129 1.2665 0.2533 | 1.89 0.378 2.568 0.5136 3.275 0.655 ‚4.007 0.8014 
0,4 0.2952 | 0.1181 0.5707 | 0.2283 | 0.865 0.346 1.175 0.47 1.499 0.5996 1.834 10.7336 
0.6 0.1938 | 0.11628 0.3747 | 0.2248 | 0.5679 | 0.3407 | 0.7715 | 0.463 0.984 0,5904 1.204 |0,7224 
0.8 0.1476 | 0.1181 0.2854 | 0.2283 | 0.4327 | 0.3462 | 0.586 0.47 0.7497 | 0.5998 0.9175 0.734 
1 0.1242 | 0.1242 | 0.24 0.24 0.3638 | 0.3638 | 0.4943 | 0.4943 0.6304 | 0.6304 0.7715 0.7715 


1,2 0.1102 | 0.1322 | 0.2131 | 0.2557 | 0.3229 | 0.3875 | 0.4388 | 0.5266 | 0.557 0.6711 0.6847 0.8216 
1.4 0.1021 | 0.1429 | 0.1974 | 0.2764 | 0.2991 | 0.4187 | 0.4064 | 0.569 0.5184 | 0.7258 0.634 0.8876 
1.6 0.0976 | 0.1562 0.1887 . 0.3019 | 0.2858 0.4573 | 0.3885 0.6216 0.4956 | 0.79296 0.6064 0.9702 
1.8 0.095589! 0.172 0.1848 03326 | 0.2802 0.5043 | 0.3805 | 0.6849 0.4855 | 0.8737 0,5939 1.069 
2 0.09494 | 0.1899 | 0.1835 | 0.367 0.278 0.5564 | 0.378 0.7558 0.4821 | 0.964 0.5899 1.1798 
2.11  0.09469| 0.1998 | 0.183 0.386 0.2747 | 0.5853 | 0,377 0.7955 0.808 1.0145 0.5883 11.2413 





28 36 45 BR 66 78 
1) v = y’ = ) I ‚S 
7 TR v6 PT An v6 v6 
7 7 Yy T Y LE Y LE Y I Y LI Y 





0.2 4.758 0.9515 5.538 1.1076 | 6.332 1.2664 | 7.141 1.428 7.965 1.593 8,805 1.761 
0,4 2.18 0.872 2.535 1.014 2.898 1.1592 | 3.269 1.308 3.647 1,459 4.031 1.6124 
0.6 1.431 0.8586 | 1.664 0.9984 | 1.902 1.1412 | 2.144 1.2864 | 2,394 1.4364 2.646 1.5876 
0.8 1.0903 | 0.8722 | 1.268 1.0144 | 1.449 1.1592 | 1.635 1.308 1.824 1.4592 2.016 1.6128 
1 0.9168 0.9168 | 1.066 1.066 1.219 1.219 1.3746 | 1.3746 | 1.534 1.534 1,695 |1.695 
1,2 0.8138 0.9766 0.9908 | 1.18897 | 1.082 1.2984 | 1.2204 | 1.4645 | 1,361 1.633 1,505 1.805 
1,4 0.7538 | 1.055 0.8765 | 1.227 1.002 1,403 | 1,13 1.582 1.261 1.765 1.394 1.0536 
1.6 0.721 1.153 0.838 1.3408 | 0.958 1.5328 | 1,081 1.72896 | 1.2055 | 1.9288 |1.3324 2.1318 
1.8 0,7058 | 1.270 0.8206 | 1.477 0.9382 | 1.6876 | 1.058 1.9044 | 1.181 2.1258 1,305 12,349 
2 0.7012 | 1.4023 | 0.8152 | 1.6304 | 0.932 1.864 1.051 2.102 1.172 | 2,344 11,295 2,591 
2.11 | 0.6992 | 1.5652 | 0.8129 | 1.7145 | 0.9294 | 1.961 1.048 | 2,211 1,165 |2,458 11.293 12,728 
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Die Lanze bestand aus einem 1000 mm langen, 20 mm breiten und 4 mm dicken, flach- 
gedrückten Metallrohr, in das am einen Ende ein 120 mm langes, vierkantiges Stück eines 
Isolierstoffes etwa 40 mm eingeschoben war. An dem herausragenden Ende war die Sonde 
befestigt. Sie bestand aus zwei Blechen (Bild 8), an deren Enden zwei Hitzdrähte weich 
angelötet waren. Die gemeinsame Lötstelle der beiden Hitzdrähte lag an dem scharfkantig 
gehämmerten Ende eines Metalldrahtes, der an der Unterseite des zwischen den beiden 
Blechen liegenden isolierenden Verbindungsstücks befestigt war. Das Verbindungsstück, das 
die Aufgabe hatte, die Stabilität der Sonde zu erhöhen, trug auf seiner Oberseite, durch 
Isolierband isoliert, das Thermoelement. In Bild 8 ist das Thermoelement, dessen Lötstelle 
2,5 mm vor den Hitzdrähten lag, schematisch eingezeichnet. An den Blechen bzw. dem Draht 
an der gemeinsamen Lötstelle der Hitzdrähte waren Drähte angelötet, die zu den Anzeige- 
geräten führten. 

zum Galvanometer Sonde 
Hitzdraht 





Hıtzdraht 







Aurbeiwiderstand 
(/) Milliamperemeter 











NW 4 | 
Regulierwiderstand N akkumulator EHE) 


Bild 7. Bild 8. Bild 9. 


Die Hitzdrähte waren Platin-Iridiumdrähte von 0,021 und 0,041 mm Durelimesser und 
4 mm Länge. Sie lagen in zwei yerschiedenen Kreisen einer Wheatstoneschen Brücken- 
schaltung. Gegenübergeschaltet waren ein Festwiderstand von 1000 2, mit dessen Hilfe das 
Brückeninstrument, ein empfindliches Ruhstrat-Galvanometer, auf 0 einreguliert wurde. 
Die Regelung des Heizstroms für die Sonde, dessen Stromstärke 2 M.A. betrug, erfolgte 
über einen Widerstand. Als Anzeigegerät für die Heizstromstärke diente ein Milliampere- 
meter mit mehreren Meßbereichen, als Stromquelle ein Akkumulator (Schaltung Bild 9). 

Um den gesamten auszumessenden Strahlraum koordinatenmäßig aufteilen zu können, 
wurde die Lanze, die eine cm-Teilung trug, auf dem senkrecht durch einen Zabntrieb verstell- 
baren Tischehen eines Stativs befestigt. Die Höhenlage konnte an einer senkrechten em- 
Teilung abgelesen werden. 

Das Thermoelement, dessen eine Lötstelle in den Luftstrom ragte, und dessen andere 
Lötstelle in einem Ölbade lag, dessen Temperatur mit einem Thermometer bestimmt wurde, 
war mit einem Galvanometer verbunden, dessen Ausschläge mit t, bezeichnet, ein Maß für 
die Temperatur waren. 

Um störende Einflüsse von außen, wie Luftzug und Temperaturschwankungen abzuhalten, 
wurden Fenster und Türen abgedichtet. 

Die Messungen wurden wegen der Empfindlichkeit der Anordnung in windstillen Nächten 
durchgeführt. 


Eichung. Vor Beginn der Messungen wurde die Zweidrahtsonde an dem Eichgebläse 
des K.W.J. mit Staurohr und zwei Kleindruckwaagen nach Reichardt mit verschiedenen 
Meßbereichen geeicht. 












































780 i Zahlentafel 7. 
I | I 
| | 
| Tacho- Ausschlag Kurbel- Heizstrom 
| er C 5 der Druck- w m/s widerstand in MA 
775 Pe | | meter waare 11/4 ind ın M.A. 
| SE | 10 1,5 0.0942 156.2 200 
770 y, t } 1 et 
r | EEE 20 3,6 0,182 157,6 200 
| 
| | | 30 13,9 0,286 160,6 200 
u | R Zn a 50 39,7 0,485 168 200 
| /- | | | 70 82,2 0.697 172,9 200 
Pe EBEE er DEE Druck waage 
| | 11/6 200 
| E07 90 19,1 0.888 174.2 200 
. l il l J . 
". = SE ZEN. 110 29,0 1,095 175,8 200 
Bild 10. 130 40,1 1.289 177,6 200 
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Zur Berechnung der Geschwindigkeiten aus den Anzeigen der Druckwaagen wurde die 
Formel 


FT} 
m =] 15,15: z ug C:a 


verwandt. Dabei ist b der Barometerstand, T die absolute Temperatur, € der Eichfaktor der 
Druckwaage und a die Anzahl der Skalenteile der Druckwaage. Das Ergebnis der Eichung 
zeigen Zahlentafel 7 und Bild 10. 

Die Eichung des Thermoelements ergab, daß in dem benötigten Temperaturbereich einem 
Skalenteil am Galvanometer 3,4° C entsprachen. Aus den mit f, bezeichneten und in Skalen- 
teilen angegebenen Ausschlägen des Galvanometers bestimmt sich © zu 


9=p-ty, x. 


Ausführung der Messungen. Zunächst wurden die Geschwindigkeiten und Tempe- 
raturen in der Strahlmitte in den Höhen 40, 50, 60, 70 und 80 cm über der Tischebene ge- 
messen, um das Potenzgesetz zu prüfen. Dann wurden in den genannten Höhen die Ge- 
schwindigkeits- und Temperaturverteilungen gemessen und zwar jeweils von der Strahlmitte 
nach dem Rand hin. 


Die Messung sowohl der Geschwindigkeiten als auch der Temperaturen gestaltete sich 
dadurch schwierig, daß infolge der Verwirbelung des Strahls die Zeiger der Anzeigegalvano- 
meter nicht zur Ruhe kamen, sondern pendelten. Es mußten also die Schwankungen beob- 
achtet und daraus die Mittel gebildet werden. Nachdem die größeren Schwankungen, ver- 
ursacht durch das Verschieben der Sonde, abgeklungen waren, wurden aus jeweils ungefähr 
12 Ablesungen für Geschwindigkeit und Temperatur die Mittelwerte gebildet. 


Auswertung der Messungen: Die erste Aufgabe bestand darin, an Hand des Ver- 
suchsmaterials das Potenzgesetz zu prüfen. Für die Geschwindigkeit gilt die Gleichung 
‚ 


u=qx ’- (q = Umrechnungskonstante). 
/ 
Dabei ist x der Abstand vom Strahlursprung. Setzt man nun 
2,=2+410, ,=Z+5%9, 2,=Z+6W, 2,=2+ 70, r,=Z+8% 
und setzt für « die gemessenen Geschwindigkeiten ein, so ergeben sich die folgenden 


5 Gleichungen: 


„ 


iR 
0,64 = q(7440)°> „vs, 


a 


ı F 
0,61 = 4 (2450) ° (0,609), 
] 


Mn 4 

0.58 — q(Z+60) 3 . (0,58), 
] 
1 F’ 


0,56 = q (Z+70) m (0,55), 
1 F’ 
0,53—=q(2+80) °- (0,535). 
/ 
Für die Temperatur gilt die Gleichung 
d=Pp- «30. 


Setzt man auch hier für & die Größen &,, &,, ..., x, und für # , ein, so erhält man folgende 
Gleichungen: 

13,95=p(Z+40)50 (142), 

105 =p(Z+50730 (10,5%), 

82 -plZ+0 0 (819), 


65 =pl2+W 30 (66). 


59=p(Z+807 90 (65.2). 
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Eliminiert man nun in je zwei Gleichungen für die Geschwindigkeit g- und in je 
. e . ! 


zwei Gleichungen für die Temperatur p:© und setzt die verbleibenden rechten Seiten gleich, 
so ergeben sich aus den möglichen 20 Kombinationen 20 Werte von Z, aus denen sich durch 
Bildung des arithmetischen Mittels als Mittelwert für Z 


2153 


ergab. Dieser Wert erwies sich für die Auswertung auch der anderen Messungen als brauchbar. 
Der theoretische Ursprung des Strahls liegt nach dieser Rechnung nicht im Kopf der 
Heizspirale, sondern um 11,1 cm tiefer. Der Grund dafür liegt in der endlichen Ausdehnung 
der Spirale; in der Theorie dagegen war eine punktförmige Wärmequelle angenommen worden. 
Setzt man nun den berechneten Wert von Z in die Geschwindigkeits- und Temperatur- 
gleichungen ein, so lassen sich q und p daraus berechnen, und zwar ist 9g=2,4 und p = 10000. 
Die mit den ermittelten Werten von Z, p und q und den theoretischen Werten von ©, 


al 
errechneten theoretischen Werte sind in Klammern hinter den Gleichungen beigefügt. 
Das Ergebnis der Geschwindigkeits- und Temperaturmessungen zeigt Zahlentafel 8. Die 
mit der Sonde gemessene Geschwindigkeit ist w = yw’+v?’ m/s; « und » lassen sich bei der 
1 

Hitzdrahtmessung nicht voneinander trennen. Durch Multiplikation mit x® und Division 
durch 2,4 erhält man die auf das Maximum 1 umgerechnete und dimensionslose Geschwindig- 

1 1 

keit w- =? 0,417 wa®. 
Die am Galvanometer gemessenen t,-Werte wurden mit Hilfe des Eichfaktors 3,4 und 
der Temperatur der 2. Lötstelle des Thermoelements von 21,7 ° C auf Celsiusgrade umgerechnet. 
Durch Multiplikation der t, Werte mit x’ erhält man die auf das Maximum 1 umgerechneten 


5 
dimensionslosen Ausdrücke #,*. 


Zahlentafel®. 























z, —=B51,1 27, = 61,1 
y |w=Yu?+v? 0,417 25 w ta 9°C 0,0001 2, 23 v=Yyu?+v? 0,417 2 w ta .,»°C 0.0001 2, 15 
0 0,64 0,992 13.95 69.1 0.982 0,61 1.01 10,5 | 57,7 0.995 
2,5 0.595 0,929 13.37 67,2 0.94 0.58 0,96 10,1 56,0 0.958 
h) 0,41 0.64 8,96 52,2 0.63 0.46 0.754 7,58 47,5 0,742 
7,5 0,27 0.42 5,12 39,1 0,36 0.28 0,47 4,37 | 36,6 0,417 
10 0,16 0.30 3,13 32,4 0,22 0.19 0.32 2,67 30,7 0.262 
12,5 0,11 0,23 1,42 26,5 0,10 0,18 0,29 1,58 | 27,1 0,15 
2. = TEI , =81,1 
0 0,58 1.004 8.2 49.2 0.997 0.56 1.002 6.5 44.07 0,988 
2,5 0,55 0,952 7.78 482 0,953 0.54 0.956 6,22 42,9 0.95 
5 0,48 0,829 6.5 43,8 0.796 0.49 0,874 5,55 40,6 0,846 
7,5 0,34 0.588 4,7 37,7 0,55 0.40 0,721 4,36 36.5 0.66 
10 0.22% 0.387 2,78 | 31,2 0,34 0.29 0,524 3.34 33,1 0,46 
12,5 0,17 0,30 2,14 | 29 0,262 0.27 0,478 2,67 30.8 0,408 
15 0,13 0,224 1,31 | 26,2 0,16 0,23 0,41 1.60 27,1 0.24 
17,5 0,17 0,30 1.54 26.9 0.236 
= 9,1 
0 0,53 1,0 5,49 | 40,2 1,012 
2,5 0,51 0,96 5.08 39 0,942 
5 0,465 0,87 4.31 | 36,4 0.8 
7,5 0,40 0,75 4 35,3 0,744 
10 0,32 0,61 2.92 | 31,6 0,54 
12,5 0.25 0,46 2,64 | 30,7 0,40 
15 0,235 0.44 2,27 | 29,4 0.346 
17,5 0.16 0,30 1,31 | 26,2 0,24 
20 0,13 0,24 0,83 | 24,5 0.15 
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Bild 11 enthält die über dem Abstand von der 
Strahlmitte als Abszisse aufgetragene Geschwindig- 
keit w—=yuwu’+v’” und die Temperatur in Celsius- 
geraden. 

Um die Gültigkeit des Breitengesetzes 
nachzuprüfen, wurden die Breiten durch Division 








Y ’ ’ 
r durch x(2=n) dimensionslos gemacht und die 





20 75 [7 5 0 
-—— Abstand ıon der Strahlmitte aan 


. y 5 u 
dazugehörigen Werte von 0,0001 3 #, und 0,417 23 w 
Bild 11. 


in Bild 12 und 13 eingetragen. 

Vergleicht man nun die aus der Theorie sich ergebenden Breiten an den Stellen, an 
denen Temperatur und Geschwindigkeit auf die Hälfte abgesunken sind, mit den entsprechen- 
den Breiten von Bild 12 und 13, so läßt sich die noch unbekannte Konstante e gemäß der 
Beziehung A, 
a—=@e) * 
berechnen. Die Ausrechnung ergibt für ce den Wert 

ce— 0,0417 
und damit für den Mischungsweg I 
1= 0,0417 x. 
Die in der Gleichung b=r-x auftretende Konstante r bestimmt sich aus der Gleichung 
ee 
na a A" Fena r = 0,366. 
Damit ist b — 0,366 «. 

Die letzte noch unbekannte Konstante &, die in der Gleichung != eb auftritt, stellt sich, 

wenn für ! und b ihre Werte eingesetzt werden, durch die beiden errechneten Konstanten e 





: e - 
und r dar, und zwar ist e = = 0,1139 und damit 
— 0,1139». 
In Bild 12 und 13 sind die aus der Theorie sich ergebenden Werte von ©,, u, und 
w,=yui+v; über n als Abszisse zum Vergleich mit den gemessenen Werten eingetragen 
(9, ausgezogen, w, gestrichelt). 
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Der Vergleich ergibt nur für die mittleren Profilhöhen Übereinstimmung. In der Strahl- 
mitte und am Strahlrand ist die Übereinstimmung mit den theoretischen Ergebnissen schlechter. 

Im unteren Teil des Profils schwanken die Meßpunkte sehr stark. Die Erklärung dafür 
ist wahrscheinlich in der Tatsache zu suchen, daß dieselben Störungen bei kleinen Geschwin- 
digkeiten größere Abweichungen hervorrufen als bei großen. 


Zusammenfassung. 


Im ersten Teil dieser Arbeit wurde auf Grund von Ansätzen von Prandtl und Taylor 
für die scheinbare Schubspannung und eines Ansatzes für die turbulente Wärmeübertragung 
die Geschwindigkeits- und Temperaturverteilung für den Fall der durch Dichtedifferenz gegen 
die Umgebung hervorgerufenen ebenen Strömung berechnet. Bei der Benutzung des Prandtl- 
schen Schubspannungsansatzes wurde die Voraussetzung gemacht, daß die Austauschgrößen 
in der scheinbaren Schubspannung und der turbulenten Wärmeübertragung identisch sind, 
während nach Taylor die Austauschgröße der turbulenten Wärmeübertragung zweimal so 
groß ist wie die der scheinbaren Schubspannung. Die Lösung erfolgte durch Reihenansätze. 
Die Ergebnisse sind den Ansätzen entsprechend verschieden. Während nach Prandtl das 
Temperaturprofil stärker abfällt als das Geschwimdigkeitsprofil, das Temperaturprofil also ganz 
innerhalb des Geschwindigkeitsprofils verläuft, ist das Ergebnis der Rechnung nach Taylor 
umgekehrt. 

Die Theorie ergab für beide Annahmen, daß die Geschwindigkeit in der Mitte des Strahls 
von der Entfernung von der heißen Linie unabhängig ist, die Temperatur dagegen dem Ab- 
stand von der heißen Linie umgekehrt proportional ist. 

Die Berechnung der in erster Näherung vernachlässigten Druckunterschiede ergab für 
die Mitte des Strahls bei der Taylorschen Annahme Gleichdruck, bei der Prandtlschen 
einen kleinen Unterdruck gegenüber dem Druck in großer Entfernung. 

Im zweiten Teil wurde das rotationssymmetrische Problem behandelt. Die Theorie 
lieferte hier das Ergebnis, daß die Geschwindigkeit mit x s und die Temperatur mit x ® 
abnimmt. Für das räumliche Problem wurde die Rechnung nur mit dem Prandtlschen An- 
satz durchgeführt. Der Taylorsche Ansatz konnte hier nicht verwandt werden, weil er 
unter der Voraussetzung ebener Turbulenz abgeleitet ist und der analoge rotationssymmetrische 
Ansatz auch bei anderen Aufgaben versagt hat. Der „modifizierte Taylorsche Ansatz“ würde 
voraussichtlich gegenüber dem Prandtlschen Ansatz ähnliche Untschiede aufweisen wie im 
zweidimensionalen Fall. 

Auch hier lieferte die Berechnung der Druckunterschiede einen Überdruck im Inneren 
des Strahls. Die bei der Berechnung der Druckunterschiede auftretenden Integrale wurden 
mit Hilfe der Simpsonschen Regel ausgewertet. 

Der rotationssymmetrische Fall, der sich am leichtesten experimentell verwirklichen 
läßt, wurde durch Versuche geprüft, bei denen die Geschwindigkeits- und Temperaturverteilung 
mit Hitzdrahtsonde und Thermoelement gemessen wurde. Die von der Theorie geforderten 
Exponenten m und p wurden bestätigt. Die Geschwindigkeits- und Temperaturprofile stimmten 
dagegen nicht sehr genau mit dem theoretischen überein. Für die in der Theorie offen- 
bleibende Zahl e wurde der Wert c=0,04165 ermittelt, was !=0,1139b und b= 0,3655 x 
liefert. 
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Das statistische Problem 
der Korrelation als Variations- und Eigenwertproblem 
und sein Zusammenhang mit der Ausgleichsrechnung. 
Von Hans Gebelein in Eßlingen. 


Von einem voll befriedigenden Korrelationsmaß ist zu verlangen, daf es 
durch sein Verschwinden mit Sicherheit anzeigt, wenn die beiden Veränder- 
lichen unabhängig voneinander sind, und daß es anderseits mit Sicherheit 
den Wert Eins annimmt, wenn zwischen ihnen vollständige Abhängigkeit 
besteht. Die Aufgabe, ein solches Korrelationsmaß zu gewinnen, führt auf 
ein Variationsproblem und kann zurückgeführt werden auf die Frage nach 
dem kleinsten Eigenwert einer homogenen Fredholm schen Integralgleichung. 
Bei der Untersuchung dieses Variationsproblems wird auch der Zusammen- 
hang mit den bisher üblichen Korrelationsmaßen und mit dem Problem 
der Ausgleichsrechnung ersichtlich. 


1. Problemstellung und übliche Korrelationsmaße. 


Einer der bemerkenswertesten Gesichtspunkte für die Beurteilung zweidimensionaler 
Wahrscheinlichkeitsverteilungen w (x, ,) ist der, ob durch Festlegung der einen Veränder- 
lichen x die Verteilung der y-Werte mehr oder weniger stark beeinflußt wird oder nicht. 
Eine solche Beeinflußbarkeit besteht bekanntlich nicht, wenn » (x, y) das Produkt einer 
Funktion w, (x) von x allein und einer Funktion w, (y) von y allein ist. Die Merkmale x 
und y sind in diesem Fall voneinander stochastisch unabhängig. Andererseits kann es vor- 
kommen, daß zu jedem bestimmten Wert für.« nur ein einziger y-Wert auftritt. Dann ist y 
eine Funktion von x und man spricht von vollständiger Korrelation zwischen den beiden 
Merkmalen. Im allgemeinen liegt der Befund zwischen diesen beiden Grenzfällen und es 
entsteht dann die Frage nach einem Maß für die Straffheit des Zusammenhangs zwischen x 
und y. Dies ist das sog. Korrelationsproblem der Statistik. 

Um die Korrelation zwischen y und x zu kennzeichnen, gibt es bekanntlich eine An- 
zahl verschiedener Rechengrößen, die wir zunächst der Reihenfolge nach aufführen. Wir 
betrachten vor allem den Fall einer sog. geometrischen Wahrscheinlichkeitsverteilung für die 
Merkmale x und y. Eine solche wird durch eine wesentlich positive Funktion w (x, ,y) mit 
der Normierung 

BE ee 


dargestellt. Wird über x bzw. y allein integriert, so entstehen die beiden Funktionen 
,()=)\w(az,dy und w,W=\w(,ydz. . . . .. . 2), 


deren Integrale über «& bzw. y wiederum den Wert Eins besitzen. Mit a und b bezeichnen 
wir im folgenden die Mittelwerte für x und y auf Grund der Verteilung w (x, y); das sind 
die Schwerpunktskoordinaten der mit der Massendichte w (x, y) belegten xy-Ebene. Es ist 
a=\)\ew(e,yydedy=\rm,(e)de, | . 
2 Bi a War PET EN  AIE R 
b=\\yw(e,ydedy=\ym,(y)dy 
Weiter bedeuten s? und #? die Streuungen der « um ihren Mittelwert a bzw. der y um 
ihren Mittelwert 5; dies sind die beiden Trägheitsmomente der mit w (&,%) belegten Ebene 
um Achsen durch den Schwerpunkt und parallel zu den Koordinatenachsen. Es ist also 
"—\\(e - a’ w(e,ydedy=|\(r —- a®’w,(e)de, (4 
= |\\(y-bYwia,dedy=\ly- bYw,(y)dy 
Als Korrelationsmaße kommen Zahlengrößen in Betracht, welche erstens für jede vor- 
gegebene zweiparametrige Verteilung sich berechnen lassen, zweitens den Wert Null an- 
nehmen, wenn x und y stochastisch unabhängig voneinander sind, und drittens gleich Eins 
werden im Falle vollständiger Korrelation zwischen x und y. 
Die erste und zweite dieser Bedingungen erfüllt der auf dem Zentrifugalmoment der 
Verteilung beruhende klassische 
Korrelationskoeffizient: 


= 


\(e—-a)(y-b)wiz,„dady. TER ER GETEBITIT,:  ° 


1 
st 














XUM 


Z. angew. Math. Mech. R R ec i r BR 
Ba. 3 Nr.6 Des. 104 Gebelein, Das statistische Problem der Korrelation 365 





Bekanntlich ist stets r?<1 (vgl. unten Gl. (13)); der Wert r=+1 tritt auf, wenn y linear 
mit & zunimmt und der Wert r=— 1, wenn y linear mit x zusammenhängt und mit wachsen- 
dem x abnimmt. Der Mangel dieses Korrelationskoeffizienten ist der, daß auf Grund von 
r =( nicht auf vollkommene Unabhängigkeit mit w (#,y) =w, (x) w,(y) geschlossen werden 
kann, während anderseits trotz vollständiger stochastischer Abhängigkeit von x und y immer 
dann, wenn x und % nichtlinear zusammenhängen, r von Eins verschieden ausfällt. 

Namentlich um dem zweiten Übelstand abzuhelfen, wurden andere Korrelationsmaße er- 
funden. Ein geeignetes Mittel zur eingehenderen Analyse der Verteilung w (x, y) sind die 
sog. Regressionslinien y(x) und @(y). Es sind dies bei mechanischer Auffassung die 
geometrischen Orte für die Schwerpunkte der Streifen in der «y-Ebene parallel zu den 
Koordinatenachsen 


= \yalz,y)dy 1 
7(&) : m iyw(a,y)dy, 
\w(z,yJ)dy w,(a@)'' N 
(6). 
= \zw(e,y)d& 1 
ziy)=° : = \zw(e,y)Jd® 
\w(z, dx  mw,(y) 


In den beiden Grenzfällen stochastischer Unabhängigkeit und vollständiger Abhängigkeit 
sind die Regressionslinien deutlich gekennzeichnet. Bei Unabhängigkeit ist wegen w (x, y) 
—=mw,(#)w,(y), y(x)=b und z(y)=a, d.h. die Regressionslinien sind die waagerechte und 
die lotrechte Gerade durch den Schwerpunkt der Verteilung. Es ist jedoch zu beachten, daß 
aus dieser Lage der Regressionslinien, d. h. der Schwerpunkte nicht etwa umgekehrt auf das 
Vorliegen stochastischer Unabhängigkeit geschlossen werden kann, denn es folgt daraus 
nicht, daß mw (z,y)=w, (x) w,(y) ist. Anderseits ist bei vollständiger Abhängigkeit, wo die 
x und % einander eindeutig durch die Funktion y(x) oder & (y) zugeordnet sind, % (x) = y (x) 
und @(y)=x(y), so daß in diesem Falle die beiden Regressionslinien zusammenfallen. 

Auf diesen Regressionslinien beruhen die von Pearson in die statistische Praxis ein- 
geführten 

Korrelationsverhältnisse: 


Baia 
Mn, ra (y(xz) - bb)’ w, (x) dx (Korr.-Verh. von y zu | 


ki s? 


EN AED | 
z\(@(y)- a)’ w,(y)dy (Korr.-Verh. von x zu y) 


Diese Korrelationsverhältnisse nehmen den Wert Null bei Unabhängigkeit und vor allem den 
Wert Eins ganz allgemein bei vollständiger Verbundenheit sicher an. Daher wäre die Messung 
der Korrelation mit Hilfe dieser Größen sehr zufriedenstellend, wenn sie eindeutig wäre. 
Handelt es sich aber nicht um einen der beiden Grenzfälle mit k,x=k%,=0 oder 1, dann 
stimmen im allgemeinen %k),, und k%, nicht miteinander überein. Nur wenn die Regressions- 
linien geradlinig sind, besitzen die beiden Korrelationsverhältnisse den gleichen Wert, welcher 
dann auch gleich dem Quadrat des Korrelationskoeffizienten r ist. 

Fbenfalls von Pearson wurde für die Messung der Korrelation das Kontingenz- 
ınaß (Mean square Contingeney) 


f- l \ m (a; ya) — mw, (a) m, (ya)y? 
vim — Ina — 1) ,(;) 0, (Yyr) 
ih 
. (8) 
1 | \tr ww) ‚| 
yım In lc w,(x;) w,(Yr) J 
ik 





vorgeschlagen. Hierfür muß eine aritlımetische Verteilung »w (x; x) vorliegen, bei welcher 
m und » die Anzahl der Zeilen und Spalten der sog. Korrelationstabelle ist. Daß bei 
stochastischer Unabhängigkeit wegen w (x; y,) =", (x;) w, (yx), f?=0 wird, ist sofort zu sehen. 
Bei vollständiger Abhängigkeit dagegen, wobei notwendig m =n ist, gilt für die n von Null 
verschiedenen Wahrscheinlichkeiten w (x; yx) 


w (2, y)=mw, (X) = mw, (Yı). 


Daher sind die » nichtverschwindenden Glieder w* (x; y,)/[w, (#;) w, (yx) alle gleich Eins und 
wegen m=n wird auch f?= 

Der besondere Vorzug dieses Pearsonschen Kontingenzmaßes ist der, daß nur die 
Häufigkeiten, aber nicht die &- bzw. y-Werte selbst vorkommen. Daher läßt es sich auch 
dann eindeutig angeben, wenn die Merkmale nicht zahlenmäßig faßbar sind, wie z. B. Farbe, 
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Geschlecht usw. Sein Nachteil besteht darin, daß f? sich nur für arithmetische Verteilungen 
bilden läßt und daß, wie sich später zeigt, eine sinnvolle Verallgemeinerung auf Verteilungs- 
funktionen »w (x, y) überhaupt nicht möglich ist. 

Jedes dieser Korrelationsmaße hat also bestimmte Vorzüge und Mängel; ihr gegenseitiger 
Zusammenhang ist nicht leicht zu übersehen. Im folgenden wird nun gezeigt, daß das 
Korrelationsproblem in seiner allgemeinsten Fassung auf ein Variationsproblem führt, daß 
dieses Variationsproblem weiter auf eine homogene Fredholmsche Integralgleichung und 
damit auf ein Eigenwertproblem zurückgeführt werden kann, und daß die bisher üblichen 
Korrelationsmaße nichts anderes sind als nach verschiedenen Methoden gewonnene erste 
Näherungen für den kleinsten Eigenwert dieser Integralgleichung. Handelt es sich um eine 
arithmetische Häufigkeitsverteilung, dann entartet die Integralgleichung und an ihre Stelle 
tritt ein homogenes, lineares Gleichungensystem. 

Die Frage nach der besten Zuordnung der y und x zueinander tritt auch in der Aus- 
gleichsrechnung auf, wenn es sich darum handelt, eine Reihe Beobachtungspunkte durch einen 
einfachen Kurvenzug möglichst gut wiederzugeben. Die neuen Gesichtspunkte, unter denen 
hier das Korrelationsproblem betrachtet wird, liefern auch eine Übersicht über die Möglich- 
keiten für die Lösung jener Aufgabe und es zeigt sich, daß die für das Korrelationsproble m 
entwickelten Verfahren sich auch, bei der Aufgabe des Kurvenausgleichs mit Vorteil an- 
wenden lassen. 

2. Variationsprobleme und Integralgleichungen. 

Da alle die bisher üblichen Korrelationsmaße bestimmte Mängel aufweisen, stellen wir 
nun ohne Rücksicht auf leichte Berechenbarkeit die Frage nach einem theoretisch voll- 
kommenen und dabei begrifflich einfachen Korrelationsmaß. In bezug auf begriffliche Ein- 
fachheit ist am befriedigendsten der Korrelationskoeffizient r. Sein Hauptmangel ist der, 
daß er gegenüber einer (nichtlinearen) Transformation der Veränderlichen x und y je für 
sich nicht invariant ist. Werden den beiden Merkmalen andere Zahlen zugeordnet, d.h. wird 
das Abszissenmaß x gemäß einer bestimmten Vorschrift f (x) verzerrt und ebenso der Ordinaten- 
maßstab entsprechend einer Funktion 9 (y), dann ändert sich im allgemeinen der Zahlen 
wert für r. Dieser Nachteil der Größe r fällt fort, wenn grundsätzlich alle Transformationen 
dieser Art berücksichtigt werden. Um ein verbessertes Korrelationsmaß zu gewinnen, fragen 
wir daher nach dem größten Wert, den r bei solehen Transformationen annehmen kann und 
kommen so auf das folgende 

Variationsproblem: 

Gegeben ist eine wesentlich positive Funktion w (x, y) mit der Normierung || w(x,J)dedy=1. 
Es sind die Funktionen f(x) und 4 (y) so zu bestimmen, daß 


© ei yay)wia,yJ)dedy ’ 
’ Fate (m Wdrdw = Mas . ..... . 
“2 if? (x) w (a, y)dx d y\ily)w(e, dx dy 
wird. Dann ist K? ein Maß für die Korrelation in der Verteilung »w (x, y). In der Ausdrucks- 


weise der Mechanik bedeutet dieses Korrelationsmaß den größten Wert, den bei allen mög- 
lichen Verzerrungen in «- und y-Richtung der Quotient des quadrierten Zentrifugalmoments 
und der beiden Trägheitsmomente der Verteilung annehmen kann. 
Man kann das Variationsproblem etwas einfacher schreiben, wenn man zur Konkurrenz 
nur normierte Funktionen f(x) und 9(y) zuläßt, welche die Bedingungen 
\Fao)w(e,„dedy=\f’(a)w (e)de=1, | 





x Er (10) 
wi, )dedy=z\fiy)w,(y)dy=1i | 
erfüllen. Wegen der Gleichungen (1) bis (4) sind solche Funktionen z. B. f(x) =1, f(x) = nen 
? b ’ 5 nt 
usw. und gy)=1, g(y) = J g usw. Infolge dieser Festsetzung lautet das Variationsproblem 
der Korrelation in seiner ersten Fassung 
K’=|\f()yg()w(e,Jdedy=Maximum . . . . 2.2... WM. 
Damit dieses Variationsproblem sinnvoll wird, ist jedoch noch eine Nebenbedingung 
erforderlich. Wegen w (x,y) 0 folgt mit Hilfe der Schwarzschen Ungleichung 
felsy)w (a, y)drdyy?=f\\ifle) Yw(z,y)-g(y) Vw (z,y)dedyy 
<s|\fla)w(e, dedy-\\wWw(ae,ydedy=1. 


Bei dieser Abschätzung gilt dann das Gleichheitszeichen, wenn die Funktionen f(x) yw (x, y) 
und g(y)yYw(x,y) zueinander proportional sind. Dies trifft für beliebiges w (x, y) zu, wenn 
f(«)=g(y)=1 ist. Dieser Fall ist aber wegen G]. (1) eine Selbstverständlichkeit und daher 
aus der Betrachtung auszuschließen. Daher sind für das Variationsproblem nur solche 
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Funktionen zur Konkurrenz zuzulassen, die für die „Belegungen“ w, (x) bzw. w,(y) zu f(x)=1 
bzw. zu 9(y)=1 orthogonal sind. Die Bedingungen hierfür lauten 


\f)w(e,)dedy=\f(e)w, («)de—=0,| .. 
2 R (12). 
Wr, y)drdy=|yW)w,()dy=® | 
Sind sie erfüllt, dann ist X? im allgemeinen von Eins verschieden und selbst verständlich 
sicher kleiner als Eins. Insbesondere aber wird dann im Falle vollständiger stochastischer 
Unabhängigkeit wegen w (x, y)= w, (x) w, (y) 
fl) y(y)w (a, y)dedy=|flx) w,(«)de\y(y)w(y)dy=0 
und daher, wie es sein muß, auch K?=0. 
Zwei Funktionen, welche auch die Bedingungen (12) erfüllen, sind nach Gl. 8) z. B. 
, z—a y—b 2 $ ? ’ 
f)= e und yay=}?" ® Für sie geht das Integral der Gl. (11) in das Quadrat des 
Korrelationskoeffizienten r nach Gl. (5) über. Daher gilt für K? die Abschätzung 
Va ee a 
und außerdem wurde damit das obengenannte Ergebnis erhalten, daß der Korrelations- 
koeffizient r absolut stets höchstens gleich Eins ist bei beliebiger Verteilung w (x, y). 
Das Variationsproblem (11) enthält zwei verfügbare Funktionen. Wir führen es nun 
in ein solches mit nur einer freien Funktion über, indem wir uns zunächst f(x) oder 9 (y) 
fest angenommen denken und dann die andere Funktion so hinzu bestimmen, daß das 
Doppelintegral seinen unter diesen Umständen größten Wert annimmt. Diese Aufgabe läßt 
sich mit Hilfe der Schwarzschen Ungleichung (\FGda=”<|F’dx\G?’dx lösen, indem 


F(x)= f(x) ymw, (x), G (2) = Vo, (@) \yyw(ie,y)dy 
gesetzt wird. Damit ist 
\ 1 
\\fa)gy)wie, drdyy?=t\fle)yw, (x) Ve, (e) \ya)wie,y)dydx)’ 
M’ 


— 


zu | 





m 58 ; 
<Z\f’(e)w,(e)dx- \ zen (a)wiz,y)dy’d 
. 1 


wr 
sen ( N h 2 r 
\ Frau E AL] E 200 YAy’de. 


Dies ist der größte Wert, den das Variationsintegral für ein bestimmtes 4(y) und für be- 
liebige f(x) annehmen kann. Die obere Schranke wird dann erreicht, wenn F(x) und @ (x) 
einander proportional sind, wenn also f(x) und 9 (y) miteinander durch die Gleichung 


G (x) \yly)wi(ae y)d y 
2 > % e « arY « =— Pr ei MR | 
F (x) w, (x) f(x) konst. 


verbunden sind. Die Konstante © selbst folgt durch beiderseitige Integration der Gleichung 
0? F? (x<)= @?(x) auf Grund der Bemerkung, daß |\F?(x)dx=1 ist. So ergibt sich für C*? 


@=|@(a)dx —\ 2 m (yy)w(a,y)dy’d, 


also genau der Wert der oberen Schranke in der Abschätzung. 

Wir erhalten damit das Ergebnis: Liegt 9 (y) fest und ist nur die eine Funktion f(x) 
frei verfügbar, dann nimmt das Variationsintegral (11) den unter diesen Umständen größt- 
möglichen Wert 


er 
Je 4 as » 2 . 
( -\,,@lsW@" N 5 
dann an, wenn für f(x) die Funktion 


TOO). (#,y)dy er N a 9. Pe 
eingesetzt wird. 

Wird bei dieser Betrachtung von der richtigen Funktion g (y) ausgegangen, welche dem 
Doppelintegral den größten überhaupt möglichen Wert K? verleiht, dann wird C?= K? und 
daher ist Gl. (14) eine neue Fassung unseres Variationsproblems, wenn man wieder g(y) als 
frei verfügbar betrachtet. Selbstverständlich können bei dieser Zurückführung des Variations- 
problems auf ein solches mit nur einer Unbekannten die Funktionen f(x) und g(y) auch 
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ihre Rollen vertauschen. Wir erhalten also als zweite Fassung für das Variationsproblem 
der Korrelation die folgenden beiden Integrale: 
. 1 3 2 i 
K?=\ — ‚( y4(y) w (x, y) dy)’ de = Maximum 
und j ’ ee en 2 


K’= (\ f(&) w (x, y) d x)” d y= Maximum 


\ w,(4) 
Um diese Integrale näherungsweise auszuwerten, liegt es nahe, für f(x) bzw. g(y) 

lineare Funktionen einzusetzen, die natürlich die Normierungsbedingung (10) und die Ortho- 
: - 2. a . en R Ä 2 —a 

gonalitätsbedingung (12) erfüllen müssen. Diese Funktionen lauten fe) =— bzw. 


Rn E 
g y)=*" r °. Mit ihnen wird wegen (2) und (6) zunächst 


sb . Ya) 
\g(y) m (a, y) d y=\! / wie, Wdy=! = m, (®), 
z(y) — a 


c—a 
\fla)w(z,J9)de= \ zw (z,yde= w,(y) 


Ss 


und daher nach (16) 


ER Be 
> { (x Ir — (2) — D’m.(x z—=k!, 
I (@) gm (e,y)dy’dxe= rE \(y(a) bw, (a)de=k),: 


| 1 

Iw,(y) 

Diese Näherungswerte sind also die beiden Pearsonschen Korrelationsverhältnisse (7). 
Durch dieses Ergebnis wird die Stellung der Korrelationsverhältnisse in der allgemeinen 

Theorie und ihr Verhältnis gegenüber den anderen Korrelationsmaßen aufgezeigt. %A},, und ki, 

sind zwei verschiedene erste Näherungen für unser allgemeines Korrelationsmaß K?. Bei k},.x ist 


A  AEERERE 
(fl) w(a,yda)’dy= = ae) a) w,(y)dy=k%y- 


—t ia . ; 4 s 
uW=", ', während für f(x) beliebige Funktionen in Wettbewerb stehen; umgekehrt ist 


FErRR z—a __ 5 ’ u r ö e 
beik,, f(&)= 2 während eine Mannigfaltigkeit von Funktionen (x) zur Konkurrenz 


zugelassen sind. Beidemal ist also das Variationsproblem nicht voll ausgeschöpft und daher 
müssen im allgemeinen beide Korrelationsverhältnisse kleiner als K? ausfallen, wobei K? ihre 
obere Grenze ist. Anderseits sind sie als Näherungen des Variationsproblems sicher besser 
als der Wert, den das Variationsintegral (11) auf Grund linearer Ansätze für f(x) und 9 (y) 
annimmt. Daher ist r? eine untere Schranke für beide Korrelationsverhältnisse. 

Die zweite Funktion kommt in den Formeln für k),, und k%,„ jeweils nicht mehr explizit 
vor. Sie ist bis auf einen Normierungsfaktor im ersten Fall gleich y(x) — b, im zweiten 
gleich &(y) — a, stimmt also im wesentlichen jedesmal mit der einen Regressionslinie über- 
ein. Auch diese Funktionen erfüllen die Orthogonalitätsbedingungen (12). Sind sie linear, 


‘ ... . IE Emuap y—b 
dann müssen sie wegen der Orthogonalität und der Normierung mit z und * j über- 


einstimmen. Es sind dann die beiden Funktionenpaare miteinander identisch und daher 
liefern auch die beiden Korrelationsverhältnisse den gleichen Zahlenwert, der überdies dann 
mit r? übereinstimmen muß. Damit erhalten wir das oben ohne Beweis genannte Ergebnis, 
daß im Falle geradliniger Regressionslinien Ay, —=k%,=r? ist. 

Denkt man sich in den Variationsproblemen (16) für f(x) und g(y) die Lösungsfunk- 
tionen eingesetzt, so daß die Integrale den Wert K? annehmen, dann bestehen, nach dem 
Vorgang der Gl. (15) für diese beiden Funktionen die simultanen Integralgleichungen 


1 en 
Ki) ed 
(17). 
EN, EEE * 
ıW=, w, gar (@,y)da 


Diese beiden Gleichungen liefern zusammen 
1 (mw(e,»y) "Ifn@y) 1f 
fa =r\ w,(x) 7 (y) d ITXK \ w,(x) . 


w(2,) 

w, (9) 

ze i \ \ w (x, 9y) w (2, y) 
K? )) mw, (2) w, (y) 


f()dzdy 


f(z)dzdy. 
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Dieses Ergebnis ist eine homogene Fredholmsche Integralgleichung für die Funk- 
tion f(x) 
' ‚(W@,z) 
f@=il se f(z)d2. ET A EEE TE EEE 
mit dem Kern 
4 mw (a, y) w (2, y) 
We2)=\ 2 


Man kann diese Gleichung auch als solche mit symmetrischem Kern hinschreiben. Sie 
lautet dann 


EA: 1 
dy und mit /= Re: 


;\ W(x 
f (x) Ymw, (©) = A| 27, Er 


Selbstverständlich gibt es auch für 9 (y) eine entsprechende Integralgleichung. Sie kann 
durch Elimination von f(x) aus den beiden simultanen Integralgleichungen (17) gewonnen 
werden. 

Die erste, selbstverständliche Eigenfunktion der Integralgleichung (18) ist f(«)=1, 
der zugehörige Eigenwert A=1. Durch die Nebenbedingung (12) für f(x) ist sie bereits aus 
der Betrachtung ausgeschlossen worden. Wir spalten sie daher folgerichtig vom Kern ab 
und erhalten dann mit dem am besten bei Gl. (18a) reduzierten Kern 


NO PS © 


W (x, 2 1 : s W (x,2) —w, (x) w, (2) 
\ -— Yw, (2) f(x) yw, (2) f (z) : 
yw()w,(e) 4 vrareym (ar yw,(x) m, (z) 
nach Beseitigung der Wurzeln die Integralgleichung 


"We, B a (x) w, (2), 


fi) =A\ Ber er De 
,(®) 
Die Frage nach dem durch das a (11) erklärten Korrelationsmaß K? 
ist wegen A=1/K? nunmehr die Frage nach dem kleinsten Eigenwert dieser Integralgleichung. 





Dieser erste Eigenwert ist nach der Theorie der Integralgleichungen auch äurch ein Variations- 
problem bestimmt und daher ergibt sich nun für das Variationsproblem der Korrelation die 
folgende dritte Fassung 


’=- =|\\XW(a,2) -w,(x)w,(2)} f(x) f()dedz= Maximum . . . . (20). 
“min 
Das Maximum entsteht, wenn für f(x) die erste Eigenfunktion der Gl. (19) eingesetzt wird. 
Verwendet man statt dessen für eine Näherungsrechnung die Funktion f (x) =, dem 
erhält man wiederum das Pearsonsche ee kz y. 

Ohne weitere Annahmen über die Verteilungsfunktion w (x, ,) läßt sich über die Inte- 
gralgleichung natürlich nicht sehr viel aussagen. Vor allem . ihr Kern sicher definit, 
denn wenn man bei (20) die Orthogonalitätsbedingung | f(x) w, («) d=—=0 berücksichtigt und 
die Reihenfolge der Integrationen vertauscht, entsteht wieder der Ausdruck (16) für K?, und 
diesem Variationsintegral ist sofort anzusehen, daß es nur positive Werte annehmen kann. 
Dagegen kann man nicht damit rechnen, daß die Eigenwerte einfach sind, und daher ver- 
sprechen Verfahren zur Ermittlung des ersten Eigenwerts, die auf eine Trennung der Eigen- 
werte hinauslaufen, hier keinen Erfolg. Um eine Übersicht über die bestehenden Möglich- 
keiten zu gewinnen, betrachten wir nun wieder die beiden Grenzfälle vollständiger Unab- 
hängigkeit und vollständiger Abhängigkeit der beiden Merkmale. 


a) Vollkommene Unabhängigkeit: 


In diesem Fall wird wegen w (x,y)=w, (x) mw,(y) das Integral auf der rechten Seite 
der Gl. (19) für eine beliebige Funktion f (x) 


| w en m, (x) w, @, „ie (x, y) w 2) e)dyda (we) fle)de 


(@) (2) x c) mw, (9) 
—=||w, a w,(y)f(@)dydz -\w,(2)f(e)dz 
w,(@)f(@)d2— \w,()f(2)d2=0 
Daher wachsen für alle als Eigenfunktionen der Integralgleichung in Betracht kommenden 
f(x) die zugehörigen Eigenwerte 4 über alle Grenzen. Vor allem ist bereits der kleinste 
Eigenwert unendlich groß und daher AK?’= ; - =0. 


"min 
24 
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Während aber aus dem Verschwinden des Korrelationskoeffizienten r und auch der 
Korrelationsverhältnisse ky,, bzw. k&, nicht folgt, daß x und y voneinander stochastisch un- 
abhängig sind, ist dieser Schluß berechtigt, wenn das neue Korrelationsmaß A? den Wert 
Null annimmt. Dies folgt am einfachsten durch Betrachtung der Bilinearreihe für den Kern 
der Integralgleichung. Sind f;(&) die sämtlichen Eigenfunktionen und 4; die zugehörigen 
Eigenwerte, die durchaus nicht alle voneinander verschieden sein müssen, dann lautet diese 
Darstellung 

W(x,2) —w, (x) w, (2 ER. 
4 @) he r ‚N = L yw()hk(a)ywl@)fi(). - - - -. (21). 


Im vorliegenden Fall sind alle 4, unendlich, und daher vereinfacht sich (21) zu W (x, z) 
w, (x) w,(2)=0 und weiter für z=x zu 


Wi ee EN I 


Mit Hilfe der Definitionsgleichung für W (#,2) und der GIn, (2) wird nun 


W (2,0) -nt(a)=\T Y) 


w,(4) 
\ era) 
\ W, (y) 


dy— 2m, (x) \w(a,y)dy+ nt (a)\w,(y)dy 


2m, (2) w (x, y) +} (x) m, 0) dy 
—— 1 Hin . . . . \2 
\ w,@) ve (a) er, le), drAay—d. 
Da der Integrand wesentlich positiv ist, kann diese Gleichung nur dann bestehen, wenn 
überall dort, wo w(x,y) und damit auch w, (x) und w,(y) von Null verschieden ist, die 
Beziehung 
vw (2,y)= mw, (x) w, (y) 
gilt. Das heißt aber, die Merkmale « und y sind in der Tat stochastisch unabhängig von- 
einander. 


b) Vollkommene Abhängigkeit: 

In diesem Fall treten nur längs einer Kurve y=y(x) in der xy-Ebene von Null ver- 
schiedene Wahrscheinlichkeiten auf. Die zweiparametrige Verteilung » (x, y) ist als Funk- 
tion von & und y ausgeartet, denn es ist 


[0 für y 3 y (x) 
w (xy) = : 
Io für y=y(«). 
Dagegen sind die Integrale 

\w(z,yJ)dy=w,(x) und \w(z,y)de=m,(y) 
definiert, und es ist wie stets \w, («)de=|m,(yJ)dy=1. 

Wir nehmen nun an, die Zuordnung zwischen den zusammengehörigen «- und y-Werten 
sei umkehrbar eindeutig, so daß y=y (x) und ebenso die Umkehrung z<=. (y) eindeutige 
Funktionen sind. Wegen der ausschließlichen Verteilung der Wahrscheinlichkeiten längs 
der Kurve y(x) besteht dann die Differentiationsformel 

w,()de=m,(y)dy, 
und außerdem ist für jede beliebige, eindeutige Funktion f (x) 
na )dr=f(ey) wa )dr=flely)-w,(y). 
Diese Tatsache ermöglicht die Berechnung des Integrals auf der rechten Seite der Gl. (19) 
für eine beliebige Funktion f(x), welche die Orthogonalitätsbedingungen (12) erfüllt. Es ist 
"W(e,2) — w, (x) w, (2) "le (a, y) w (2, y) 
z2)d2= — — ——f(z2)dydz— \w, (2) f(2)dz 
\ w, (e) rade-\\ am, [Air Imre) 


FAR, 
nl) 


ron @naz- me nay 


or 
=..@ End 
=f(z(y(e))) =f(@). 


In diesem Fall sind also alle zu w, («) orthogonalen Funktionen Eigenfunktionen der 
Integralgleichung oder Linearkombinationen solcher. Die zugehörigen Eigenwerte aber sind 
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alle gleich Eins, d.h. 4=1 ist unendlichfacher Eigenwert. Eine normierte Eigenfunktion 
ist z. B. 
a 


fd=——. 


Um vollständige Abhängigkeit nachzuprüfen, genügt es daher, mit diesem einfachsten Ansatz 


zu rechnen, was auf die Ermittlung des Korrelationsverhältnisses k% „hinausläuft. Jedenfalls ist 
also zum Nachweis der Korrelation Eins die Lösung der Integralgleichung nicht erforderlich. 

Es ist nun noch die Frage zu klären, ob auch umgekehrt beim Vorliegen des Korre- 
lationsmaßes K’=1 mit Sicherheit auf vollständige stochastische Abhängigkeit der beiden 
Merkmale, d. h. auf eineindeutigen, funktionsmäßigen Zusammenhang zwischen x und y ge- 
schlossen werden kann. Im angenommenen Fall gibt es eine Funktion f(x), welche zu- 
sammen mit der nach Gl. (17) zugeordneten Funktion g (y) dem Variationsintegral für K? den 
Wert Eins verleiht. Wir betrachten nun den Ausdruck für AK? in seiner ursprünglichen 
Fassung nach Gl. (11). Wegen der Normierungen (10) ist allgemein 


K’=\\f()y We (a, drdy=z\\tf’(a)+ Pl) — If) - gr we (ae, )drdy 
ih u 
=1-Z\\f@&) - WW) vw @p)drdy. 
Damit sich X’=1 ergibt, muß also 


IV) -sWfw a dzdy=0. . ... ...% 0.1) 


sein. Dies ist nur dann der Fall, wenn überall dort, wo w(x#,,) nicht verschwindet, f (x) 
— 9 (y) = const ist. 

Umgekehrt ergibt sich K’=1 immer dann, wenn »(x,y) nur in solchen rechteckigen 
Bereichen von Null verschieden ist, wo Streifen parallel zur y-Achse und zugeordnete 
Streifen parallel zur x-Achse sich überdecken. Zu allen Streifen eines einzelnen solchen 
Streifensystems (Bild 1) gehört der gleiche Wert der Konstanten für f(«)=9g(y). Die Funk- 
tionen f(x) und g(y) können, weil sie orthogonal zu w, (x) bzw. w,(y) sein müssen, nicht 
für alle x und y konstant sein, sondern sie sind stückweise konstante Treppenkurven, die 
mindestens je einen positiven und einen negativen Wert annehmen. Daher besteht der Ver- 
teilungsplan in der zy-Ebene aus mindestens zwei solchen beschriebenen Streifensystemen. 
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Bild 1. Einzelnes Weil 
Streifensystem. Bild 2a. Bild 2b, Bild 2e. 


Die Bilder 2a bis e zeigen «rei Beispiele für derartige Verteilungspläne. Nur in deu 
schraffierten Feldern, die zu waagerechten und senkrechten Streifen mit gleichen Ziffern ge- 
hören, darf w («,y) von Null verschieden sein. Wird die Streifeneinteilung unbegrenzt ver- 
feinert, dann entsteht der oben betrachtete Fall, daß nur längs einer Kurve der x y-Ebene 
die Wahrscheinlichkeiten nicht verschwinden. Diese Kurve braucht aber nicht, wie oben für 
die Rechnung angenommen wurde, die x- und y-Werte einander eineindeutig zuordnen, denn nur 
beim ersten Beispiel ist dies der Fall, während die anderen beiden Beispiele auf mehrdeutige 
Funktionen führen. Es ist für die Verteilungspläne kennzeichnend, daß nicht punktweise 
die x und y einander zugeordnet werden, sondern daß Klasseneinteilungen für die Abszissen 
und ÖOrdinaten auftreten und eine eineindeutige Zuordnung nur zwischen den x-Klassen und 
den y-Klassen mit gleichen Ziffern besteht. Für den stochastischen Zusammenhang von & 
und y entsteht dadurch ein gewisser Spielraum, den auch das Korrelationsmaß A? nicht zu 
erkennen gibt. 

Das Korrelationsmaß A? ist gemäß seiner Definition unabhängig von der Zuordnung 
der x- und y-Werte zu den beiden Merkmalen der Verteilung. Bei den bisher entwickelten 
Formeln für K? kommen aber die beiden Zuordnungsfunktionen f(x) und 9(y) einzeln oder 
zusammen explizit vor. Besonders wünschenswert wäre dagegen ein Ausdruck für den 
kleinsten Eigenwert der Integralgleichung (19), welcher weder f(x) noch 9 (y), sondern nur 


24° 
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die Verteilungsfunktion ı»* (x, y) enthält. Um einen solehen Ausdruck zu gewinnen, liegt es 
nahe, von der Entwicklung des Kerns durch eine Bilinearreihe nach Gl. (21) auszugehen, 


W (x, 2 m, (x) mw, 2) _ \ . w,(& )f; (x) j m, (2) fi 3 Go BE - ; ; \ 


vyr,(ae)w, (2) — 


Wir betrachten diese Gleichung nur formal, ohne nach ihrer Konvergenz und nach den dazu 
notwendigen Voraussetzungen zu fragen. Für 2=.r liefert die Integration über & 


Me It 


Der Integrand erlaubt folgende beiden Umformungen 


W(e,2) — m; (a) [ w’(e,y) ua. wi „) (2) w, Di 
w, (x) = ) m, PL, ” w,()=\ Er (y) 


Daher folgt für J das rg 


J= ea Ei nr. Wr dedy=|\\ ni} dedy Pr a: (24). 
3. ı (8) w, (9) we, (@) m, (y) h hr 

Dieser Ausdruck stimmt in seinem Aufbau mit dem Pearsonschen Kontingenzmaß (8) 
für arithmetische Verteilungen überein. Er wäre vielleicht geeignet, um Schlüsse auf den 
ersten Eigenwert 4, und damit auf die Stärke der Korrelation zu ziehen, wenn dieser Eigen- 
wert sicher einfach wäre und wenn die Reihe der reziproken Eigenwerte konvergieren würde. 
Beides muß nicht der Fall sein, denn wir haben gesehen, daß im Falle vollkommener Ab- 
hängigkeit 4,—=1 sogar ein unendlichfacher Eigenwert ist, was natürlich auch das Unend- 
lichwerden des Kerns der Integralgleichung zur Folge hat. Man hat also beim Korrelations- 
problem sehr mit singulärem Verhalten der Integralgleichung zu rechnen, und daher lassen 
sich die aus der Theorie der Integralgleichungen bekannten Verfahren zum großen Teil nicht 
anwenden. Ein klares, aber selbstverständliches Ergebnis läßt sich nur im Falle vollständiger 
Unabhängigkeit von « und y feststellen. In diesem Fall ist bereits /, unendlich groß und 
ebenso alle übrigen Eigenwerte; daher ist nach (24) J=0. 

Während bei geometrischer Verteilung w (x, y) die Größe J sich demnach nicht als 
brauchbar zur Messung der Korrelation erweist, gibt sie für arithmetische Verteilungen') 
gewisse Aufschlüsse über die Korrelationsverhältnisse. Bei arithmetischen Verteilungen 
treten an die Stelle der Integralgleichungen lineare Gleichungensysteme, und die Anzahl der 
Eigenwerte ist endlich und gleich der Zeilen bzw. Spalten der Korrelationstabelle. 

Ist »w (&; 9x) die Wahrscheinlichkeitsverteilung für m mögliche «-Werte x; und » mög- 
liche y-Werte y;, dann lautet das Gegenstück zur Kernfunktion W(x,z) der Integral- 
gleichung (18) 

. 
Ya e) u mit TIER, 


pe 0, (Mr) i 
k 


W (ai a) = 


und der Integralgleichung selbst entspricht das lineare Gleichungensystem mit m Gleichungen 





OR ac, Ya uni. -  . \ w ri TIYR) 2;. 
f(&;) " w, (z;) f(&ı) =A — , v,(; ) wo, (ya) f(«&ı). 


Die Spur der Determinanten dieses Gleichungensystems liefert die Summe na m reziproken 
Eigenwerte 


Er »? (2; %% 1 | l 
Da ir 
(2) m; (Yyr) ko hy Am-ı 


Dabei ist 4,—=1 der selbstverständliche Eigenwert, der zur Eigenfunktion f(x) =1 gehört 
und beim Übergang von Gl. (18) zu Gl. (19) abgespalten wurde. Wir erhalten damit das 
Ergebnis 


I > - Bd) _5_ tw lin) —mı (a) m, (yn)y? 
PAR we ) m, (yr) ds w,(2;) w, (yxr) 
ik ik 
(24a). 
1 
+ we: 
*n- 


1) Für aritnmetische Verteilungen läßt sich die hier vorgetragene Korrelationstheorie mit elementaren alge- 
braischen Hilfsmitteln entwickeln. Sie ist ausführlich dargestellt in einem demnächst erscheinenden Buch des Ver- 
fassers über mathematische Statistik. (Verlag Quelle u. Meyer, Hochsehulwissen in Einzeldarstellungen.) 
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Ebenso könnte man die Rechnung auch mit Hilfe eines linearen Gleichungensystems für 
die 9 (yx) durchführen, welches der Integralgleichung für g (y) entsprechen und » Gleichungen 
enthalten würde. Das Ergebnis (24a) würde dadurch sich nur insofern ändern, daß die 
Reihe der reziproken Eigenwerte mit dem Glied 1/A„_, abbrechen würde. Da aber der Aus- 
druck auf der linken Seite in beiden Fällen übereinstimmt, weil er in bezug auf die x und y 
symmetrisch ist, und da, wie sich leicht zeigen läßt, die beiden Darstellungsweisen auf 
identische Systeme endlicher Eigenwerte führen’), so müssen, falls m z£n ist, alle nur in 
der einen Entwicklung vorkommenden, überschüssigen Eigenwerte unendlich sein. 

Es sei nun n<m. Die Zahl der von Null verschiedenen reziproken Eigenwerte ist 
dann n—1. Der kleinste dieser Eigenwerte, A,, liefert unser Korrelationsmaß K’—=1/},, für 
alle übrigen aber gilt 1/4, =K?’. Also ist J<S(n—1)-K? und erst recht, wenn eine in m 
und » symmetrische Abschätzungsformel bevorzugt wird, 

J 
y(m —- 1) (n —1)' 
Der Ausdruck auf der rechten Seite dieser Formel ist das Pearsonsche Kontingenz- 


2} 


maß f? nach Gl. (8). Wir erhalten damit das Ergebnis 


J<ym—-Dm-IK: oder K° 


= 1 \ tm (a ya) — mw, (a) m, (yr)y” 
Van —-YDn—-D_ w, (x) w, (Yx) 
ik 
2 1 | \Y _ "a _,\ - 
 VYlm —1)(n—1) — w, (2) ©, (Yr) J EEE a. 
l FR 
Y(m Banictt R 





Es ist also f? sicher kleiner als unser Korrelationsmaß K°? und damit kleiner als Eins. 
Seinem Aufbau nach ist f? ein Mittelwert aus den (m — 1) bzw. (n—1) reziproken Eigen- 
werten des aus der Korrelationsaufgabe entspringenden Eigenwertproblems. 

Diese Zusammenhänge lassen erkennen, daß erstens eine Anwendung des Pearsonschen 
Kontingenzmaßes auf geometrische Verteilungen grundsätzlich nicht möglich ist, denn dieser 
Mittelwert der reziproken Eigenwerte ist bei einer Integralgleichung entweder unbestimmt 
oder Null, also keine für die Messung der Korrelation geeignete Größe; zweitens, daß das 
Kontingenzmaß mit dem Korrelationskoeffizienten r und den beiden Korrelationsverhält- 
nissen Ay, und kyx Sich nicht gut vergleichen läßt, weil der Zusammenhang jener Größen 
mit K? von dem in Gl. (25) dargestellten allzu verschieden ist. 

Besonders bemerkenswert ist wiederum der Fall vollkommener stochastischer Abhängig- 
keit. In diesem Fall ist m=n, und für die allein vorkommenden, zusammengehörigen 
Werte x; und 9%; gilt 

"(y)ew, ()=R, (Yyr): 


Das Gleichungensystem, das an die Stelle der Integralgleichung getreten ist, reduziert sich 
unter diesen Umständen auf die n Gleichungen f(#;)=4f(«;), welche alle die Lösung /=1 
besitzen. Daher ist A=1 n-facher Eigenwert des vollständigen Kerns und damit (rn — 1) 
facher Eigenwert des reduzierten Kerns. Es wird J=n—1 und f?=1. Das Kontingenzmaß 
läßt also bei arithmetischen Verteilungen vollständige, stochastische Abhängigkeit sicher 
erkennen. 


3. Reihenentwicklungen für die Größen der Korrelationstheorie. 

Die oben entwickelte Korrelationstheorie ist nichts anderes als eine allgemeine Theorie 

der zweidimensionalen Verteilung w (x, y). Um die Funktion ® (x, y) nach Orthiogonalfunktionen 
zu entwickeln, braucht man zwei Funktionensysteme, nämlich eines mit Funktionen von x und 

eines mit Funktionen von y allein. Wir verwenden zu diesem Zweck die normierten 

ÖOrthogonalpolynome steigenden Grades 9„(x) zur Belegung m, («)—=\w(x,y)dy und y, (4) 
zur Belegung w,(yJ)=\w(a,y)d&. Für diese Polynome bestehen die Gleichungen 

(0 für iS£j _j0 für dZ£j 


\p; la) Pla) w, (x) de= ein ; und \ vl) yzly) mw, (y) d y= 1 ein (26). 


Die ersten dieser Polynome lauten 
x l 


a Yy ) 
9)=1; p,(e) = — usw, »y)=l; y,(yW= z usw. 


2) Der Beweis findet sich im Buch des Verfassers. 
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Für weitgehend beliebige Funktionen f(x) bzw. 4(4y) bestehen die Reihenentwicklungen 
ao 
I@)= ee % a; 9; (x) mit ;=\fla)gp;le)w,(e)d« 
i—0 
Es 
Y . 
smw= 2 by) mit eig) yiinm,(Wdy 
j=0 
Wir wollen im folgenden annehmen, daß die vorkommenden Reihen dieser Art gleichmäßig 
konvergieren, wenn sie nicht sogar nach endlich vielen Gliedern abbrechen, so daß stets 
gliedweise Integration gestattet ist. Unter dieser Voraussetzung ist z. B. 


\ f’ (x) w,(«)de= \ | ap; («)) w,()de= Rn a; a; \ Pilz) ya (a) w, (x) da = .E a; 
J\>B i,k=\ i=( 


und entsprechend für g(y), so daß für f(x) und g(y) die sog. Vollständigkeitsrelationen gelten 


L an 
if (&)w, ()de= Na;; \eiy)a,(y)dy= 3 b; RS Pr 
i—lV j=V 


Für die Wahrscheinlichkeitsverteilung (x, y) bestehe eine (ebenfalls gleichmäßig kon- 
vergente) Entwicklung der Form 
L 
\r 
w(a,y)= z GjPila)yv;ly)w, (ae)w,(y). 
i,j=0 
Durch Multiplikation beider Seiten mit dem Produkt zweier bestimmter Polynome 9;(x) und 
y;(y) und Integration über x und y ergibt sich dann wegen der Orthogonalität der y; (x) 
und der y;(y) für die Entwicklungskoeffizienten e;; die Formel 


elle) y;p)w(e,y)daedy. 
Diejenigen Koeffizienten e;;, für welche ö=0 oder j=0 ist, lassen sich sofort angeben. Es ist 
für i=j=0 elle, ydedy=l, 
füri=0, j7V0 ey; ()wle,y)dedy=\y;ly)a,(y)dy=V 


und 
für 1320, j=0 c=l\lpile)a (a, y)dedy=\gp;(r)w,(e)de—=0. 


Damit vereinfacht sich die Reihenentwicklung für »w (x, y) zu 
Io) 
(2, y=m,(z)w,(y)+ . CGz;pila) yily)w, (x) w,(y) 
ij 
oder 
1 + he} ey) mt eigeyz)wle,ydedy . . (29. 
ij 


wir, My) 
w,(e)m,(y) 
Wir vermerken noch die der Vollständigkeitsrelation (28) entsprechende Formel für 


mix, 9). Sie lautet 


vw? (se, 9) 4/2 = 
Nom, **ds- \\( 2 eier) wa), (y)drdy 


iR it 
on ee) 
Y 2 2 2 Y 2 
= pe. ee y)w,(y)dy= = Ci; 
i,j=0 ij N 


oder, weil e,,=0 und Gi =V ist für 5 bzw. j0, 


BE hy ie 


i,j=1 


\ \ m’ (a, 9) 
w,(z)w,(y) 


Dies ist der Ausdruck J von Gl. (24), der dem Pearsonschen Kontingenzmaß f? nach 
Gl. (25) für arithmetische Verteilungen zugrunde liegt. 

Es ist noch zu erwähnen, daß im Falle einer arithmetischen Wahrscheinlichkeitsver- 
teilung w (#; ya) mit endlich vielen Zeilen und Spalten der Korrelationstabelle die Anzahl der 
Funktionen 9; (x) mit der der Spalten und die Anzahl der Funktionen y;(y) mit der der 
Zeilen übereinstimmt. Die betrachteten Reihenentwicklungen brechen daher in diesem Falle 
alle mit endlich vielen Gliedern ab. 
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Im Mittelpunkt unserer Korrelationstheorie steht der Ausdruck (11) 
K'=| \ \f(«) g(y)w(x,y)dxd yy- - Maximum, 
wobei die zur Konkurrenz zugelassenen Funktionen f(x) und g(y) noch die Nebenbedingungen 
\f(a)w, («)de=0, \F(ia)w(s)de=1, 
\yy)w,Wy)dy=®, \W)r,J)dy=1 


erfüllen müssen. Nach (28) und (29) besagen diese Orthogonalitäts- und Normierungs- 
bedingungen in bezug auf die Reihenentwicklungen für f(x) und g(y), daß 


n 2 
u.=b,=0 und da= N 1 ER 2 2 RE 
R s i—1 j1 
sein muß. Nun wird 
\ \fa)g()wie,yJ)dxdy 
u be) n 
||) a;y;(2) 2 b; vi) L Ckı Pr) yvıly)w, ()n,(J)dazdy 
il] j=!1 k,l=0 
ee) 
= Z 9b;0;;- 
ijel 


Das Variationsproblem der Korrelation lautet also in der neuen Fassung mit Hilfe der Ent- 
wicklungskoeffizienten e;; 


er —r w . 
RK? =| R.; a;b; Cij Maximum 


ij 
L ee 
. I ß WR, or VE 
unter den Nebenbedingungen 2, a;=1 und Ze b;=1. 
i=1 j=]1 
Der einfachste Näherungswert für K*? entsteht, wenn man a,—=b,=1 und alle übrigen 


a; und 5; gleich Null setzt, d.h. wenn man mit den Prüffunktionen f(&)=g, (x) = s u und 
yy)=y,(y)= y j : rechnet. Das Ergebnis dieser Näherungsrechnung ist an Stelle von K? das 
(Juadrat des Korrelationskoeffizienten r. Also ist nach (32) 

he VE Fe ee er Tr En  _° 


In der zweiten Fassung des Variationsproblems (16) und in den Integralgleichungen der 
Korrelation (17) kommen die Funktionen 


\fi«) w(z,9)dz und \gap)wi(a,y)dy 


vor. Mit Hilfe der Reihenentwicklungen wird 


\ fia)w(e,J)de= \ = a; Pr (x) i + e% G;Ypila)wy; u w, (2) w,(y)da 
ki ijj=]l 
— N a;6;y;(y) w, (y) 
ij! 
und ebenso 
\sWw iz, y)dy= 2 bc ;Wi(@)w, (@). 
ij] 


Die Integralgleichungen sagen aus, daß diese Funktionen mit 


Ky(y)n,(y)=Kw,(y Lb;y;(y) 
j—=1 
und mit : 
Kf(&) w, («)= Km, (x) N a;9;(«) 
i=1 











. ng y : Z. angew. Math. Mech. 
376 Gebelein, Das statistische Problem der Korrelation Bd.9r Nr.6 Dez. 1941 





übereinstimmen. Durch Vergleich ergeben sich daher die folgenden, an die Stelle der simul- 
tanen Integralgleichungen tretenden, unendlichen linearen Gleichungensysteme für die Koeffi- 
zienten a; und b;: 


1 — N = 
d; m \ b;c;;, = \ d,Cij e > R i ° . r R z (34). 
31 I3=1 


Die Frage ist die nach dem größtmöglichen Wert für K. Man kann K isolieren, indem man 
beide Seiten dieser Gleichungen quadriert, addiert und die Vollständigkeitsrelationen (31) be- 
achtet. Auf «diese Weise ergibt sich 


x -: \?2 ) 
K’= #2 | z b; a Maximum, 
i=1 el 
‚ (35). 
K’= 3 | v 


a; ci;) —= Maximum 
jet \i=1 


Dies sind die beiden Variationsintegrale in Koeffizientenfassung, wie man auch durch un- 
mittelbares Ausrechnen der beiden Ausdrücke (16) sofort bestätigen kann. 


Setzt man in der ersten dieser Formeln 4, =1,b,=b,=-:—0, dann ergibt sich das 
Pearsonsche Korrelationsverhältnis %k)., und ebenso ergibt sich k%,, wenn man in die 
zweite Formel a,=1 und a,=a,=--=0 einsetzt. Wir erhalten damit für die beiden Korre- 
lationsverhältnisse die Darstellungen 

nn be) 
TO, 20 TESERR, 3 or 
kyx= —_— GC ku=. Ci; . . . . . . . . . . . (36) 
i-1 jm1 


Man ersieht an diesem Ergebnis sofort, daß beide Korrelationsverhältnisse größer als r?=e}, 
ausfallen müssen, und man überzeugt sich auch leicht, daß sie nur dann untereinander und 
mit r? übereinstimmen, wenn die oben besprochenen Regressionslinien geradlinig sind. 
Für die Regressionslinien 4(x) und = (y) gelten nämlich, wie oben im Anschluß an Gl. (16) 
ausgeführt wurde, die Gleichungen 

y (a) —b 4 h l 

Re = randye, Rdn nile), 


i=!| 
z(y) —a IE ER EET WEHT. RR CE  - 
Fu) ehe er ade Zu. 


’ => 
Sie sind demnach dann und nur dann geradlinig, wenn die Reihen rechts sich auf die linearen 
Funktionen e,,9, (x) bzw. e,,y, (y) reduzieren. Die Bedingungen dafür aber lauten 


ur ei hauen... . ER, 


und dies sind die gleichen Bedingungen, unter denen ky.=k%,=r? ist. 

Durch Elimination der Funktion g (y) aus den beiden simultanen Integralgleichungen (17) 
wurde oben die Fredholmsche Integralgleichung (18) für f(x) allein gewonnen. Ebenso 
entsteht durch Elimination der b,; aus den Gln. (34) ein lineares Gleichungensystem für die 
Größen a; allein. Es ist 





0. ne L L 
De a S I 4 $ 
2 TEE EEE ' . ' A 
H=R bj ckj= KR: Heim e$ Cr mit Cm I Gjchj- 
oc N man il 
j=]1 j=1 i=l i=] we 
Das der Integralgleichung (18) entsprechende Gleichungensystem lautet also 
u 
1 a ” 
(6. = K: \ a4;C;; mit Cr d- 2 C;jChj (k Te ARE: a an 
— j 1 


i=1 
Für das Maß der Korrelation kennzeichnend ist der größte Wert K?. Er ist zu er- 
mitteln aus der „charakteristischen Gleichung“ 


C, K? C. Cs 
0, GC; K:’ Us 0 39 
N ans 5 
( 31 ( 32 ( 33 K 
Eine erste Näherung für K? ist K’=0,= aj;=kz,. Dieser Wert ist nichts anderes als 
FE 
das auf Grund von f(x) =, (x) sich ergebende Korrelationsverhältnis. 
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Bemerkenswert ist noch die Tatsache, daß die Größen C;, auch ohne Zurückgreifen 
auf die Entwicklungskoeffizienten c;; berechnet werden können. Wegen der Entwicklung 
ee) 
” v , 
\f(#) w (x,y) de= R. a; C;;y;(y) w,(y) 
| ij=1 
ist nämlich mit f(x) = 9; (x) bzw. f(x) = yı; (x) 


2 


\gila)w (e,y)de= F Gjy;(y)w,(y) und \pr (er) w (ae, yde= P.> CHı yıly) w,(y). 

Daher ist 

Bi | 

\ wg) \p9:(a)w (a, y)da \or()w (a, y)dx-dy 

=| 2 Gjyzly) N crıyı W)e,)dy= 3 Gjekj=Cik: 
jel 1-1 j= 
Das Ergebnis lautet also 
Y { 1 
C;rk= \ w,(y) \pi: (a) w (x, y)dx \ Pr (a) w (a, 9) he RE 


Dem Aufbau nach sind diese Größen Verallgemeinerungen des Ausdrucks für das Korrelations- 
verhältnis k%,. Dieses ergibt sich für i=k=1. 

Die Ergebnisse dieses Abschnitts lassen erkennen, daß alle Größen der Korrelations- 
theorie sich ziemlich einfach mit Hilfe der Entwicklungskoeffizienten e;; der Verteilung w (x, ») 
ausdrücken lassen, und daß diese Reihenentwicklung es gestattet, die Korrelationseigen- 
schaften gut zu übersehen. Das Korrelationsproblem wird damit auf das Problem der Reihen- 
entwicklung von w(x,y) nach den zu den Verteilungen w, (x) und w,(y) gehörigen Ortho- 
gonalpolynomen zurückgeführt. Es dürfte sich lohnen, diesem allgemeinen Entwicklungs- 
problem Aufmerksamkeit zu schenken und für seine praktische Durchführung Rechenver- 
fahren zu entwickeln, welche die Lösung mit erträglichem Zeitaufwand ermöglichen. 

Die hier dargestellten Überlegungen beziehen sich alle auf die Beurteilung einer vor- 
gelegten Verteilung w(x,,y). Handelt es sich darum, die Ergebnisse auf umfassendere 
statistische Gesamtheiten zu übertragen, dann sagt die Wahrscheinlichkeitsrechnung aus, wie 
sehr damit zu rechnen ist, daß die Verteilung w (x, y) bei Erweiterung des Beobachtungs- 
materials sich ändert. Die Folgen für den Korrelationsbefund lassen sich mit Hilfe des neuen 
Korrelationsmaßes K infolge seiner Extremumseigenschaft besonders gut beurteilen, wie im 
genannten Buch des Verfassers gezeigt wird (Ergebnis von Teil E Gl. (50)). In diesem Werk 
wird auch die Theorie in ihren Einzelheiten an vollständig durchgerechneten Beispielen 
(Beispiel 34 und 35) erläutert, so daß hier von der Wiedergabe von Beispielen abgesehen 
werden kann. 


4, Zusammenhang zwischen Korrelationstheorie und Ausgleichsrechnung. 


Dadurch, daß in der hier entwickelten Korrelationstlieorie «ie Zuordnungsfunktionen f(x) 
und g9(y) in jedem Fall erst so bestimmt werden, daß der stochastische Zusammenhang 
zwischen x und y möglichst klar zutage tritt, besteht ein enger Zusammenhang mit der be- 
kannten Aufgabe der Ausgleichsrechnung, eine Reihe von Beobachtungspunkten durch einen ein- 
fachen, glatten Kurvenzug möglichst gut wiederzugeben, Bild 3. 

Die Beobachtungspunkte mit den jeweils zusammen- ty Y=Flx) 
gehörigen Koordinaten x und y stellen, unter dem Gesichts- 
punkt der Statistik betrachtet, eine zweidimensionale stati- 
stische Gesamtheit mit der (arithmetischen) Häufigkeitsver- 
teilung w (x, y) dar. Der stochastische Zusammenhang zwischen 
x und y muß sehr eng sein, wenn Kurvenausgleich überhaupt [_ 
Sinn haben soll; das Maß der Korrelation ist also nahezu A 
gleich Eins. Die Aufgabe der Ausgleichsrechnung ist die, eine Bild 5.2 Amen einer 
Kurve mit der Gleichung ; 








—> I 


„=F«(«) 
so zwischen den Beobachtungspunkten hindurchzulegen, daß das mittlere Abweichungsquadrat 


9=|\(y—- Fe)’ w(e.Jdedy=Minimum . . . .... (4) 
wird. 
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Wir arbeiten wieder mit den Orthogonalpolynomen g,(2) und yw;(y). Demnach ist 


ı l 
wegen y, (1) ne: umgekehrt y=b+ty,(y), und für F(x) machen wir den Ansatz 


t 

F (x) - b+d,+d,9,()+d,9,(x)+ :- ge +2 d; 7;(&). 
i—l 

Damit wird 


nz 2 


9-\\ (b Hm, Wb- N age) wa Wdzdy 
2 i=0 


E. 


\ \ (‘ v, (m) Be % d;y; @) vw (z,y)dedy 


s—=® 


Je 
Es 


=’ \ | v(y)w(a,y)dedy-- 2 \ \yv, 2 d;pi(«) w (x, J)daedy 
i—V 


\2 


F \ \ | > d; pi (@)) w(e,y)dedy 
Jo i—=\ } 
Pu, Wiy-2t I die) do yWdrdy 
i=0 


+\ | I dig; (o)) w, (x) dx 


i=V 


# E 
Q Y N , » 
= 2: I do, +Ld. 
i—0 i=u 
Es sind nun die Koeffizienten d; so zu wählen, daß 


E 


9=rt-2t do, + N &=Minimuim ..... 0... (da) 
i=0 i—V 


wird.. Die Bedingungen dafür lauten 


do ' 


a a de = ia ale an 
dd; 2t-c,+24=0, oder d;=t-e;, für i=0,1.... 


Insbesondere ist d,—=0, weil stets c,,—=0 ist. Damit wird die bestmögliche Ausgleichs- 
funktion nach dem Ansatz (41) 


be) 


F (se) ——— b + t = Cjı Pi (x), 
i—1 
für das zugehörige @ aber ergibt sich 
o=r[i I a,)= nee era ni See 
i—1 


und daher ist die mittlere, quadratische Abweichung in Ordinatenrichtung 
Be ee 
Die Koeffizienten der Ausgleichsfunktion ergeben ausgerechnet 
d-t.0, tl (@)vi@y)w (a, y)dady 
= \\p (le) y—-b)w(a,y)dydz 
— \9 (2) (ya) -b)w, (a) dr —\p;(r)-yla)m,(a)de. 


Die Größen d; sind also gleichzeitig auch die Entwicklungskoeffizienten der Regressions- 
linie y (x), d.h. die beste Ausgleichsfunktion F (x) ist diese Regressionslinie. 


Wir erhalten damit das folgende Ergebnis: Ist eine Reihe von Beobachtungspunkten 
durch eine Funktion F(x) auszugleichen, dann wird das mittlere Abweichungsquadrat Q zu 
einem Minimum, wenn man als Ausgleichsfunktion die Regressionslinie y(x) für die Ver- 
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teilung w(x,y) der Beobachtungspunkte wählt. Der Wert von @ wird dabei gleich der 
Abweichung des Korrelationsverhältnisses k,, von Eins, noch multipliziert mit dem Streuungs- 
quadrat #? der Ordinaten. Nur im Falle vollständiger stochastischer Abhängigkeit wird mit 
kyx=1 das Abweichungsquadrat Q@=0. Die Aufgabe reduziert sich dann auf ein Inter- 
polationsproblem. ’ 

Wird die beschriebene Aufgabe der Ausgleichsrechnung unter Vertauschung der Ab- 
szissen und Ordinaten durchgeführt, dann tritt an die Stelle von 9 (x) die andere Regressions- 
linie x (y) und an die Stelle von k,, das zweite Korrelationsverhältnis k%,. Beide sind im 
allgemeinen nicht identisch, und daher ergeben sich voneinander verschiedene Ausgleichs- 
kurven und im allgemeinen verschiedene @. Ebenso wird die Ausgleichskurve etwas ver- 
ändert, wenn man vor Durchführung der Ausgleichsrechnung in Abszissen- oder Ordinaten- 
richtung Verzerrungen vornimmt, indem man z. B. die Beobachtungspunkte logarithmisch auf- 
trägt. Die Aufgabe der Ausgleichsrechnung ist also nicht eindeutig. Eindeutigkeit könnte 
sich grundsätzlich nur erreichen lassen, wenn man statt mit der einen Funktion F(x) mit 
zwei Verzerrungsfunktionen f(x) und 4 (y) arbeiten würde, wobei an die Stelle des Variations- 
problems (47), welches zum Korrelationsverhältnis führt, das allgemeine Variationsproblem 
für K? treten würde. Jedoch dürfte sich dieser Mehraufwand in der Praxis kaum lohnen. 

Die Praxis stellt aber für die Aufgabe der Ausgleichsrechnung neben der Bedingung (41) 
noch eine zweite, nicht streng faßbare Forderung, daß nämlich die Ausgleichskurve ein mög- 
lichst einfacher und glatter Kurvenzug sein soll. Die Regressionslinie erfüllt diese Bedingung 
nicht immer, vor allem dann nicht, wenn die Zahl der Beobachtungspunkte verhältnismäßig 
gering ist. Am besten eignen sich meistens als Ausgleichskurven Polynome nicht zu hohen 
Grades. Solche sind die Abschnitte der Reihenentwicklung (42), denn 

n 
ATES EN WATRE " © 
i—1 


ist ein Polynom »-ten Grades. Es hat überdies den Vorzug, daß man den zugehörigen Wert 
/ n 
) = e(i ee | REED ET 


sofort angeben kann, und daß man nur weitere Entwicklungsglieder hinzuzurechnen hat, 
wenn die Genauigkeit durch Erhöhung des Grades n verbessert werden soll. 

Die Orthogonalpolynome 9,„(x) zur Belegung w, (x) wären somit ein ausgezeichnetes 
Hilfsmittel auch für die Aufgabe der Ausgleichsrechnung. Ihre Anwendung begegnet noch 
Schwierigkeiten, solange für die Arbeit mit diesen allgemeinen Polynomen die notwendigen 
Rechenverfahren nicht durchgebildet sind. Es ist hierzu bemerkenswert, daß für den Fall 
einer Belegung w, (x), bei welcher m äquidistante Stellen alle mit dem gleichen Gewicht 1/m 
versehen sind, diese Rechnung mit den ÖOrthogonalpolynomen zum Zwecke des Kurvenaus- 
gleichs in der Praxis Eingang gefunden hat, und zwar bei den Arbeiten des Instituts für 
Konjunkturforschung. Es handelt sich dabei um die Berechnung des sog. Trends, d.h. des 
zeitlichen Verlaufs bestimmter Wirtschaftsgrößen. wenn von kurzdauernden Schwankungen 
abgesehen wird. Das Rechenverfahren wurde von Dr. Paul Lorenz im Auftrag von Prof. 
Dr. Ernst Wagemann entwickelt und istim Sonderheft 9 der Vierteljahrshefte für Konjunktur- 
forschung, Berlin (1928) mitgeteilt und in neuer Bearbeitung in Sonderheft 21, Berlin (1931). 

Als Beispiel betrachten wir noch den einfachsten Fall des Kurvenausgleichs mit Hilfe 
einer Geraden. Hier wird die Entwicklung (42a) mit dem linearen Glied abgebrochen. Es 
ist wegen c,,=r 


t . 
F(i&)=b-+tec,,9,(&2)= b+— (2 —a). 
Die Ausgleichsgerade hat die Gleichung 





tr 
: (2 — a); 


y—-b= 


’ 2 = ir 
sie geht also unter der Richtung a gemäß der Formel tg a = - durch den Schwerpunkt z= a, 


b 


y=b der Verteilung. Die mittlere quadratische Abweichung ist 


Ay=yl-r*. 
Beim Ausgleich durch eine Gerade tritt also an die Stelle des Korrelationsverhältnisses k},, 
der Korrelationskoeffizient r. 304 
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KLEINE MITTEILUNGEN 


Allseitigkeit des Druckes. Die folgenden 
Zeilen beziehen sich auf den Druck und die Ber- 
noullische Gleichung. Die Allseitigkeit des Druckes 


wird gewöhnlich durch eine Gleichgewichtsbetrach- 
tung an einem infinitesimalen Tetraeder nach- 
gewiesen, die Bernoullische Gleichung. in welche 
u. a. dieser allseitige Druck eingeht. auf dem Weg 
über die Eulersche Bewegungsgleichung äufgestellt. 
Eine elementare Darstellung. welche direkt sowohl 
die Bernoullische Gleichung wie die Allseitigkeit 
des Druckes erbringt. mag daher nicht unerwünscht 


sein. Wir betrachten der Reihe nach die ruhende 
und bewegte „ideale Flüssigkeit“ (Reibungswir- 


kung vernachlässigt); o sei die Dichte, g die Erd- 
beschleunigung. 0:9 =y das spezifische Gewicht. 
und ce die Geschwindigkeit (der Größe nach). 

x 





Bild 1. 


} A302.1 


I. Ruhende Flüssigkeit. Ausgehend von 
der Vorstellung, daß die Flüssigkeit außer der 
Schwerkraft nur senkrechten Druckkräften unter- 
liert. wollen wir an Hand von Bild 1 zeigen. daß 
der Druck auf ein Flächenelement df» in einem 
beliebigen Punkte P» unabhängig von der speziel- 
len Lage des Elementes ist. Indem wir also vor- 
läufie noch zulassen. daß der Druck von der Lage 
abhängt, sei 92 der Druck für eine bestimmte Lage 
von dfs in Ps. und analog pı der Druck für eine 
bestimmte Lage eines Elementes dfi in dem von 
Ps verschiedenen Punkt P.. Alsdann konstruieren 
wir (was auf mannigfache Weise möglich ist) in 
der Figur einen .„Stromfaden“, welcher durch die 
Berandung von d/> und dfı zeht. und der im übri 
gen geschlossen sei, wobei in der Figur von dieser 
„Ergänzung“ nur ein Teil gestrichelt angedeutet 
ist. Nehmen wir nun als Verrückung jene vor, bei 
der nur die Flüssigkeit innerhalb des geschlossenen 





Fadens etwas verschoben wird. wodurch derse!be 
in Flüssirkeitsteilchen vom zwleichen Inhalt dr 


aufgeteilt wird. derart. daß bei der Verrückung 
jedes Teilchen an die Stelle seines Vorgängers 


tritt und der Faden von der ausgezozenen in die 
strichpunktierte Lage übergeht. und bedenken. daß 
zu der Verrückunz offensichtlich auch die ent- 
verengesetzte existiert. so muß, da die Flüssigkeit 
sich im Gleichgewicht befindet. die 

{Gesamtarbeit der beteilizten Kräfte} = 0 
sein. Dieselbe setzt sich zusammen 

a) aus der Arbeit der Druckkräfte ı dv — pdv: 

b)aus der Arbeit der Schwerkraft I ydröh 

‚dv Nöh. 

Schreiben wir jetzt (1) an. heben d®v weg, und 
machen noch den Grenzüberganz drv>0. wodurch 
aus Döh das entsprechende Integral idh=kh— ha 
wird. so erhalten wir 


#3) 


(M—p)+rih, 


a (9) 
en te 
oder 


„.t+rhzntrk: - - 


Ersetzen wir in der ganzen Überlegung df.» durch 
den (in der Figur nicht gezeichneten) Querschnitt d /! 
mit dem Druck »%. so liefert uns derselbe Strom- 
faden die zu (3) analoge Gleichung 


DIESES TR sea A, 


woraus durch Vergleich mit (3) folgt 
0 SE a :  % 


von der 
> 


aus, daß der Druck unabhängig 
Schnittrichtung oder allseitig ist. Sodann «ibt (3) 
die Verteilunz des Druckes im Schwerefeld an). 

2. Bewegte Flüssigkeit. Ist die Flüssig- 
keit bewegt, so gehen wir ganz ähnlich wie vorhin 
vor. Der Unterschied ist nur der, daß der frühere 
Stromfaden jetzt der eindeutig zu der festen Lage 
von d/» gehörende Stromfaden ist. von dem dfı 
irgendeinem festen Querschnitt in einem von Ps» 
verschiedenem Punkt Pı bedeutet. und die gestri- 
chelten Fortsetzungen sich i. a. nicht schließen. 
Als Verrückung betrachten wir die durch die Strö- 


(D) sagt 


munz in dem kleinen Zeitintervall 0... dt voll- 
zogene, so daß jetzt de=dYV-df wird. wo dV 
das kleine Durchflußvolumen des Fadens ist. Und 


schließlich tritt an Stelle der früheren Gleich- 
zewichtsbedingung (1) die Arbeitsgleichung 

fGesamtarbeit der\ __fZunahme der ki- 6 

\beteilieten Kräftef ° \netischen Ensegio} ai 
Alsdann geht. wie man sieht. unsere neue Über- 
legung in dem Spezialfall, daß die Flüssigkeit ruht. 
vollständig in die frühere über, und es wird daher 
die neue aus der Gleichung (6) sich ergebende Be- 
dingung die frühere (3) als Spezialfall in sich ent- 
halten. Die linke Seite von (6) berechnet sich wie 
oben im Falle der ruhenden Flüssirkeit. Was die 
rechte Seite anbelangt. so führen wir für die Ge 
schwindigkeiten in Pı und Ps die Bezeichnungen 
Cı. C2 ein und deuten. wo dies nötig ist. durch den 
oberen Zeiger (0) bzw. (d?) noch den Zeitpunki 
an. auf den sich die jeweilige Geschwindigkeit be 
zieht, so daß also z. B. cs(dt) die Geschwindierkeit c 
in Ps zur Zeit d? bedeutet. 

Bilden wir nun die Differenz der Bewegrungs- 
energie des strichpunktierten Fadens zur Zeit d/ 
und des auszezorenen Fadens zur Zeit o. so setzt 
sich diese aus den beiden folgenden Gliedern zu- 
sammen: 


1. aus einem Glied. welches herrührt von dem 
beiden Fäden nicht gemeinsamen Teil: 
er 2 r lo 
= odr ([e,14t | — [c, >) 
2. aus einein Glied. welches herrührt von den 
beiden Fäden gemeinsamen Teil und das wir an 


der durch die Verrückung in der Fizur gerebenen 
Einteilung in Raumzellen berechnen: pro Zelle be- 
trägt dann der Unterschied 


1 odv sfetd n] 2. [e®]2\ : 


hier ist die geschweifte Klammer offenbar die zeit- 
liche Änderung dc? (für die betreffende Zelle) oder 


dc? c 
c c . 
de cedt und also wegen cdf=Öös mit- 
of Öt j 

. „Öc : 2 

hin 2 ös, so daß man insgesamt erhält 


. 


\ "dc 
C 
odı e. 34 08- 


weg und 


d. h. 


1) Wie ieh nachträglich bemerkte, findet sieh dieser Teil 
der Überlegung im wesentlichen schon bei A. Föppl, 
Vorlesungen über Teehnische Mechanik, i. Bd. 1938, S. 
371358, 


Schreiben wir jetzt (6) an. heben dı 
machen noch den Grenzübergang dt, 
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x v N s 
. s oc 
speziell auch do—>0. so wird aus \ ös das 


Öt 
"de i 
entsprechende Integral 1 s, und wir erhalten 
c 
P >. 
0 ‚m 2. | "de » (7) 
9, —- Pe.) Tr h,—h, =. ıC: — CH TO ri 
P, 
oder 
P» 
ER 
Me R Be ( 
Ptrk za=mntikhryan oe [5,45 ®. 


« 
> 
1 


wo nunmehr alle Größen sich auf den Zeitpunkt o 
beziehen. Ersetzen wir in der ganzen Überlegung 
df, durch den (in der Figur nicht gezeichneten) Quer- 
schnitt d73 mit dem Druck 2%, s0 liefert uns der- 
selbe Gedankengang die zu (8) analoge Smichung 


0 © oc 
” SEE u ee > (( 
mn tr > ci=p ya ( o| 5,48 (9). 
Pı 
woraus durch Vergleich mit (8) folgt 
1 


(10) sagt aus, daß der Druck unabhängig von der 
Schnitttläche oder allseitig ist, wobei dieser Druck 
stets auf das mitbewegte Flächenelement zu be- 
ziehen ist. Sodann gibt (9) die allgemeine Bernoul 
lische Gleichung für instationäre Strömungen. 

3. Stellt man sich übrigens den Druck als eine 
Folge der mehr oder weniger engen Packung der 
Moleküle vor. so ist seine Allseitigkeit evident. 
Ist die Flüssigkeit reibend, so führt man den 
Mittelwert der Normalspannung als Druck ein, der 
dann definitionsgemäß allseitig ist. Im Falle der 
tuhe gilt auch jetzt wieder (3), während für be- 
wegte Flüssigkeiten auf der rechten Seite von (8) 
ein .„Verlustglied“ von der Dimension eines 
Druckes tritt. welches daher rührt, daß mecha- 
nische Energie in Wärme übergegangen ist. 

4. Spezialisieren wir auf den Fall. daß die Strö- 
mung aus der Ruhe heraus erfolgt (unter Einfluß 
der Schwere), d. h. die Stromfäden zu jeder Zeit 
aus einem Bereich kommen, wo die Flüssigkeit 
praktisch ruht. so erhalten wir. wenn wir den 
Punkt Pı jeweils in diesem Bereich wählen. aus (8) 

P» 


m 
pn t7h=m+trh,+5c’+o | 5,45 Kar) 
Pı 
Hier ist nun aber die linke Seite nach 1. für alle 
Stromlinien dieselbe. d.. h. eine nur noch von der 
Zeit abhängende Konstante ( (f). Solche Strömun- 
zen sind nun bekanntlich drehungsfrei. d. 08 ist 
> 


> 


(zu einer festen Zeit) für jede Stromlinie \ cds 
P 


1 
einfach das Geschwindigkeitspotential ® im Punkte 
Pz und daher das in (11) auftretende Integral 
dessen zeitliche Änderung 

Ps Ps 


de En _d» = 
[3 yfeıs- ee Ze |. 


P, Pı 
womit wir für solche Strömungen die Gleichung 
erhalten (wenn jetzt der Zeiger 2 überall weg- 
gelassen wird) 
o RX 
+ yh+-—c?-+o _ 
pP gro 
Die letzte Beziehung gilt nun für alle Potential- 
strömungen, da, wie bekannt, jede Potentialströ- 
mung durch plötzliche Stoßdrucke auf eine die 


BE 46 U 


Flüssigkeit umgebende Membran, deren Inneres 
zuächst ruht. erzeugt werden kann. 
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Zur Berechnung des kleinsten Eigenwerts 
von y"—+ip()y=d. 

Es sei gegeben y„”+Ap(z)y=0 mit den Rand- 
werten (0) =y(1)=0. Im Intervall UuS7SI1 
sei » positiv und endlich. 

Zur Berechnung des ersten Eigenwerts 4, ist in 
der Praxis folgendes Verfahren seit langem üblich. 

Man geht von einer Funktion yı aus, welche die 
Randbedingungen erfüllt und sich dem vermuteten 
Verlauf der ersten Eigenfunktion möglichst an- 
nähert: also keine Nullstelle im Intervall O<r<1 
hat. 

Nach dem Schema 

Yntı: Pundadz.... W 
leitet man durch zeichnerische oder numerische 
Integration aus der „Ausgangsfunktion „1“ eine 
Folge von Funktionen %». %s. . Yn+ı ab, wobei 
die Integrationskonstanten bei den Integrationen 
so bestimmt werden, daß die Funktionen dieser 
l’olge die Randbedingungen erfüllen. 

Setzt man Q (r) — Y„/Y„+, und ist Qı der größte, 
Q>: der kleinste Wert, den Q annimmt, so ist 
Qı >4ı >Qe und bei hinreichend großem rn wird 
Qı — Q> beliebig klein. Der erste Eigenwert läßt 
sich mithin nach diesem Verfahren mit beliebiger 
(Genauigkeit berechnen. Ein kurzer Nachweis für 
die Konvergenz dieser .„Iterationsmethode“ wird 
nachstehend gegeben. 

Es seien 9; die Eigenfunktionen und 4, die Eigen 
werte der Randwertaufgabe, wobei O<A, <<... 
sein soll. 

Man denke sich nun die Ausgangsfunktion yı 
nach den Eigenfunktionen entwickelt 


Yı = ı Ta TCGPT a; R 


Da jede Eigenfunktion der Gleichung 9%’ +4;,p9; 
0 und den Randbedingungen zenügt. ist auch 


P; )4,|Spy;drdr und es wird 
Cı Tı | CaP2 , Cs Ps 
0-1 i.n- r 
ai Yn _ A a1 el Q 
) = BR = (oO) 
Yn+ı Cı Yı , CaP2 |, Cs Pa | 
an Te Ton 


wenn man die Integration (1) mit der Entwicke 
lung (2) durehführt. Multipliziert man den Bruch 
rechts mit /% im Zähler und Nenner, so wird 


’ 3, \n-ı 7, \nmi 
amtaWel; ) +69; : 
Q=/ÄÜ. > > 
? k i .\n i.\n (4). 
: IR hy B hy N 
ayıt Pal; a 3 Y3 7 
| ke} 


Da die Klammerpotenzen mit wachsendem n be- 
liebig klein werden, folgt aus (4), daß @ bei hin- 
reichend großem n beliebig wenig von 4, abweicht. 

Es bleibt noch zu zeigen, daß Q@ Werte annimmt, 
die oberhalb und unterhalb 4, liegen. Bezeichnet 
man mit B den Faktor von A, in (4), so genügt es, zu 
zeigen, daß die Funktion (B—1) im Intervall 
0<r<1 das Vorzeichen wechselt. Es ist 
B—-i= 

[ A, n 1 hı n 3 n 
9; (2 tanz) E 

(A, n A, \" 
HyMTrT 2r;) I rl; 
Der Nenner dieses Ausdrucks hat nach den über 
Yı gemachten Voraussetzunzen keine Nullstelle. 
Dagegen wechselt der Zähler Z sein Vorzeichen. 
Wegen der Orthogonalität der Eigenfunktionen ist 
nämlich 

1 

\Zpgy,de=0. 

a 
Da pY, sein Vorzeichen nicht wechselt, muß dies 
bei Z der Fall sein. Damit ist bewiesen, daß 
erstens Qı >4, > Qs ist. und zweitens, daß diese 
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}, eingerenzenden Werte bei hinreichend großem n 
beliebig nahe zusammenrücken. 

Aus den vorstehenden Entwickelungen ist auch 
ersichtlich, wie die Wahl der Ausgangsfunktion 
yı die Konvergenz beeinflußt. Je besser yı ge- 


troffen wird, um so kleiner werden die Beiwerte 
Ca, (3, in der Entwickelung (2). Andererseits 
erkennt man auch, daß diese Methode nur für den 
ersten Eigenwert brauchbar ist. 


Hannover. H. v. Sanden. 311 


BUCHBESPRECHUNGEN 


Dr. WALTHER LIETZMANN, Frühgeschichte 


der Geometrie auf germanischem 
Poden. 94 S. m. 91 Abb. Breslau 1940. Verlag 


Ferdinand Hirt: Preis geb. 3.50 M. 

Das .Buch macht den Versuch, „die Entstehung 
geometrischer Vorstellungen aus der Freude an 
schmückender Formgestaltung“ in der germani- 


schen Frühzeit nachzuweisen; die dafür herange- 
zogenen Funde entstammen dem Bereiche von 


Norddeutschland, Dänemark, Südschweden und 
Südnorwegen. Verf. betont aber selbst (S. 15), dab 
„der hier behandelte große Kulturkreis der Vor- 
geschichte zur Entwicklung unserer wissenschaft- 
lichen Geometrie im Gegensatz zu der babylo- 
nischen, ägyptischen, griechischen Frühgeometrie 
kaum Wesentliches beigetragen hat“. In der Tat 
eeht aus den Funden wohl kaum mehr hervor als 
ein ästhetisches Gefühl für die Schönheit gewisser 
Linienführungen und Formen. Weitergehende 
Schlüsse auf das Vorhandensein geometrischer Be- 
eriffe oder gar auf die Kenntnis geometrischer 
Sätze zu ziehen. wäre wohl allzu kühn und ge- 
sucht. Auch Lietzmann selbst drückt sich in dieser 
Beziehung nur sehr vorsichtig oder mit Frage- 
zeichen aus. Man kann daher wohl auch weniger 
von einer „Frühgeschichte als vielmehr von einer 
„Vorgeschichte“ der Geometrie sprechen, die ohne 
weitere Folgen für die Wissenschaft geblieben ist, 
wenn sie auch vom kunstgeschichtlichen Stand- 
punkte aus, namentlich für die Geschichte der 
Ornamentik, Interesse bieten mag. 

Berlin. E. Mosch. 275 

Dr.-Ing. BRUNO ECK, Technische Strö- 
mungslehre. VIL+ 252 S. m. 278 Abb. Ber- 
lin 1941, Verlag Julius Springer. Preis geb. 13,50 M. 

An Stelle der 1936 erschienenen „Einführung in 
die Technische Strömungslehre“, die in 2 Bände, 
„Theoretische Grundlage“ und „Strömungstechni- 
sches Praktikum“, geteilt war, läßt der Verfasser 
nunmehr ein Buch „Technische Strömungslehre“ 
erscheinen, dem ein weiteres über die Versuchs- 
und Laboratoriumstechnik der Strömungslehre 
folgen soll. 

Das vorliegende Buch ist eine Erweiterung des 
ersten Bandes des früheren Werkes). Diese Er- 
weiterung besteht zum Teil in einer breiteren Dar- 
stellung des Inhaltes des früheren Werkes, zum 
Teil in der Hinzufügung neuer Kapitel (z. B. über 
Kavitation und Gasdynamik). Besonders stark er- 
weitert sind die Ausführungen über die Wirkung 
der Reibung in strömenden Flüssigkeiten und alles 
was damit zusammenhängt: Rohrströmung, Kör- 
perwiderstand, Strömungsablösung, Maßnahme zur 
Widerstandsverminderung usw. Hier findet man 
eine gute Darstellung des heute gesicherten Stan- 
des der Erfahrung und Theorie, während die zur 
Zeit auch schon behandelten schwierigeren Pro- 
bleme der Turbulenzforschung zweckmäßigerweise 
noch nicht erwähnt sind. Sehr schöne Bilder, größ- 
tenteils vom Verfasser selbst aufgenommen, und 
die Darlegung vieler Versuchsergebnisse unter- 
stützen und ergänzen hier den verständlich und 
anschaulich geschriebenen Text: auch ist genügen- 
des Zahlenmaterial zur unmittelbaren praktischen 


1) Siehe die Besprechung in Z. angew. Math. Mech. 


Bd. 16 (1936) S. 188. 


Verwertung durch den Ingenieur beigefügt. Der 
Tragflügelströmung ist eine recht hübsche, durch 
markante Bilder unterstützte Darstellung gewid- 
met. Das Kapitel über Gasdynamik bringt das 
Wesentlichste aus diesem Gebiet in gut verständ- 
lichen Ausführungen. Im Schlußkapitel findet man 
nützliche Bemerkungen über strömungstechnische 
Messungen und Beschreibungen von Meßapparaten. 

All die bisher erwähnten Kapitel geben einen 
schönen Überblick über die Probleme und bespre- 
chen auch in anregender Weise eine Reihe tech- 
nischer Anwendungen. 

Dagegen sind gegen die Darstellung der Haupt- 
gesetze der Strömungslehre, die dem oben er- 
wähnten Kapitel vorausgeht, erhebliche Bedenken 
zu erheben. Jeder, der auf diesem Gebiet lehrend 
tätig ist, wird die Erfahrung machen, daß gerade 
bei einer Einführung in diesen Gegenstand 
äußerste Genauigkeit und Schärfe in Ausdruck 
und Formulierung und im Abwägen des Wesent- 
lichen und Unwesentlichen notwendig ist. Mit Ver- 
einfachungen und Vernachlässigungen ist häufig 
gerade dem Anfänger nicht gedient, insbesondere 
wenn sie zu Scheinbeweisen und ausgesprochenen 
Unrichtigkeiten führen. Er wird dadurch nicht 
„schonend“, wie sich der Verfasser im Vorwort 
ausdrückt, auf weitergehende Darstellungen vor- 
bereitet. 

Durch eine entsprechende Umarbeitung dieses 
Teiles würde das Buch nur gewinnen. 

Karlsruhe. Spannhake VDI. 318 


Dr.-Ing. habil. FRANZ KARAS VDI, Elasti- 
sche Formänderung und Lastvertei- 
lung beim Doppeleingriff gerader 
Stirnradzähne (VDI-Forschungsheft 406, Bei- 
lage zu „Forschung auf dem Gebiete des Ingenieur- 
wesens“, Ausgabe B, Bd. 12, Januar/Februar 1941). 
23 S. m. 13 Bildern u. 4 Zahlentaf. Berlin 1941, 
VDI-Verlag G.m.b.H. Preis brosch. 5 M. 

Bei den üblichen Zahnrädern sind abwechselnd 
ein Zahn, dann zwei Zähne im Eingriff. Von Be- 
lang ist die Kenntnis der Formänderung, um hier- 
aus die Lastverteilune und die Profilkorrektion für 
den stoßfreien Eingriffsbeginn eines neuen Zahnes 
zu bestimmen. Dr. Karas berechnet die Form- 
änderung unter Annahme eines starren Zahnrad- 
körpers aus den drei Anteilen durch Biegung, 
Schub und Walzenpressung, die sich von gleicher 
Größenordnung ergeben und die Größe der Her- 
stellungsfehler erreichen. 


Dresden. GC. Weber. 296 


Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 


bleibt vorbehalten): 


Dr.-Ing. WILHELM CAUER, B. Apl. Prof. a. d. 
Techn. Hochschule Berlin u. Leiter des Labora- 
toriums der Mix & Genest AG.. Theorie der 
linearen Wechselstromschaltungen, 


1. Bd. XIl+614 S. m. 426 Abb. Leipzig 1941, 
Akademische Verlagsgesellschaft Becker & Erler 


Kom.-Ges. Preis geb. 45 M. 


VIII + 300 S. 
Oldenbourg. 


HEINRICH DÖRRIE, Vektoren. 
München u. Berlin 1941, Verlag R. 
Preis geb. 13,50 M. 
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- NACHRICHTEN 


Professor Dr. Kurt Beyer 60 Jahre alt. 

Am 27. Dezember 1941 vollendete Prof. Dr. Ing. 
Kurt Beyer das 60. Lebensjahr. Beyer ist Mit- 
glied des Beirates der Zeitschrift für Angewandte 
Mathematik und Mechanik und gehört zu den 
Gründern der dieser Zeitschrift eng verbundenen 


Gesellschaft für Angewandte Mathematik und 
Mechanik. 
Nach mehrjähriger Assistententätigkeit an der 


T. H. Dresden, wo er auch studiert und promoviert 
hatte, ginz er 1908 als Sektionsingenieur der 
Staatsbahnen nach Siam, wo unter seiner Leitung 
zahlreiche Brücken und andere Bauten entstanden. 
Kurz vor Beginn des Weltkrieges kehrte er nach 
Deutschland zurück und trat 1915 freiwillig in das 
Heer ein. Nach Verwendung in verschiedenen an- 
deren Stellen war er zuletzt beim deutschen Feld- 
eisenbahnchef in Kleinasien tätig. 1919 wurde Beyer 
dann die ord. Professur für Statik der Baukon- 
struktionen und technische Mechanik an der T. H. 
Dresden übertragen. Auch weiterhin hielt er die 
enge Verbindung mit der Praxis aufrecht. Als 
technischer Berater großer Firmen war er führend 
insbesondere an der Entwicklung der Förderan 


lagen des Braunkohlentagebaus beteiligt. Groß- 
bagger und Abraumförderbrücken größten Aus- 


maßes entstanden unter seiner Mitwirkung. 

Seine reichen Erfahrungen und seine Unter- 
suchungen legte er in einer Reihe wissenschaft- 
licher Veröffentlichungen nieder, von denen das 
-weibändige Handbuch „Die Statik des Eisenbeton- 
baus“ (1933/34) besonders hervorgehoben werden 
muß. Für die Zeitschrift für Angewandte Mathe- 
matik und Mechanik schrieb er zahlreiche wertvolle 
Referate von Büchern seines Fachgebietes. 


ZUSCHRIFTEN AN 


H. Okubo: The stress distribution in a semi- 
infinite domain having a plane boundary and 
compressed by a rigid body (Z.angew. Math. Mech. 
Bd. 20 (1940) S. 271» 276). I received from Pro- 


Neben Praxis und Forschung entfaltete Beyer an 
der T. H. Dresden eine umfangreiche Lehrtätigkeit. 
Zahlreiche tüchtige, von den Baufirmen gesuchte 
Ingenieure gingen aus seiner Schule hervor. Mit 
E. Trefftz, ©. Weber und anderen Dresdner Hoch- 
schullehrern hält er seit Jahren ein Seminar über 
Festigkeitslehre, in dem manchem seiner Schüleı 
die Anregung zur Dissertation gegeben wurde. 

Zu Beginn dieses Krieges stellte sich Beyer er- 
neut der Wehrmacht zur Verfügung. Er machte 
den Polenfeldzug mit, wurde aber Anfang 1940 zur 
Betreuung wehrwirtschaftlich wichtiger Aufgaben 
in die Heimat zurückgerufen. 

Dem auf der Höhe seines Wirkens stehenden 
Mann der Praxis, dem Forscher, dem Lehrer, seien 
noch viele Jahre erfolgreicher Arbeit beschieden. 
Möge er weiter wie bisher dieser Zeitschrift ver 
hunden bleiben. Willers. 333 


Neu berechnete Tafeln. 
F. Tricomi, Atti R. Ace Se. Torino, vol. 76 (1941 
gibt bei n=1 2... 10 die Werte der Laguerre 
schen Funktionen 


t 
e’ an 
BE 45 ni dem (E 


t gm, 

mit 4 Dezimalen für t—=01 02... 09 und mit 
5 Dezimalen für t = 1,00 1.25... 2,75: 3.0 35... 5,5: 
6.T.., 5 U Bi 


Stuttgart. Fritz Emde. 335 
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DEN HERAUSGEBER 


P. 273: In the expression of 3 e the 3rd. term 
s 


an 
2 en D„e®@"rtD« is added. Egs. (13) and (14) 
n=0 
are replaced 


with (+2u) Bun — Du=0 and 











fessor Dr. Göhner some remarks on ıny paper N SEE et a Au 
published in this magazine and will make the And ee umpacksety. E. DUB: " 2 vn 
corrections as follows: term nn C„ is replaced with — - ; and 
Table I. 

n An Bı Cy D, 

v - 0.740 X 10 P — 1,856 X 105 P 

1 0.507 X 10° P 0.069 — 0,416 X 105 P 0.173 

2 0.055 0.025 — 0.045 0.063 

3 0.016 0.013 - 0.013 0.033 

4 0.006 0.007 0.005 0.019 

5 0.004 0.005 0.003 0.013 

6 0.001 0.003 - 0.001 0.008 

7 0.001 0,002 0.001 0.006 

8 0.000 0.002 - 0.000 0.004 

9 — 0.001 0.001 0.001 0.002 

10 0.001 0.001 - 0.001 0.002 

11 -— 0.000 0.000 
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Stresses in unit P 





Fig. 4. 


Ay with — And the numerical calcula- 


u 
A+2u' 
tions in $ 6 are corrected as follows: 

the simultaneous eqs. to determine the constants 


An, Bn,Cn and D„ are A=— 1210 X 10°C, 
Bu = 0.3989 X 10° D,„, 
nP 
ER ) - >% u Wi 
On = 03774 %X 10° 1,3 
Din 


— 9547 \' 
2.547n E@m+1?-4®' 


y' 1.783 n — 0.8333 m —04 17 
eg (2m+1?—4An? 


m 


e Din 


nP 


= 0.1441 X 10% 5, P+ 1.153 X 10 +; 


where 


ml, 1, 2, he Re 11. 

The values of the constants calculated from the 
above simultaneous egs. are tabulated in Table l. 
In this case the vertical displacement of the con- 
tact surface is U, = — 0,516 X 10° P. To compare 
the displacement of the contact surface in the case 
treated in $ 4 with the result obtained here, the 
displacements of both cases are plotted in Fig. 3. 
Next we will calculate the stresses on the contact 
surface. From eq. (8). at a=(, we have 


il u} > 
2(A+Uu U, F 
— >’ C„2ncos2nB-— - 
sind 2 
n=] 
7; 2(itu)u 
ap r 
A; +2 u) sin d 


ni 
aa 


2 
P% Du(2n+1)cos(2n—+1)2 
n=V 





The values of stresses calculated from the above 
eqs. are plotted in Fig. 4°). The second equation 
in (21) does not hold at d=0 (or d=n), so uz in 
(20) must vanish at $=0 if the numerical solution 
obtained here is continuous. While u3 is rot 
exactly zero but is + 0.051 X 10° PP. The discere- 
paney is caused by the neglection of the higher 
terms of the trigonometrie series which give some 
influence on the result near the point d=0. 


Yonezawa, Japan. H. Ökubo. 319 


zur 
all- 


Willers: Benutzung projektiver Skalen 
Unterteilung von Skalen anderer Funktionen. 7. 
zsew. Math. Mech. Bd. 20 (1940) S. 291 bis 292. 

Herr Dr. Heinrich-Breslau macht mich 
freundlicherweise auf zwei Versehen in der Durch- 
rechnung des Beispieles in obiger Notiz aufmerk- 
sam, deshalb sei die Rechnung hier nochmals kurz 
wiederholt. Die Abweichung der durch projektive 
Imterteilung zwischen den drei Punkten 


4. sT OT 

r, zu ei I. = Ri I — N 

18 BET 18 

einer te-Skala 
exakten wird 


gewonnenen Teilpunkte von den 
nach Schwerdt abgeschätzt durch 
Ss 2 a-eı) 3 
j} 2.3 E P% | £i< R 
RR u 4 (+ a) | 36) 
hier R < 0,001. 
Formel eibt 


Das eibt 
zeleitete 


Die in der Mitteilung ab 


12sin£ RE: 6 
18 (e 1.785) + : 


3, \> 


[6--4 cos? 


cos? £ cos? &) 


COS 


Die eckige Klammer ist in dem Bereich von 40° bis 
zu 50° positiv und kleiner als 4,2, so daß hier 
R < 0,0002 wird. Soll die Abweichung der ge- 
näherten Teilung von der exakten kleiner als 
0,01 em werden, würde man nach der ersten Ab- 
schätzung Maßeinheiten bis zu 10 em, nach der 
zweiten bis zu 50 cm als zulässig erhalten. An dem 
wesentlichen Ergebnis, daß die zweite Formel zu 
einer brauchbareren Abschätzung führt als die erste, 
ändert sich also nichts. 


Dresden. Willers. 354 


H. Görtler, Instabilität laminarer Grenzschichten 
an konkaven Wänden gegenüber gewissen drei- 
dimensionalen Störungen. Z. angew. Math. Mech. 
Bd. 21 (1941), S. 250 bis 252. 

Auf S. 250 ist ein Druckfehler übersehen worden. 
In (2) muß die letzte Zeile natürlich lauten: 

— 2a? ö? ( °) 2 Uuw. 


v 


Der Buchstabe U war vergessen worden. 


Göttinzen. H. Görtler. 332 
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